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Abstract 


This  is  the  third  and  final  annual  technical  report  of  a  three  year  research  program 
to  study  both  experimentally  and  theoretically  the  propagation  of  intense  charged 
particle  beams  into  vacuum.  This  program  was  initially  funded  by  the  Air  Force  Office 
of  Scientific  Research  for  the  period  April  1,  1985  to  April  30,  1987  under  Grant 
No.  AFOSR-84-0091.  The  new  three  year  program  is  a  continuation  of  that  research 
under  the  same  Grant  Number.  The  goal  of  this  research  is  the  generation  of  nearly 
charge-neutral,  current-neutral  charged  particle  beams  (plasmoids)  which  are  capable 
of  propagating  beam  energy  in  vacuum  at  a  significant  fraction  of  the  velocity  of  light. 


1  Introduction  and  Synopsis 

Studies  of  the  propagation  of  intense  charged  particle  beams  into  vacuum  have  been  pursued 
at  the  University  of  Maryland  under  AFOSR  sponsorship  since  1984.  The  principal  goal  of 
this  program  has  been  to  investigate  under  what  conditions  intense  charged  particle  beam 
energy  may  be  propagated  into  free  space  without  the  benefit  of  confining  magnetic  fields 
or  conducting  boundaries.  Initial  experimental  and  theoretical  work  has  centered  around 
the  propagation  of  intense  electron  beams  into  vacuum  after  passage  through  an  ion  source 
localized  to  the  injection  point.  In  such  systems  propagation  can  occur  when  ions  are  accel¬ 
erated  from  the  source  region  into  vacuum  by  the  space  charge  fields  of  an  electron  beam 
injected  at  a  current  level  above  the  space  charge  limit.  The  resultant  neutralization  of  the 
electron  beam  by  the  accelerated  ions  can  allow  effective  propagation  in  the  vacuum  region. 

The  research  progress  during  this  final  year  is  contained  in  Section  2.  Studies  on  beam 
propagation  with  no  applied  fields  have  involved  our  Laser  Controlled  Beamfront  Accelerator 
concept.  Both  experiment?!  <lIlU  retical  insults  are  summarized.  A  second  exp<-ri. cental 
effort  on  the  production  of  high  quality  electron  beams  via  a  “pseudospark”  device  has  been 
initiated  in  addition  to  general  studies  of  pulsed  power  systems  and  associated  diagnostics. 
The  major  achievements  are  as  follows: 

A.  In  the  laser-controlled  beamfront  accelerator  experiments,  an  ionization  channel  is  gen¬ 
erated  by  time-sequenced  formation  of  plasma  along  the  drift  tube  wall  that  contributes 
to  effective  beam  propagation  into  vacuum.  This  is  accomplished  by  sequential,  time- 
delayed  interaction  of  laser  beams  with  targets  on  the  drift-tube  wall.  The  intense 
electron  beam  injected  into  this  system  will  propagate  with  a  beamfront  velocity  Vf(t) 
that  is  determined  by  the  timing  of  the  ionization  channel  in  the  drift  tube.  Ions  are 
pulled  in  from  the  plasma  near  the  wall  to  essentially  space  charge  neutralize  the  elec¬ 
tron  beam.  Some  ions  trapped  in  the  space-charge  well,  also  known  as  the  “virtual 
cathode,”  of  the  electron  beamfront,  will  be  accelerated  to  a  velocity  vj(t).  These 
trapped  ions  are  collectively  accelerated  to  high  energies  and  aid  in  the  formation  of 
a  charge  and  current  neutralized  beam  system  for  propagation  into  free  space.  The 
key  requirement  is  that  the  ions  remain  trapped  in  the  well  or,  alternately,  that  the 
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laser-controlled  beamfront  velocity  is  synchronized  with  the  changing  accelerating  gra¬ 
dient  associated  with  the  well.  Our  first  experiment  with  a  900  keV,  20  kA,  30  ns 
electron  beam  was  based  on  a  gradient  of  40  MV/m  and  achieved  successful  proton 
acceleration  to  about  18  MeV  in  a  distance  of  40  cm.  Attempts  to  increase  this  energy 
by  increasing  the  laser-controlled  beamfront  velocity,  consistent  with  an  accelerating 
gradient  of  90  MV/m,  were  unsuccessful.  This  result  is  explained  by  numerical  studies 
that  showed  the  degradation  of  the  beamfront  electric  field.  Higher  iou  energies  can 
only  be  achieved  by  increasing  the  electron  beam  energy,  say  from  900  keV  to  about 
1.5  MeV.  A  second-generation  experiment  based  on  these  theoretical  results  has  been 
designed  and  built.  It  uses  a  gradient  of  60  MV/m  and  a  distance  of  100  cm.  Control 
of  the  beamfront  over  the  100  cm  distance  has  proved  successful,  achieving  a  0.3  c 
beamfront  velocity. 

B.  Theoretical  research  has  been  concerned  with  four  problems:  the  initial  phase  of  ion 
production  and  deration,  the  simulation  of  cur  laser  controlled  beamfront  acceler¬ 
ation  experiment,  the  propagation  of  a  partially  neutralized  electron  beam  across  a 
transverse  magnetic  field,  and  the  ion  hose  instability  of  a  partially  neutralized  beam 
with  Bennett  profile.  Each  of  these  four  activities  produced  excellent  results  that  have 
been,  or  will  be,  published.  Particularly  noteworthy  is  the  good  agreement  with  the 
laser  controlled  beamfront  acceleration  experiments,  as  already  mentioned  in  A  above. 
The  successful  modeling  of  the  rather  complicated  physics  including  the  ionization  by 
the  electron  beam  is  very  helpful  not  only  in  the  design  of  new  experiments,  but  also 
in  the  overall  evaluation  of  this  type  of  accelerator. 

C.  During  this  contract  year  we  started  a  new  experiment  called  the  "pseudospark  dis¬ 
charge.”  A  novel  device,  invented  in  Germany,  the  pseudospark  has  many  interesting 
applications  such  as  pulsed-power  switching  and  electron  beam  generation.  We  are 
particularly  interested  in  the  latter  application.  Our  first  experiments  have  demon 
strated  that  the  pseudospark  discharge  with  suitable  electrode  configuration  produces 
an  electron  beam  of  extremely  high  quality.  Measurements  showed  a  beam  of  25  kV. 
100  A,  10  ns  with  a  normalized  brightness  of  10u  A/(m-rad)2.  This  is  considerably  bet- 
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ter  than  the  quality  of  beams  from  thermionic  cathodes  and  makes  the  pseudospark 
source  very  attractive  for  advanced  accelerator  applications  (linear  colliders,  FELs, 
etc.).  The  physics  of  the  pseudospark  device  is  not  yet  well  understood,  and  we  believe 
that  we  can  make  an  important  contribution  in  this  area.  From  preliminary  scaling 
considerations,  we  expect  even  higher  current  densities  [>  103  A/cm2]  and  brightness 
values  [>  10u  A/(m-rad)2j  at  increased  voltages. 

A  complete  list  of  papers  and  presentations  that  have  resulted  from  our  research  during 
the  past  several  years  is  given  in  Appendix  A,  with  copies  of  more  recent  publications  included 
in  Appendix  B. 
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2  Research  Progress 

2.1  Experimental  Research 

2.1.1  Experiments  on  the  DRAGON  Pulse  Line  Accelerator 

A  systematic  experimental  measurement  of  currents  in  a  low  pressure,  gas  filled  drift  tube, 
into  which  a  high-current  electron  beam  is  injected,  was  carried  out.  In  order  to  produce 
high-current  electron  beams  of  good  reproducibility,  we  developed  a  new  diode  system  in 
which  an  anode  with  a  specially  designed  foil  changer  is  employed.  A  paper  detailing  this 
diode  with  foil  changer  is  published  in  the  Rev.  Sci.  Instrum.  £0,  3556  (1989),  and  enclosed 
in  Appendix  B.  We  observed  that  as  gas  pressure  in  the  drift  chamber  is  increased,  the 
magnitude  of  the  net  current  increases  up  to  a  peak  value,  and  then  decreases.  The  pressure 
at  which  this  peak  value  is  achieved  was  interpreted  as  the  critical  pressure  required  for 
full  charge  neutralization  of  the  beam.  Above  this  critical  gas  pressure,  the  waveform  of 
the  current  decays  nearly  exponentially  with  a  slow  decay  time  constant.  The  decay  time 
constant  was  attributed  to  L/R  of  the  drift  tube  and  plasma  column,  and  the  resistivity  of 
the  plasma  can  be  inferred  from  the  time  constant.  A  paper  detailing  this  study  is  published 
in  the  Proc.  1989  IEEE  Particle  Accelerator  Conference,  p.  1052  (1989). 

2.1.2  Laser  Controlled  Beamfront  Experiments 

The  basic  concept  behind  the  Laser  Controlled  Beamfront  Experiment  is  shown  in  Fig.  1. 
An  intense  relativistic  electron  beam  is  injected  into  an  evacuated  drift  tube  at  a  current 
level  several  times  the  vacuum  space  charge  limit,  given  approximately  by 

=  17,000(7o2/3-l)3/2 

1  (1  +  '2£nb/a)(l  —  /)  1  ’ 

where  b  is  the  drift  tube  radius,  a  is  the  beam  radius,  7o  is  the  relativistic  mass  ratio  for 
the  electrons  at  injection,  and  /  =  n,/ne  represents  any  charge  neutralization  provided  by 
positive  ions.  As  indicated  in  Fig.  la,  a  virtual  cathode  forms  at  the  injection  point  with 
a  depth  approximately  equal  to  the  anode-cathode  potential  difference.  The  axial  position 
of  the  virtual  cathode  downstream  of  the  anode  plane  is  usually  on  the  order  of  the  anode 
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cathode  gap,  so  that  megavolt  potentials  at  the  virtual  cathode  are  formed  only  millimeters 
away  from  the  grounded  anode.  It  is  this  very  high  electric  field  on  the  upstream  side  of 
the  virtual  cathode  that  collective  accelerators  usually  seek  to  exploit  for  ion  trapping  and 
acceleration. 

If  the  beam  is  injected  into  a  localized  gas  cloud  (Fig.  lb),  ionization  processes  can 
quickly  build  up  sufficient  ion  density  to  neutralize  the  electron  beam  space  charge  and  the 
virtual  cathode  can  move  downstream  to  the  edge  of  the  gas  cloud.  This  motion  can  result  in 
the  acceleration  of  a  few  ions  to  energies  considerably  higher  than  the  depth  of  the  potential 
well  at  the  virtual  cathode. 

In  order  to  control  the  motion  of  the  virtual  cathode  over  distances  greater  than  a  few 
centimeters,  however,  a  means  of  providing  an  ionization  channel  whose  axial  extent  can  be 
controlled  as  a  function  of  time  is  required  (Fig.  lc).  In  the  present  experiments,  shown 
schematically  in  Fig.  2,  the  ionization  channel  is  generated  by  time-sequenced  laser-target 
interactions.  A  laser  pulse  is  separated  into  many  approximately  equal  energy  beams  which 
are  then  optically  delayed  over  different  path  lengths.  The  laser  light  then  vaporizes  and 
ionizes  a  target  material  on  the  drift  tube  wall,  and  ions  drawn  into  the  beam  by  the  electron 
space  charge  at  the  beamfront  provide  the  required  time-sequenced  channel  of  ionization  to 
control  beamfront  motion. 

During  the  last  year,  a  complete  1  meter  Laser  Controlled  Beamfront  Experiment  has 
been  designed,  constructed,  and  tested  (see  Figs.  2  and  3).  This  new  experiment  is  a  signif¬ 
icant  upgrade  of  the  previous  50  cm  experiment,  in  which  controlled  beamfront  motion  was 
first  demonstrated.  Injected  electron  beam  characteristics  were  1  MeV,  27  kA,  30  ns.  Tests 
of  the  optical  system,  detailed  in  Table  1,  indicate  that  optical  system  losses  can  be  kept  to 
acceptable  levels  even  in  a  40  mirror  system  if  sufficient  care  in  system  alignment  is  taken. 

In  order  to  measure  the  propagation  of  the  beamfront  down  the  drift  tube,  a  series  of 
current  collection  probes  were  installed  in  the  wall  of  the  drift  tube.  As  the  beamfront 
approaches  these  probes,  the  current  measured  should  rise  gradually  until  the  beamfront 
passes  by  the  probe  location,  at  which  beam  current  at  the  wall  should  fall  rapidly  (see  Fig. 
lc).  Sample  oscilloscope  waveforms  from  the  current  collecting  probes  at  the  drift  tube  wall 
are  shown  in  Fig.  4,  and  show  these  expected  characteristics  clearly. 
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Figure  2:  Experimental  configuration  for  the  beamfront  propagation  study  on  the  laser 
controlled  collective  ion  accelerator. 


Figure  3:  Photograph  of  the  laser-controlled  collective  ion  accelerator  with  the  100  cm  drift 
tube. 


8 


Table  1:  Summary  of  optical  system  tests 


Channel 

Axial 

Position 

(cm) 

Partial 

Mirror 

Reflectivity 

— 

Theoretical 

Laser  Energy 

( Joules ) 

1 

5 

0.95 

0.300 

2 

10 

0.95 

0.285 

3 

15 

0.95 

0.270 

4 

20 

0.95 

0.257 

5 

25 

0.95 

0.244 

6 

30 

0.925 

0.348 

7 

35 

o.}25 

0.322 

8 

40 

0.925 

0.298 

9 

45 

0.925 

0.276 

10 

50 

0.90 

0.339 

11 

55 

0.90 

0.306 

12 

60 

0.89 

0.303 

13 

65 

0.88 

0.294 

14 

70 

0.86 

0.301 

15 

75 

0.83 

0.315 

16 

0.80 

0.307 

17 

85 

0.75 

0.307 

18 

0.67 

0.304 

19 

95 

0.50 

0.309 

20 

— 

0.309 

I 

I 

I 
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Probe  #1:  11.5  cm 
2  volts/ div 
10  ns/div 


Probe  #2:  36.5  cm 
1  voit/div 
10  ns/div 


Probe  #3:  615  cm 

500  mV/div 

10  ns/div 


Probe  #4:  66.5  cm 

500  mV  dn- 
10  ns/div 


Figure  4:  Typical  traces  of  the  charge/current  probe  used  to  measure  the  beamfront  propa¬ 
gation. 
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A  series  of  studies  were  made  of  the  beamfront  propagation  velocity  down  the  drift  tube 
for  various  delays  between  the  firing  of  the  laser  and  the  injection  of  the  electron  beam. 
Studies  were  made  for  two  acceleration  gradients.  Data  for  the  first  gradient,  a  two  part 
gradient  designed  to  take  the  beamfront  velocity  smoothly  from  0.05c  to  about  0.3c  in  a 
distance  of  100  cm,  are  shown  in  Figs.  5-7.  Good  agreement  is  seen  between  the  beamfront 
velocity  and  that  associated  with  the  laser  timing  only,  as  expected,  when  the  laser  is  fired 
just  before  the  beam  is  injected.  If  the  laser  is  fired  much  earlier,  the  drift  tube  becomes 
filled  with  plasma  and  the  beam  propagates  unencumbered  by  space  charge  effects.  Data 
for  a  second  gradient,  a  slower  gradient  with  a  final  velocity  of  only  about  0.2c  is  shown  in 
Figs.  8-10. 

2.1.3  Pseudospark  Discharge  Experiment 

We  have  continued  a  new  experimental  program.  Pseudospark  discharge.  Although,  the 
pseudospark  discharge  has  interesting  switching  characteristics  which  are  being  studied  by 
other  laboratories,  our  current  interest  is  centered  around  the  study  of  electron  beams.  A 
simple  discharge  chamber  of  modular  type  was  constructed  to  investigate  this  new  electron 
beam  source.  The  chamber  was  operated  at  low  voltage  ~  25  kV,  producing  ~10  Hz, 
electron-beam  pulses  of  — - 100  A,  ~10  ns.  The  rms  emittance  of  the  beam  was  measured  and 
found  to  be  ~50  mm-mrad,  which  corresponds  to  the  normalized  brightness  of  the  beam  of 
~  10u  A/(m2rad2).  Papers  detailing  this  study  are  published  in  Proc.  1989  IEEE  Particle 
Accelerator  Conference,  p.  316,  and  AIP  Conference  Proceedings  195.  Dense  Z-Pinches,  p. 
418  (1989),  Appl.  Phys.  Lett.  56,  1746  (1990),  and  to  appear  in  Phvs.  Fluids  B  2,  (10) 
(1990)  and  Proc.  of  Beams  ’90;  copies  are  enclosed  in  Appendix  B. 

2.2  Theoretical  Research  on  Beam  Propagation  Systems 

We  will  discuss  the  progress  made  during  the  past  year  in  four  separate  but  related  topics 
that  were  under  investigation.  The  four  topics  are:  “Initial  Ion  Production  and  Acceleration 
Phase,”  “Laser  Controlled  Beamfront  Accelerator,”  “Intense  Beam  Propagation  Across  a 
Magnetic  Field,”  and  “Ion  Hose  Instability  of  the  Bennett  Profile  Beam.” 
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Figure  5:  Wall  probe  measurements  of  beamfront  propagation  versus  time  for  different  laser 
firing  delay  times  and  without  the  laser  ( “two-part-’  gradient). 
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Figure  9:  \oitage observed  at  the  wall  probes  for  the  different  delay  times  (  "slower  "  gradient). 
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Distance  downstream  in  meters 


Figure  10:  Average  velocity  versus  distance  downstream  for  different  laser  tiring  delay  time 
and  without  the  laser  (“slower"  gradient). 


2.2.1  Initial  Ion  Production  and  Acceleration  Phase 

During  the  past  year,  we  have  completed  numerical  studies  related  to  the  initial  phase  of  the 
ion  production  and  acceleration  process  in  a  localized  gas  cloud  system.  These  results  are 
published  in  the  Journal  of  Applied  Physics,  February  1990  issue  and  a  copy  is  included  in 
Appendix  B. 

In  experiments  in  which  an  intense  relativistic  electron  beam  is  injected  into  an  evacuated 
drift  tube  with  a  localized  gas  cloud  located  near  the  anode,  that  serves  as  the  source  of  ions 
for  effective  beam  propagation,  a  large  fraction  of  the  injected  electron  current  is  observed 
downstream.  These  experiments  have  been  simulated  using  a  particle- in-cell  code  which 
realistically  models  ionization  of  the  gas.  It  was  found  that  when  the  injected  electron  beam 
current  exceeds  the  space-charge  limiting  current,  the  majority  of  the  ions  produced  achieve 
energies  of  the  order  of  the  beam  energy  and  provide  for  an  effective  channel  to  space  charge 
neutralize  the  electron  beam.  There  are  a  few  ions  that  are  accelerated  to  energies  several 
times  the  electron  Hearn  energy  by  the  coherent  motion  of  the  ions  and  the  intense  virtual 
cathode  electric  fields.  The  majority  of  the  ions  allow  for  the  total  beam  current  to  propagate 
to  the  downstream  surface  once  the  ion  channel  has  also  propagated  to  this  location.  The 
dependence  of  the  peak  ion  energy  on  the  system  parameters  as  observed  in  the  simulations 
was  also  examined.  For  the  parameter  regimes  investigated  with  beam  energies  up-  to  3  MV, 
beam  currents  up  to  35  kA,  gas  pressures  up  to  600  mTorr,  and  gas  cloud  widths  up  to  6  cm, 
peak  ion  energies  of  5-6  times  the  electron  beam  energy  have  been  observed  in  the  numerous 
simulations. 

2.2.2  Laser  Controlled  Beamfront  Accelerator-Simulation 

We  have  written  an  electrostatic  particle-in-cell  code  called  LCA  to  simulate  the  laser- 
controlled  acceleration  experiments.  A  copy  of  the  paper  describing  this  code  is  included  in 
Appendix  B.  In  the  simulations,  an  electron  beam  of  voltage  current  /o,  and  radius  Rf, 
is  injected  into  a  grounded  cylindrical  drift  tube  of  length  d  and  radius  Rw  along  the  axis 
of  the  drift  tube.  The  region  of  the  drift  tube  extending  from  the  anode  to  a  distance  r0 
downstream  is  filled  by  hydrogen  gas  at  a  constant  pressure  p0  [see  Fig.  11].  The  electron 
beam  is  assumed  to  be  focussed  by  an  infinitely  strong  guide  magnetic  field,  so  that  particles 
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Figure  11:  Simulation  model. 


in  the  simulation  move  only  along  the  axis  of  the  drift  tube.  The  beam  radius  is  also  assumed 
to  be  much  smaller  than  the  wall  radius  so  that  the  charge  and  current  density  and  the  axial 
electric  field  are  approximately  uniform  across  the  beam  cross-section. 

The  macroparticles  in  the  simulation  obey  the  relativistic  equations  of  motion,  the  electric 
field  Ez  on  axis  is  computed  by  first  computing  the  potential  <f>  at  equally  spaced  grid  points 
on  the  axis,  numerically  computing  the  derivative  at  points  lying  halfway  between  successive 
grid  points,  and  linearly  interpolating  to  obtain  the  value  of  Ez  at  other  points.  Ionization 
of  the  neutral  gas  is  modeled  by  dividing  the  gas  region  into  grid  cells  and  monitoring  the 
amount  of  ionization  in  each  grid  cell  which  is  produced  by  impact  ionization.  When  the  total 
number  of  ions  produced  in  a  grid  cell  exceeds  the  number  of  ions  in  an  ion  macroparticle, 
an  ion  macroparticle  and  an  electron  macroparticle  are  created  at  the  center  of  the  grid.  The 
actual  details  have  been  explained  elsewhere. 

The  time  required  for  the  laser  beam  to  create  plasma  after  striking  the  target  and  the 
time  required  for  the  plasma  once  it  is  created  to  travel  from  the  wall  to  the  center  of  the 
drift  tube  are  assumed  to  be  constant  quantities  which  are  left  unspecified.  The  front  of 
the  laser-produced  plasma  is  assumed  to  sweep  smoothly  from  one  end  of  the  drift  tube  to 
the  other  and  the  only  effect  of  the  plasma  in  the  simulation  is  to  completely  neutralize  any 
space-charge  in  the  drift  tube  behind  the  plasma  front. 

The  first  set  of  results  was  obtained  for  a  900  keV,  20  kA,  1  cm-radius  electron  beam 
which  is  injected  into  a  5  cm-radius,  50  cm-long  drift  tube  with  a  2  cm-wide.  100  mTorr  cloud 
of  hydrogen  gas  located  next  to  the  injection  plane.  The  front  of  the  laser-produced  plasma 
is  assumed  to  travel  down  the  drift  tube  at  a  velocity  which  increases  linearly  from  0.04c  to 
0.2c  over  a  distance  of  45  cm.  These  parameters  are  the  same  as  those  in  the  experiment 
described  above.  Figure  12  shows  the  peak  proton  energy  measured  45  cm  downstream  from 
the  injection  plane  as  a  function  of  the  time  delay  between  the  start  of  the  beam  pulse  and 
the  start  of  the  laser  pulse.  For  these  runs  it  was  assumed  that  the  laser-produced  plasma 
was  created  10  ns  after  the  laser  beam  struck  the  target  and  that  the  time  required  for  the 
plasma  to  travel  from  the  wall  to  the  center  of  the  drift  tube  was  equal  to  the  time  required 
for  a  proton  to  travel  from  the  wall  to  the  surface  of  the  electron  beam,  assuming  that  the 
potential  depression  produced  by  the  beam  was  equal  to  the  beam  energy  V0.  The  figure 
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shows  that  the  design  energy  of  18.8  MeV  is  attained  over  a  broad  range  of  time  delays.  If 
one  compares  the  phase-space  trajectory  of  the  peak-energy  proton  macroparticle  and  the 
phase-space  trajectory  of  the  laser- produced  plasma  front,  we  find  for  this  design  that  the 
accelerated  proton  tracks  the  laser  beam  trajectory  closely. 

When  the  length  of  the  drift  tube  is  doubled  from  50  cm  to  100  cm  (with  all  other 
system  parameters  the  same)  and  the  same  velocity  gradient  is  used,  the  peak  proton  energy 
measured  at  the  end  of  the  drift  tube  falls  short  of  the  design  value,  i.e..  the  original  velocity 
gradient  cannot  be  extended  to  longer  distances.  The  reason  the  original  velocity  gradient 
cannot  be  scaled  to  longer  distances  can  be  seen  as  follows.  The  equation  of  motion  of  a 
proton  which  is  being  accelerated  by  an  electric  field  E:  is  dv/dt  =  (e/m)Ez.  which  can  be 
rewritten  as  Ez  —  (me2 j e)vdv / dz.  If  the  velocity  gradient  dv/dz  used  in  the  above  runs  is 
substituted  into  the  preceding  expression,  we  find  tha*  electric  field  needed  to  accelerate 
the  proton  at  the  desired  velocity  gradient  is  approximately  Ez  =  334 3  MV/m,  where 
3  =  v/c,  e.g.,  for  3  =  0.3,  an  electiic  field  of  more  than  100  NlV/m  is  needed.  In  Fig.  13,  we 
have  plotted  the  magnitude  versus  the  location  o>.  peak  electric  field  tor  the  run  with  a 
100  cm-long  drift  tube  in  which  the  greatest  peak  proton  energy  was  measured.  Notice  that 
as  the  beam  front  moves  downstream,  the  peak  electric  field  tends  to  fall  until  it  is  no  longer 
large  enough  to  continue  accelerating  the  proton  at  a  constant  velocity  gradient. 

The  tendency  of  the  peak  electric  field  to  fall  as  the  beam  front  moves  downstream 
suggests  that  it  may  be  better  to  accelerate  the  protons  with  a  steep  velocity  gradient  at  the 
start  and  then  taper  the  gradient  as  the  beam  front  moves  downstream.  In  the  second  series 
of  runs,  we  therefore  chose  to  do  a  set  of  runs  in  which  the  beam  front  velocity  increases 
linearly  with  distance  (i.e.,  dv/dz  =  constant)  until  it  reaches  a  transition  point,  after  which 
the  beam  front  velocity  increases  linearly  with  time  (i.e.,  dv/dt  =  constant). 

The  plasma  front  velocity  at  the  anode  plane  was  chosen  to  be  0.04c,  the  plasma  front 
velocity  at  the  downstream  end  of  the  drift  tube  was  chosen  to  be  0.4c.  and  the  velocity  3tc 
at  the  transition  point  zt  was  varied.  Figure  14  shows  the  peak  proton  energy  measured  at 
90  cm  as  a  function  of  3t  for  three  different  values  of  zt.  Note  that  Vo  =  1.5  MV  in  this  run. 
The  figure  shows  that  by  adjusting  3t ,  energies  in  excess  of  the  design  energy  of  60.7  MeV 
can  be  achieved. 
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Figure  13:  Magnitude  versus  location  of  peak  electric  Held. 
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Figure  14:  Peak  proton  energy  measured  at  90  cm  versus  transition  velocity  A 


It  should  be  noted  that,  as  might  be  expected,  the  peak  proton  energy  depends  strongly 
on  the  beam  energy,  since  the  peak  electric  field  increases  with  increasing  beam  energy  In 
Fig.  15,  we  plot  the  peak  proton  energy  measured  at  the  downstream  end  of  a  50  cm-long 
drift  tube  as  a  function  of  the  electron  beam  energy.  In  each  case  the  same  velocity  profile 
is  used.  For  beam  energies  greater  than  1.4  MeV,  peak  proton  energies  exceed  the  design 
value.  As  the  beam  energy  decreases,  the  peak  proton  energy  decreases  sharply,  e.g.,  at  1 
MeV,  the  peak  proton  energy  is  already  less  than  20  MeV. 

Optimal  Acceleration-New  Design.  The  largest  proton  energies  achievable  in  any 
given  system  can  be  calculated  by  integrating  the  peak  electric  fields  obtained  from  the 
simulations  as  a  function  of  distance  (assuming  that  the  magnitude  of  the  peak  electric  field 
depends  only  on  the  position  of  the  plasma  beam  front  and  not  on  its  velocity).  The  required 
variation  of  the  plasma  front  position  with  time  can  be  obtained  by  integrating  the  equations 
rdv/dz  =  eE2%max(z)  and  dz/dt  =  v[z ),  where  Ez<max(z)  is  the  peak  electric  field  at  r.  Figure 
16  shows  3  vs.  r  and  t  vs.  r  for  a  1.2  MeV,  20  kA  electron  beam  which  is  injected  into  a  50 
cm-long  drift  tube. 

We  have  developed  the  electrostatic  PIC  code  LCA  to  simulate  beam  propagation  and 
collective  ion  acceleration  in  a  laser- controlled  accelerator.  Although  the  code  only  crudely 
models  the  experiment,  it  has  been  able  to  reproduce  some  of  the  experimental  results  and 
may  be  useful  as  a  design  tool  for  future  laser-controlled  acceleration  experiments. 


2.2.3  Intense  Beam  Propagation  Across  a  Magnetic  Field 

We  completed  the  studies  on  an  initial  model  of  the  propagation  of  an  intense  electron-ion 
beam  across  a  transverse  magnetic  field.  The  results  are  published  in  the  1989  Particle 
Accelerator  Conference  Proceedings,  and  a  copy  of  the  work  is  included  in  Appendix  B. 

Previous  theoretical  studies  have  shown  the  existence  of  a  self-consistent  downstream 
Bennett  equilibrium  for  the  electrons  and  ions  when  no  applied  magnetic  field  is  present 
[C.D.  Striffler.  R.L.  Yao,  X.  Zhang,  Proc.  of  the  1987  IEEE  PAC.  page  975].  We  have  related 
these  downstream  properties  to  the  diode  voltage,  the  transmitted  electron  beam  current, 
and  the  ion  properties  in  the  localized  gas  cloud  region.  For  our  experimental  parameters. 
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Figure  15.  Peak  proton  energy  measured  at  50  cm  versus  electron  beam  energy 
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Figure  16:  Optimal  acceleration  profile. 

\ o  =  1  MV.  7  =  5  kA,  the  downstream  self-pinched  equilibrium  state  is  composed  of  cold 
ions  with  low  axial  speed  (<  0.05  c).  and  electrons  with  a  temperature  of  about  60  keV  and 
an  axial  speed  of  about  0.S4  c.  The  nearly  charge- neutral  beam  system  has  been  shown  to 
effectively  propagate  up  to  the  diode  current  of  20  kA  where  the  electron  temperature  is 
predicted  to  be  about  200  keV.  The  5  kA  case  was  chosen  for  examining  the  eifects  of  a 
transverse  magnetic  field  on  the  electron-ion  beam  system,  mainly  because  of  experimental 
reproducibility.  We  considered  the  propagation  of  this  intense  electron-ion  beam  across  an 
applied  magnetic  field.  We  find  that  in  the  intense  beam  regime,  the  propagation  is  limited 
due  to  space-charge  depression  caused  by  the  deflection  of  the  electron  beam  by  the  transverse 
field.  This  critical  field  is  of  the  order  of  the  peak  self-magnetic  field  of  the  electron  beam 
which  is  substantially  higher  than  the  single  particle  cutoff  field. 

2.2.4  Ion  Hose  Instability  of  the  Bennett  Profile  Beam 

The  self-consistent  Vlasov  equilibrium  that  leads  to  the  Bennett  radial  density  profile  for  the 
downstream  electron-ion  beam  system  lias  been  analyzed  relative  to  the  ion  hose  instability. 
Specifically,  we  consider  the  limiting  case  where  the  ion  beam  species  is  cold,  stationary,  and 
provides  charge  neutralization.  The  most  simple  model  assumes  that  each  beam  is  rigid.  We 
then  investigate  the  effects  due  to  the  anharmonic  nature  of  the  potential  well  using  a  '  pread 
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density-1  model  and  the  effects  due  to  an  axial  thermal  velocity  spread  of  the  electron  beam. 
In  addition,  nonlinear  effects  due  to  large  beam  displacements  are  numerically  investigated. 
For  the  rigid  beam  model,  the  linearized  equations  for  transverse  motion  of  the  beams 


are 
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and  ze  =  u0,  i,  =  0.  We  have  defined  u>e  =  sje^rito/ -  \Je2ne 0/eom,, 
nt0  =  neo  =  n0,  and  the  displacements  x  are  assumed  to  be  much  less  than  the  beam  cross- 
sectional  dimensions  a.  For  harmonic  displacements,  we  find  that  resonant  growth  occurs 
for  wavelengths 
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where  the  growth  rate  is  /mu.  The  spectrum  of  growth  rates  is  shown  in  Fig.  17  by  the 
solid  line  (R.B.),  where  the  strong  resonant  growth  is  clearly  seen.  The  system  parameters 
for  this  case  are  V0  =  1  MV,  I0  =  20  kA,  30  =  v0/c  =  0.65.  a  =  1  cm,  and  Rw/a  =  10. 

The  effects  due  to  the  Bennett  profile  which  results  in  particles  oscillating  at  different 
transverse  frequencies  are  modelled  by  annular  layers  with  different  densities.  After  weighting 
by  density  and  averaging  to  obtain  the  net  beam  displacement,  we  obtain  a  dispersion  relation 
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1  hese  results  are  displayed  in  Fig.  17  by  the  dashed  curve  labelled  S.D.  We  see  that  the  spread 
in  transverse  oscillation  frequency  (weighted  by  density)  results  in  an  unstable  spectrum  up 
to  the  peak  on  axis-value  of  versus  up  to  the  average  value  of  uJ,/\/E  for  the  rigid 

beam  (R.B.)  model. 

The  effects  due  to  an  axial  thermal  velocity  spread  of  the  electron  beam  species  is  mod¬ 
elled  by  the  one-dimensional  relativistic  Maxwellian  distribution.  The  beam  is  assumed  to 
be  composed  of  many  rigid  disks  having  the  same  density  profile  but  different  axial  velocities. 
These  results  are  shown  by  the  dashed  curve  in  Fig.  17  labeled  T.Y. 
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Figure  17:  Ion  hose  instability  growth  rates  versus  wavenumber.  System  parameters  are: 

\n  =  1  MV.  I0  =  20  kA.  R.B.  -  Rigid  beam  model:  S.D.-  Spread  density  model  (anharmonic 
potential  well):  T.V.  -  Thermal  velocity  spread  model. 


In  summary,  our  studies  of  the  ion  hose  instability  of  a  Bennett  profile  beam  indicate: 

•  From  the  rigid  beam  model,  the  instability  exists  in  the  long  wavelength  region  and  is 

weakly  absolute.  The  resonant  growth  occurs  at  A  =  . 

•  Due  to  the  anharmonic  potential  as  well  as  axial  velocity  spread,  the  maximum  growth 
is  reduced,  the  unstable  region  is  broadened,  and  the  instability  becomes  convective. 

•  For  the  case  of  large  beam  displacements,  the  nonlinear  equations  of  motion  for  the 
ion  hose  oscillations  have  been  derived  for  the  Bennett  profile.  As  the  beam  separation 
becomes  of  the  order  of  the  beam  radius,  the  instability  growth  decreases  and  saturates. 

2.3  Studies  of  Pulsed  Power  Systems  and  Diagnostics 

2.3.1  Circuit  Analysis  of  Pulsed  Power  Systems 

We  have  carried  out  a  comprehensive  circuit  analysis  of  basic  pulsed  power  systems.  The  elec¬ 
trical  pulsed  power  ch^uas  aie  classified  into  two  main  categories,  capacitive-  and  inductive- 
energy  storage  systems.  Four  basic  circuits  are  found  in  each  category,  and  they  are  further 
classified  into  lumped-  and  distributed-circuits.  It  is  noted  that  the  inductive  energy  sys¬ 
tems  are  the  duals  of  the  capacitive  energy  systems,  and  vice  versa.  We  proposed  two  new 
inductive  energy  pulsed  power  systems,  which  are  <  he  duals  of  LC  generator  and  Blumleiu 
line.  Two  papers  detailing  this  study  are  published  in  Proc.  1989  IEEE  Particle  Accelerator 
Conference,  p.  1954  (1989),  and  Journal  of  Applied  Physics  £7,  4333  (1990);  copies  are 
enclosed  in  Appendix  B. 

2.3.2  Beam  Emittance  in  a  Linear  System  and  a  Nonlinear  Lens 

In  this  study,  we  have  shown  that  for  a  beam  of  arbitrary  distribution  function,  the  rms 
emittance  is  invariant  under  any  linear  transformation  if  the  determinant  of  the  transfer 
matrix  is  unity,  and  the  normalized  rms  emittance  is  invariant  with  a  nonunity  determinant, 
which  is  associated  with  an  axial  acceleration  of  the  beam.  An  exact  rms-emittance  increase 
when  the  beam  passes  through  a  nonlinear  lens  with  a  spherical  aberration  is  derived.  Two 
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papers  detailing  this  study  are  published  in  the  Physics  of  Fluids  29,  3495  (1986),  and  Physics 
of  Fluids  B  2,  452  (1990);  copies  are  enclosed  in  Appendix  B. 
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Detailed  experimental  studies  of  a  collective  acceleration  experiment  in  which  a  time- 
sequenced  laser-generated  ionization  channel  is  used  to  control  the  propagation  of  an  intense 
relativistic  electron  beamfront  are  presented.  Ions  trapped  in  the  potential  well  at  the 
beamfront  are  accelerated  as  the  velocity  of  the  beamfront  is  increased  in  a  manner  controlled 
by  the  time-dependent  axial  extent  of  the  ionization  channel.  Beamfront  propagation  data  for 
two  different  accelerating  gradients  are  presented,  together  with  results  of  ion  acceleration 
studies  for  both  gradients. 


I.  INTRODUCTION 


II.  THEORETICAL  DISCUSSION 


Interest  in  collective  acceleration  has  arisen  because  the 
relatively  high  electric  fields  ( 100-1000  MV/m)  that  can  be 
generated  in  relativistic  electron  beam  systems  are  higher 
than  can  be  easily  maintained  in  conventional  accelerator 
structures.1 2  A  typical  collective  accelerator  involves  the  in¬ 
jection  of  an  intense  relativistic  electron  beam  into  a  source 
of  positive  ions  to  be  accelerated.  Ion  sources  used  have  in¬ 
cluded  neutral  gas-filled  drift  tubes,3  4  a  neutral  gas  cloud 
localized  to  the  injection  region,5  7  or  a  dielectric  anode  ma¬ 
terial  (the  Luce  diode8  <) ) .  In  each  case,  if  the  beam  is  injected 
at  a  sufficiency  high  current  level,  a  virtual  cathode  forms 
immediately  downstream  of  the  injection  point  and  ions  can 
be  produced  and  trapped  by  the  potential  well  associated 
with  the  virtual  cathode.  Ion  acceleration  to  energies  in  ex¬ 
cess  of  the  electron  beam  energy  can  occur  as  the  virtual 
cathode  moves  downstream  of  the  injection  point,  usually  as 
a  result  of  a  gradual  neutralization  of  the  electron  beam 
space  charge  by  increasing  numbers  of  ions  generated  by 
electron-impact  and  ion-ion  avalanche  ionization  processes. 

Although  considerable  progress  has  been  achieved  in 
the  study  of  such  “naturally  occurring”  collective  accelera¬ 
tion  processes  in  recent  years,3"10  the  goal  of  controlling  the 
acceleration  process  over  longer  distances  has  proven  an  elu¬ 
sive  one.  To  date,  the  most  successful  experiments  in  this 
regard  have  been  the  IFA-1  and  IFA-2  ionization  front  ac¬ 
celerator  experiments  of  Olson."  12  Recently,  our  group  re¬ 
ported  initial  successful  operation  of  a  laser-controlled  col¬ 
lective  ion  acceleration  experiment,  based  on  the  same 
fundamental  concepts  as  the  IFA  experiments  but  signifi¬ 
cantly  different  in  both  design  and  operation. 13  In  this  paper, 
detailed  experiments  conducted  in  an  attempt  to  understand 
the  potentials  and  limitations  of  this  new  collective  accelera¬ 
tor  configuration  are  reported,  including  detailed  measure¬ 
ments  of  controlled  beamfront  motion  in  experimental  con¬ 
figurations  designed  for  two  different  accelerating  gradients. 
Results  of  ion  acceleration  expenments  utilizing  both  gradi¬ 
ents  are  reported  and  compared  with  theoretical  expecta¬ 
tions.  In  Sec.  II  of  this  paper  the  accelerator  concept  is  de¬ 
scribed  and  basic  design  considerations  are  discussed. 
Expenments  conducted  to  investigate  the  concept  are  de- 
scnbed  in  Sec.  Ill,  and  conclusions  are  drawn  in  Sec.  IV. 


The  basic  concept  behind  the  experiment  is  shown  in 
Fig.  1.  An  intense  relativistic  electron  beam  (900  keV,  20 
kA,  30  ns)  is  injected  into  an  evacuated  drift  tube  at  a  cur¬ 
rent  level  several  times  the  vacuum  space-charge  limit,  given 
approximately14  by 

17  000  (y273-  1)3/2 

I,  = - — - - -  (A),  (1) 

[1  +  21n(6/fl)](l -j) 

where  b  is  the  drift  tube  radius,  a  is  the  beam  radius,  y0  is  the 
relativistic  mass  factor  for  the  electrons  at  injection,  and 
f =  n,/ne  represents  any  charge  neutralization  provided  by 
positive  ions.  As  indicated  in  Fig.  1(a),  a  virtual  cathode 
forms  at  the  injection  point  with  a  depth  approximately 
equal  to  the  anode-cathode  potential  difference.  The  axial 
position  of  the  virtual  cathode  downstream  of  the  anode 
plane  is  usually  on  the  order  of  the  anode  cathode  gap,15  so 
that  megavolt  potentials  at  the  virtual  cathode  are  formed 
only  millimeters  away  from  the  grounded  anode.  It  is  this 
very  high  electric  field  on  the  upstream  side  of  the  virtual 
cathode  that  collective  accelerators  usually  seek  to  exploit 
for  ion  trapping  and  acceleration. 
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FIG.  1 .  Conceptual  design  of  laser-controlled  accelerator. 
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If  the  beam  is  injected  into  a  localized  gas  cloud  [  Fig. 
1(b)],  ionization  processes  can  quickly  build  up  sufficient 
ion  density  to  neutralize  the  electron  beam  space  charge  and 
the  virtual  cathode  can  move  downstream  to  the  edge  of  the 
gas  cloud.  This  motion  can  result  in  the  acceleration  of  a  few 
ions  to  energies  considerably  higher  than  that  associated 
with  the  depth  of  the  potential  well  at  the  virtual  cathode. 

In  order  to  control  the  motion  of  the  virtual  cathode 
over  distances  greater  than  a  few  centimeters,  however,  a 
means  of  providing  an  ionization  channel  whose  axial  extent 
can  be  controlled  as  a  function  of  time  is  required  [Fig. 
1(c)].  In  the  present  experiments,  shown  schematically  in 
Fig.  2,  the  ionization  channel  is  generated  by  time-sequenced 
laser-target  interactions.  A  laser  pulse  is  separated  into  ten 
approximately  equal  energy  beams  which  are  then  optically 
delayed  over  different  path  lengths.  The  laser  light  then  va¬ 
porizes  and  ionizes  a  target  material  on  the  drift  tube  wall, 
and  ions  drawn  into  the  beam  by  the  electron  space  charge  at 
the  beamfront  provide  the  required  time-sequenced  channel 
of  ionization  to  control  beamfront  mot  ion. 

This  basic  concept  of  operation  results  in  a  number  of 
constraints  on  the  design  of  the  experiment  The  require¬ 
ment  that  the  injected  beam  current  be  several  times  the 
vacuum  space-charge  limit  imposes  a  lower  bound  on  the 
drift  tube  radius  of  about  twice  the  beam  radius.  An  upper 
bound  for  the  drift  tube  radius  is  set  by  the  requirement  that 
ions  created  at  the  drift  tube  wall  be  drawn  into  the  beam  in  a 
time  less  than  or  equal  to  the  delay  time  between  two  adja¬ 
cent  laser  pulses. 

If  the  virtue!  cltl;  ode  is  accelerated  continuously  from  a 
low  initial  velocity  to  a  higher  final  one,  ions  will  be  trapped 
within  the  potential  well  at  the  beamfront  only  if  the  relative 
velocity  of  the  well  and  the  ions  is  less  than  a  value  given  by 

".<(2  eVo/m,)in  (2) 

where  Vn  is  the  depth  of  the  potential  well  associated  with 
the  virtual  cathode.  This  result  yields  vu,  —  v,  <  0.44  c  for  a 
900-keV  electron  beam.  Another  constraint  upon  the  accel¬ 
eration  process  is  that  the  electric  field  associated  with  the 
virtual  cathode  at  the  beamfront  must  be  large  enough  that 
ions  remain  trapped  during  the  acceleration  process.  This 
imposes  a  limit  on  the  accelerating  gradient  given  by 

v ,  rse.  <.eE  /  m, .  (3) 
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TIG  2.  Experimental  configuration. 


Although  electric  fields  associated  with  the  virtual  cathode 
at  the  injection  point  can  be  as  high  as  100-1000  MV/m, 
beam  dynamics  at  the  beamfront  as  the  beam  propagates 
downstream  may  reduce  these  field  strengths  considerably. 
In  addition,  oscillations  of  the  trapped  ions  in  the  potential 
well  at  the  beamfront  and  fluctuations  in  the  depth  of  the 
potential  well  should  also  be  considered. 

In  addition  to  these  considerations,  the  laser-target  ge¬ 
ometry  must  be  carefully  designed  such  that  electron  beam 
ionization  of  the  target  is  negligible.  Only  if  this  requirement 
is  satisfied  will  effective  control  over  the  beamfront  motion 
by  the  laser  be  possible. 

III.  EXPERIMENT 
A.  Apparatus 

The  basic  experimental  configuration  is  shown  in  Fig.  2. 
An  intense  relativistic  electron  beam  ( 900  keV,  20  kA,  30  ns 
FWHM)  is  emitted  from  a  4-mm-diam  tungsten  cathode 
located  1  cm  upstream  of  a  stainless-steel  anode.  A  1 4-mm- 
diam  hole  in  the  anode  plate  allows  almost  all  of  the  beam 
current  in  the  diode  to  pass  into  the  downstream  drift 
chamber.  A  drift  chamber  diameter  of  10  cm  was  chosen  on 
the  basis  of  both  theoretical  calculations  described  in  the 
previous  section  and  experimental  measurements  of  beam 
propagation  in  dielectric-lined  drift  tubes  of  various  diame¬ 
ters.  16  Seed  protons  for  acceleration  were  produced  by  elec¬ 
tron-impact  and  ion-ion  avalanche  ionization  of  a  localized 
hydrogen  gas  cloud  produced  by  a  fast  gas  puff  valve.  At  the 
time  of  electron  beam  injection,  the  localized  gas  cloud  wi, 
measured  to  extend  no  more  than  2  cm  into  the  vacuum  drift 
tube.  The  laser  target  geometry  chosen  consisted  of  a  5-mm- 
wide  strip  of  0.8-mm-thick  polyethylene  (CH2)  sheet 
mounted  in  a  stainless-steel  bracket  that  allowed  for  axial 
movement  of  the  target  if  necessary  to  provide  a  fresh  target 
spot  for  the  laser  pulse.  Experimental  tests  have  confirmed 
that  electron  beam  ionization  of  the  laser  target  material  is 
not  sufficient  to  allow  beam  propagation  downstream.16  A 
^-switched  ruby  laser  (6  J,  15  ns)  was  divided  into  ten  ap¬ 
proximately  equal  energy  beams  and  optically  delayed  to 
provide  a  time-sequenced  source  of  ions  down  the  50  cm 
length  of  the  drift  tube.  Optical  tests  of  the  system16  have 
demonstrated  that  actual  laser  energy  on  the  target  varied 
from  a  maximum  of  0.55  J  at  the  first  target  spot  to  a  mini¬ 
mum  of  0.32  J  at  the  last,  with  differences  primarily  due  to 
optical  losses  in  the  system.  Each  laser  beam  was  focused  to  a 
diameter  ofless  than  1  mm  on  the  target  surface  by  a  focus¬ 
ing  lens,  ensuring  adequate  laser  power  density  (about  1-2 
GW/cm:)  for  plasma  production. n 

The  following  diagnostics  were  employed  to  monitor 
the  control  of  the  electron  beamfront  and  the  subsequent 
acceleration  of  ions  in  the  iccelerator.  A  Faraday  cup  cur¬ 
rent  collector  was  installed  at  the  downstream  end  of  the 
drift  tube  to  measure  the  total  beam  current  propagated  to 
this  point.  Five  current  collecting  wall  probes  were  also  in¬ 
stalled  at  axial  positions  10.  20.  30.  40.  and  50  cm  down¬ 
stream  of  the  anode  plane  to  measure  beam  current  deflected 
to  the  drift  tube  wall  at  the  beamfront  as  a  function  of  time. 
The  current  probes  consisted  of  open-ended  50-H  semirigid 
coaxial  cables  with  an  effective  collection  area  of  about  2 
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mnr.The  energy  of  accelerated  protons  was  determined  us¬ 
ing  stacked  foil  nuclear  activation  techniques.  Titanium 
[Ti47(p,/i)V47]  and  copper  [Cu',3(p,n)Zn63]  reactions  were 
used  having  threshold  energies  of  3.7  and  4.2  MeV,  respec¬ 
tively. 

Neutrons  produced  by  the  interaction  of  accelerated 
protons  with  the  target  foils  and  with  the  drift  tube  walls 
were  also  measured  using  a  silver  activation  neutron  detec¬ 
tor. 

B.  Results 

Experiments  were  conducted  using  two  sets  of  optical 
delays  designed  to  smoothly  accelerate  the  electron  beam- 
front  from  0.05  to  0.2  c  and  from  0.05  to  0.3  c  over  a  distance 
of  45  cm.  These  corresponded  to  accelerating  gradients  of 
approximately  40  and  90  MV/m.  respectively.  Results  of 
measurements  of  the  beamfront  motion  down  the  drift  tube 
as  a  function  of  time  as  determined  from  the  arrival  time  of 
electron  current  pulses  at  the  wall  probes  are  shown  in  Figs. 
3  and  4.  In  Fig.  3  is  also  shown  typical  beamfront  propaga¬ 
tion  data  for  the  case  where  the  laser  is  fired  200  ns  in  ad¬ 
vance  of  the  beam,  and  the  beam  propagates  freely  down¬ 
stream  neutralized  by  a  preformed  plasma  channel,  and  for 
the  case  where  the  laser  is  not  fired  at  all,  in  which  case  beam 
propagation  downstream  can  only  occur  as  rapidly  as  accel¬ 
erated  ions  from  the  gas  cloud  propagate  downstream.  Re¬ 
sults  for  the  case  where  the  laser  is  swept  along  the  target  just 
as  the  electron  beam  is  injected  into  the  drift  tube  indicate 
that  effectiv#  control  is  exercised  by  the  laser  in  good  agree¬ 
ment  with  design  values  for  both  accelerating  gradients 
(Fig.  4).  Peak  currents  measured  by  the  wall  probes  both 
with  and  without  the  laser  are  plotted  in  Fig.  5,  and  indicate 
that,  as  expected,  the  laser  allows  more  effective  current 
propagation  downstream  than  is  possible  in  vacuum,  even 
when  wall  currents  are  used  as  an  indicator  of  beam  trans¬ 
port. 

Data  from  Faraday  cup  measurements  of  the  total  elec- 
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FIG  3  Wall  probe  measurements  of  bcamiront  propagation  vs  time  for 
optimum  laser  finng.  laser-tired  ?00  ns  in  advance  of  electron  beam  injec¬ 
tion.  and  without  laser  ( 40  MV/m  accelerating  gradient  l 
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FIG  4.  Wall  probe  measurements  of  beamfront  propagation  vs  time  for 
optimum  laser  firing  delay.  40  MV/m  accelerating  gradient  and  dQ  MV/m 
accelerating  gradient 

tron  beam  current  reaching  the  downstream  end  of  the  drift 
tube  as  well  as  neutron  counts  and  measured  ion  energies  for 
the  two  different  accelerating  gradients  are  shown  in  Figs.  6 
and  8  as  a  function  of  the  delay  time  between  the  finng  of  the 
laser  and  the  firing  of  the  electron  beam.  It  is  easily  seen  from 
the  data  that  firing  the  laser  too  early  in  time  simply  allows 
all  of  the  beam  current  to  flow  rapidly  downstream  as  indi¬ 
cated  from  the  wall  current  probes.  Finng  the  laser  too  late, 
which  is  equivalent  to  not  finng  the  laser  at  all,  results  in 
greatly  reduced  currents  collected  by  the  Faraday  cup,  as 
expected.  It  is  clear  that  the  firing  delays  corresponding  to 
the  transition  from  total  beam  current  propagation  (ions 
produced  too  early  with  regard  to  beam  injection )  to  little 
beam  propagation  ( ions  produced  too  late  to  provide  neu¬ 
tralization)  are  those  that  should  result  in  control  of  the 
virtual  cathode  motion  down  the  dnft  tube.  It  is  therefore 
expected  that  high-energy  ions  should  be  only  observed 
within  this  range  of  firing  delays.  Data  for  the  40  MV/m 
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FIG  5  Peak  wall  current  pulse  amplitudes  vs  axial  position  for  both  opti¬ 
mum  laser  finng  delay  and  no-laser  expenments  Collection  area  of  wall 
probes  is  about  2  mrrr 
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experiments  (Fig.  6)  clearly  show  an  increase  in  observed 
peak  proton  energies  at  this  point  to  values  consistent  with 
design  expectations.  About  10s  protons/pulse  were  acceler¬ 
ated  to  the  peak  energy  of  about  18  MeV.  Although  the  foil 
activation  results  do  not  easily  yield  an  ion  energy  spectrum, 
the  data  are  consistent  with  a  spectrum  having  a  substantial 
peak  at  about  18  MeV,  as  shown  in  Fig.  7.  Data  for  the  90 
MV/m  experiments  (Fig.  8),  however,  show  some  enhance¬ 
ment  of  accelerated  proton  energies  over  the  no-laser  results, 
but  not  to  the  designed  value  of  about  41  MeV. 


IV.  CONCLUSIONS 

Effective  control  of  a  relativistic  electron  beamfront  by  a 
laser-produced  time-sequenced  ionization  channel  has  clear- 
!y  bc~n  demonstrated  for  two  different  accelerating  gradi¬ 
ents.  In  addition,  beamfront  motion  without  the  laser-pro¬ 
duced  ionization  channel  is  consistent  with  theoretical 
expectations,  as  is  the  rapid  propagation  observed  when  the 
ionization  channel  is  produced  well  in  advance  of  electron 
beam  injection. 

Controlled  collective  acceleration  of  ions  at  a  rate  of  40 
MV/m  over  a  distance  of  50  cm  has  also  been  demonstrated, 
but  acceleration  at  the  higher  gradient  (90  MV/m)  has  not 
been  demonstrated  to  date.  This  is  undoubtedly  attributable 
to  the  difficulty  in  maintaining  the  required  >90  MV/m 
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FIG  6  Propagated  current,  neutron  production,  and  peak  proton  energy  as 
J  function  of  firing  delay  between  the  laser  pulse  and  first  rise  of  electron 
beam  pulse  in  the  diode  in  nanoseconds.  40  MV/m  accelerating  gradient 
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FIG.  7.  Counts  from  activated  2-mil  titanium  foils  vs  foil  number  in  the 
target  stack,  40  MV/m  accelerating  gradient. 


electric  fields  at  the  virtual  cathode  over  the  entire  50-cm 
accelerating  distance.  Higher  accelerated  ion  energies  will 
require,  therefore,  either  longer  accelerating  distances  and 
therefore  longer  electron  beam  pulse  durations,  or  alterna¬ 
tively  higher  injected  beam  currents  and  energies  to  main¬ 
tain  higher  electric  fields  at  the  virtual  cathode.  Future  ex¬ 
periments  are  planned  to  investigate  both  of  these 
possibilities. 
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FIG.  8.  Propagated  current,  neutron  production,  and  peak  proton  energy  as 
a  function  of  firing  delay  between  the  laser  pulse  and  first  rise  of  electron 
beam  pulse  in  nanoseconds.  90  MV/m  accelerating  gradient 
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In  experiments  in  which  an  intense  relativistic  electron  beam  is  injected  l.uo  an  evacuated  drift 
tube  with  a  localized  gas  cloud  located  near  the  anode,  ions  with  energies  several  times  the 
electron  beam  energy  have  been  observed.  These  experiments  have  been  simulated  using  a 
panicle-in-cell  code  which  realistically  models  ionization  of  the  gas.  It  was  found  that  when 
the  injected  electron  beam  current  exceeds  the  space-charge  limiting  current,  ions  are 
accelerated  to  energies  several  times  the  electron  beam  energy  by  coherent  motion  of  the  ions 
and  the  intense  virtual  cathode  electric  fields.  The  dependence  of  the  peak  ion  energy  on  the 
system  parameters  as  observed  in  the  simulations  is  also  discussed.  For  the  parameter  regimes 
investigated  with  beam  energies  up  to  3  MV,  beam  currents  up  to  35  kA,  gas  pressures  up  to 
600  mTorr,  and  gas  cloud  widths  up  to  6  cm,  peak  ion  energies  of  5-6  times  the  electron  beam 
energy  have  been  observed. 


I.  INTRODUCTION 

Because  the  electric  fields  produced  by  an  intense  rela¬ 
tivistic  electron  beam  (IREB)  are  much  larger  than  those 
that  can  be  attained  in  conventional  accelerators,  consider¬ 
able  work  has  been  devoted  to  using  IREBs  to  accelerate  ions 
to  high  energies.  In  1 970,  Graybill  and  Uglum 1  reported  that 
ions  with  peak  energies  of  f ?  ur  times  the  electron  beam  ener¬ 
gy  were  produced  when  an  IREB  was  injected  into  a  gas- 
filled  drift  tube.  In  1973,  Luce.  Sahlin.  and  Crites2  reported 
that  ions  with  energies  several  times  the  electron  beam  ener¬ 
gy  were  produced  when  an  IREB  was  injected  into  an  evacu¬ 
ated  drift  tube  using  a  diode  with  a  dielectric  insert  in  the 
anode.  Since  these  results  were  published,  collective  ion  ac¬ 
celeration  has  been  reported  for  many  other  systems,  includ¬ 
ing  those  in  which  the  ions  are  produced  by  beam  ionization 
of  dielectric  drift  tube  walls1  or  by  laser  ionization  of  a  solid 
target.4 

Recently,  Destler,  Floyd,  and  Reiser  reported  that  in 
experiments  performed  at  the  University  of  Maryland  col¬ 
lective  ion  acceleration  was  observed  when  an  IREB  was 
injected  into  an  evacuated  drift  tube  after  gas  had  been 
puffed  into  the  drift  tube  through  a  nozzle  located  near  the 
anode.  A  diagram  of  the  University  of  Maryland  experiment 
is  shown  in  Fig.  1.  The  system  consists  of  a  foilless  electron 
beam  diode,  a  grounded  cylindrical  drift  tube,  a  nozzle  locat¬ 
ed  near  the  anode  through  which  gas  is  puffed  into  the  drift 
tube,  and  magnetic  field  coils,  which  produce  a  guide  field 
for  the  electron  beam.  In  a  typical  experiment,  the  drift  tube 
is  first  pumped  down  to  vacuum  (5  \  10 ' '  Torr).  The  puff 
value  connected  to  the  nozzle  is  then  fired,  injecting  neutral 
gas  into  the  drift  tube  through  the  nozzle  and  producing  a 
cloud  of  gas  extending  no  more  than  a  few  cm  from  the  an¬ 
ode.  The  diode  is  then  fired,  injecting  the  beam  into  the  drift 
tube  through  the  cioud  of  gas.  Current  collectors  placed  at 
the  downstream  end  of  the  drift  tube  measure  the  ion  current 
and  stacked-foil  neutron  activation  analysis  is  used  to  deter¬ 
mine  the  ion  energy  and  flux. 

In  the  experiments,  it  was  found  that  when  the  injected 


current  Jn  exceeded  the  space-charge  limiting  current  1 L  for 
the  system,  ions  were  accelerated  to  energies  several  times 
the  electron  beam  energy  F0.  For  example,  using  a  1.5-MeV. 
35-kA  electron  beam,  protons  were  accelerated  to  energies 
of  greater  than  8  MeV  (5.3  Fn/nucleon)  and  xenon  ions 
were  accelerated  to  energies  of  greater  than  638  MeV  (4.9 
Fr/nucieon )  b  It  was  also  found  that  the  maximum  ion  ener¬ 
gy  was  essentially  independent  of  the  applied  magnetic 
field.7 

The  space-charge  limiting  current  for  a  system  is  the 
maximum  current  that  can  be  propagated  through  the  sys¬ 
tem  before  the  potential  depression  produced  by  the  space- 
charge  exceeds  the  beam  energy.  When  the  injected  current 
exceeds  the  space-charge  limiting  current,  a  virtual  cathode 
forms  where  the  potential  equals  the  beam  energy.  The  frac¬ 
tion  of  the  injected  current  which  exceeds  the  space-charge 
limiting  current  is  reflected  by  the  virtual  cathode  back 
through  the  anode  or  deflected  to  the  drift  tube  wall.  For  a 
solid  beam  of  radius  Rh  in  a  cylindrical  drift  tube  of  radius 
R„  with  a  large  focusing  magnetic  field,  the  space-charge 
limiting  current  is  given  by  the  formula'1 


4 


HOOPOE1-  D1  : 
1  -t-  21n (R„/Rh) 


—  (Amps), 


(1) 


where  is  the  Lorentz  factor  for  the  injected  electrons  and 
where  it  is  assumed  that  Rh 

Many  explanations  have  been  offered  for  collective  ac¬ 
celeration  in  both  gas-filled  systems  and  evacuated  systems. 
Among  the  first  models  to  be  proposed  for  collective  ion 
acceleration  in  a  gas-filled  drift  tube  was  that  of  Olson.4  Ac¬ 
cording  to  this  model,  which  is  based  on  work  by  Poukev  and 
Rostoker, 10  collective  ion  acceleration  occurs  only  whe.:  the 
current  exceeds  the  space-charge  limiting  current.  Soon 
after  the  electron  beam  is  injected,  the  beam  front  stops  near 
the  anode,  forming  a  stationary  potential  well  with  a  depth  of 
two  to  three  times  the  electron  beam  energy.  As  the  space 
charge  near  the  anode  is  neutralized  by  ionization  of  the  gas, 
the  beam  front  and  the  potential  weii  suddenly  start  to  prop- 
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FIG  1  Schematic  of  experimental  setup  for  University  of  Maryland  codec 
me  acceleration  experiments 


agate  downstream.  Ions  are  accelerated  by  the  deep,  station¬ 
ary  potential  well  and  by  the  moving  potential  well.  This 
theory  accounts  for  many  of  the  experimental  observations 
of  collective  ton  acceleration  in  gas-filled  systems  and  is  con¬ 
sistent  with  some  2-D  particle  simulations  of  IREB  injection 
into  a  neutral  gas." 

The  energies  that  are  typically  observed  in  evacuated 
systems  such  as  the  Luce  diode  configuration  are  much  too 
high  to  be  explained  by  the  theory  of  Olson.  Furthermore, 
Adier,  Nation,  and  Serlin12  have  also  reported  detailed  ob¬ 
servations  of  collective  ion  acceleration  in  a  Luce  diode  con¬ 
figuration  which  indicate  that  ion  acceleration  was  not  pro¬ 
duced  by  a  deep  potential  well  or  by  a  moving  potential  well 
at  the  beam  front.  They  suggested  instead  that  ions  were 
accelerated  by  electrostatic  waves  which  are  driven  by  ihe 
two-stream  instability.  Destler  et  al  *  suggested  that  in  evac¬ 
uated  systems  collective  acceleration  occurs  by  a  bootstrap 
process.  The  early  part  of  the  beam  pulse  forms  a  dense  plas¬ 
ma  near  the  anode.  If  the  beam  current  exceeds  the  space- 
cnarge  limiting  current  and  the  plasma  density  is  sufficiently 
high,  a  virtual  cathode  forms  at  the  downstream  edge  of  the 
plasma.  The  intense  virtual  cathode  electric  fields  draw  ions 
out  of  the  plasma,  the  ions  neutralize  the  electron  beam 
space  charge,  and  the  virtual  cathode  moves  further  down¬ 
stream.  As  this  process  is  repeated,  the  beam  electrons  trans¬ 
fer  momentum  to  the  ions  until  the  ions  attain  the  same  ve¬ 
locity  as  the  beam-front  electrons.  A  similar  model  was 
proposed  earlier  by  Ryutov  and  Stupakov1  ’  for  collective 
ion  acceleration  produced  by  injection  of  an  IREB  through  a 
dense  plasma  layer.  Mako  and  TajimaM  modified  the  model 
using  data  obtained  from  experiments  with  a  gas-filled  drift 
tube,  and  were  able  to  recover  the  experimentally  measured 
peak  ion  energy  (ten  times  the  electron  beam  energy)  and 
ion  distribution.  They  also  used  a  particle  code  to  simulate 
injection  of  an  IREB  into  an  evacuated  drift  tube  with  a 
dense  plasma  layer  located  next  to  the  injection  plane.  The 
simulations  showed  that  the  peak  ion  energies  were  limited 
by  phase  instability  of  the  ions  and  the  accelerating  electric 
fields.  However,  using  an  artificially  low  ion  mass,  they  also 
found  that  the  peak  ion  energies  never  exceeded  three  times 
the  electron  beam  energy  and  that,  when  the  ion  mass  was 
doubled,  the  ion  velocity  decreased  inversely  as  the  square 
root  of  the  mass.  Using  a  2-1 /2-D  electromagnetic  particle 
code  with  a  realistic  ion  mass,  Faehl  and  Peter"  and  Peter  er 
a/.1'1  simulated  the  same  system  and  found  that  the  virtual 


cathode  formed  at  the  downstream  edge  of  the  plasma  layer 
soon  after  the  beam  was  injected  and  that  no  subsequent 
motion  of  the  virtual  cathode  occurred.  They  also  pointed 
out  that  no  such  motion  should  be  expected  because  the 
maximum  ion  current  that  can  be  extracted  from  the  plasma 
is  never  large  enough  to  neutralize  the  electron  beam  space 
charge.  They  suggested  that  collective  acceleration  is  the  re¬ 
sult  of  ambipolar  expansion  of  the  plasma  rather  than  mo¬ 
tion  of  the  virtual  cathode. 

In  simulations  of  injection  of  an  IREB  from  a  plasma 
into  free  space,  Chang  and  Reiser1  found  that  as  ions  were 
extracted  from  the  surface  of  the  plasma  by  the  virtual  cath¬ 
ode  electric  fields,  the  location  of  the  virtual  cathode  mini¬ 
mum  shifted  downstream  from  its  initial  location.  They  also 
observed  ions  which  had  been  accelerated  to  energies  ap¬ 
proximately  twice  the  electron  beam  energy  and  proposed 
that  the  coherent  motion  of  the  virtual  cathode  and  the  ions 
as  the  virtual  cathode  shifts  downstream  to  its  final  position 
was  responsible  for  accelerating  the  ions. 

In  this  paper  we  present  a  detailed  picture  of  the  accel¬ 
eration  mechanism  in  the  University  of  Maryland  experi¬ 
ments,  in  which  an  intense  relativistic  electron  beam  is  in¬ 
jected  into  an  evacuated  drift  tube  with  a  localized  gas  cloud 
located  near  the  anode.  We  have  obtained  this  picture  using 
particle-in-cell  (PIC)  code  simulations  which  realistically 
model  ionization  of  the  neutral  gas. 

Our  simulations  differ  from  previous  simulations  of  col¬ 
lective  acceleration  in  the  following  respects.  In  our  simula¬ 
tions,  an  IREB  is  injected  into  an  evacuated  drift  tube  with  a 
well-localized  gas  cloud  located  next  to  the  anode,  as  op¬ 
posed  to  a  gas-filled  drift  tube  or  an  evacuated  drift  tube  with 
a  plasma  layer  located  next  to  the  anode.  Ions  are  produced 
in  the  drift  tube  by  realistically  modeling  ionization  of  the 
gas  by  impact  ionization,  as  opposed  to  assuming  an  initial 
distribution  of  ions  or  injecting  ions.  Realistic  ion  masses  are 
also  used.  The  finite  radius  of  the  beam  is  also  taken  into 
account,  so  that  the  simulations  are  more  realistic  than 
strictly  1-D  simulations.  Finally,  the  IREB  is  injected  into  a 
grounded  drift  tube  rather  than  free  space  so  that  inductive 
effects  can  confidently  be  neglected.  (We  also  note  that  our 
code  is  simple  enough  to  permit  a  detailed  investigation  of 
the  acceleration  mechanism. ) 

In  our  simulations,  we  find  that  when  the  injected  elec¬ 
tron  current  exceeds  the  space-charge  limiting  current,  a  vir¬ 
tual  cathode  forms  near  the  anode  plane  soon  after  the  beam 
is  injected.  Neutralization  of  the  electron  beam  space  charge 
by  impact  ionization  of  the  neutral  gas  then  causes  the  vir¬ 
tual  cathode  to  move  downstream  to  the  end  of  the  gas  re¬ 
gion.  where  the  virtual  cathode  stops.  As  the  virtual  cathode 
moves  through  the  gas.  ions  are  accelerated  to  energies  sev¬ 
eral  times  the  electron  beam  energy  by  the  coherent  motion 
of  the  ions  and  the  intense  virtual  cathode  electric  fields 
( This  picture  is  similar  to  the  mechanism  observed  by  Chang 
and  Reiser. ) 

Before  proceeding,  we  would  like  to  point  out  that  the 
general  features  of  the  above  picture  are  similar  to  results 
which  have  been  obtained  earlier  at  the  University  of  Mary¬ 
land"'  The  present  work  has  also  been  extended  to  simulate 
enhanced  collective  ion  acceleration  in  a  laser-controlled 
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FIG  2.  Schematic  of  simulation  model 


collective  accelerator.19  Applications  of  the  system  to  in¬ 
tense  beam  propagation  are  summarized  elsewhere.20 

In  the  following  sections,  we  first  describe  the  computer 
code  which  we  have  written  to  simulate  collective  ion  accel¬ 
eration  in  an  intense  electron  beam-localized  gas  cloud  sys¬ 
tem.  We  then  discuss  the  acceleration  mechanism  which  we 
observe  in  the  simulations.  Finally,  we  discuss  the  depend¬ 
ence  of  the  peak  ion  energy  on  the  system  parameters. 

II.  DESCRIPTION  OF  COMPUTER  COOE 

In  our  simulations,  a  beam  of  current  /„,  energy  K,  and 
radius  Rh  is  injected  into  a  grounded  cylindrical  drift  tube  of 
length  d  and  radius  Ru.  along  the  axis  of  the  drift  tube.  A 
schematic  of  the  model  is  shown  in  Fig.  2.  The  region  of  the 
drift  tube  extending  from  the  anode  to  a  distance  z0  down¬ 
stream  is  filled  by  neutral  hydrogen  at  a  constant  pressure p0. 
The  electron  beam  is  assumed  to  be  focused  by  an  infinitely 
strong  guide  magnetic  field,  so  that  particles  in  the  simula¬ 
tion  move  only  along  the  axis  of  the  drift  tube.  The  beam 
radius  is  also  assumed  to  be  much  smaller  than  the  wall  radi¬ 
us  so  that  the  charge  and  current  density  and  the  axial  elec¬ 
tric  field  are  approximately  uniform  across  the  beam  cross- 
section. 

The  macroparticles  in  the  simulation  obey  the  relativis¬ 
tic  equations  of  motion 

—  =  —  (2) 

dt  m 

dz  , ,, 

dt 

where  z  is  the  axial  position  of  the  center  of  the  macroparti¬ 
cle,  u  is  the  four-velocity  of  the  macroparticle, 
v  =  u/\  1  (u/cV ,  q  is  the  charge  and  m  is  the  mass  of  the 
macroparticle,  and  E.  is  the  electric  field  on  the  axis  of  the 
drift  tube  at  z.  The  code  integrates  the  equations  of  motion 
forward  in  time  using  a  leapfrog  scheme. 

The  electric  field  Ez  on  axis  computed  by  first  comput¬ 
ing  the  potential  6  on  axis  and  then  numerically  computing 
the  derivative  E.  =  —  d<t>/dz.  The  potential  d>  is  computed 
by  finite  differencing  Poisson's  equation  in  z.  taking  the  finite 
sine  transform  of  the  finite-differenced  equation,  analytical¬ 
ly  solving  the  resulting  system  of  ordinary  differential  equa¬ 
tions  in  r.  and  taking  the  inverse  finite  sine  transform. 

Ionization  of  the  neutral  gas  is  modeled  by  monitoring 


the  amount  of  ionization  in  each  grid  cell  that  is  produced  by 
impact  ionization.  When  the  total  number  of  ions  produced 
in  a  grid  cell  exceeds  the  number  of  ions  in  an  ion  macropar¬ 
ticle,  an  ion  macroparticle  and  an  electron  macroparticle  are 
created  at  the  center  of  the  grid  cell.  The  total  ionization 
produced  in  a  grid  cell  in  one  time  step  is  computed  using  the 
formula 

AN=2,N,S,(KE,)v,p0At,  (4) 

where  the  sum  is  taken  over  all  macroparticles  in  the  grid 
cell.  A',  is  the  number  of  real  particles  in  the  ith  macroparti¬ 
cle,  S,  is  the  ionization  cross  section  in  events/m/Torr  for 
each  particle  in  the  ith  macroparticlc,  KE.  is  the  kinetic  en¬ 
ergy  of  each  particle  in  the  ith  macroparticle,  v ,  is  the  veloc¬ 
ity  of  the  macroparticle  in  m/s,  p0  is  the  pressure  of  the  neu¬ 
tral  gas  in  Toit  at  the  center  of  the  grid  cell,  and  At  is  the  time 
step  in  seconds.  The  ionization  cross  sections  S,  used  in  the 
simulation  are  piecewise-linear  approximations  of  published 
ionization  cross  sections  for  molecular  hydrogen.9  In  units 
of  events/cm/T orr,  the  cross  section  for  electron  impact  ion¬ 
ization  used  in  the  simulation  is  given  by 

'0,  KE  <  20  eV, 

_  ylO,  20  eV<AT£  <  yTOx  1CP  eV, 

10 >/KE,  \J0X  10:  t\<.KE<  105  eV, 

.10- 2.  105  eV <KE. 

where  KE  is  the  electron  kinetic  energy  in  eV.  The  cross 
section  for  proton  impact  ionization  used  in  the  simulation  is 
given  by 


fo. 


5,  = 


exp(  2.383  log,n  KE  -  8.5771, 
exp(  -  1.839  log,,,  KE  +  11.036), 


KE  <  10  eV, 

10  eV<AT£ <  4.3  X  104  eV, 
4.3  x  104  eV kKE, 


where  KE  is  the  proton  kinetic  energy  in  eV. 

In  the  simulations,  we  chose  to  run  at  currents  less  than 
35  kA.  and  gas  cloud  pressures  less  than  600  mTorr  because 
we  found  that  at  high  currents  and  high  pressures  and  long 
pulse  lengths,  the  number  of  particles  in  the  simulation  be¬ 
came  too  large 

Finally,  we  note  that  when  no  gas  is  present  the  current 


that  is  observed  downstream  agrees  with  the  limiting  current 
formula  Eq.  (  1 ). 

III.  PICTURE  OF  COLLECTIVE  ACCELERATION 
MECHANISM 

We  have  run  the  code  for  many  different  values  of  the 
electron  beam  current,  electron  beam  energy,  drift  tube  radi- 


1652 


J.  Appi  Phys  ,  Vol.  67.  No.  4,  15  February  1990 


R  L  YaoandC  D  Strtffler 


1652 


us.  drift  tube  length,  neutral  gas  pressure,  and  gas  cloud 
width.  We  present  a  picture  of  the  collective  acceleration 
mechanism  that  we  have  observed  in  all  of  the  simulations  in 
which  collective  acceleration  occurred.  To  do  so,  we  will  use 
the  results  of  a  single  simulation.  (We  note  that  the  param¬ 
eter  values  differ  from  typical  parameter  values  in  the  Uni¬ 
versity  of  Maryland  experiments.  However,  simulations  per¬ 
formed  with  the  experimental  parameter  values  show  that 
the  acceleration  mechanism  is  the  same  as  described  here. ) 

In  the  simulation,  a  15-kA,  500-keV,  0.5-cm-radius 
electron  beam  is  injected  into  a  4-cm-radius,  40-cm-long 
drift  tube  with  a  1-cm-wide,  450-mTorr  hydrogen  cloud  lo¬ 
cated  next  to  the  anode  plane.  The  beam  current  and  beam 
energy  are  assumed  to  rise  instantaneously.  According  to 
Eq.  ( 1 ),  the  space-charge  limiting  current  for  the  system  is 
approximately  1.4  kA.  The  time  step  is  approximately  0.75 
ps  and  the  grid  cell  width  is  approximately  0.2  mm. 

Figure  3  shows  the  time-integrated  proton  energy  spec¬ 
trum  7.5  cm  from  the  anode  plane  1 1  ns  after  the  start  of  the 
run.  The  energy  spectrum  can  be  roughly  divided  into  two 
parts — the  bulk  of  the  spectrum,  which  is  centered  about  0.6 
MeV,  and  a  tail  that  extends  out  to  1 .5  MeV  ( three  times  the 
electron  beam  energy ) .  We  ignore  the  ions  in  the  bulk  of  the 
spectrum  and  focus  on  the  ions  in  the  high-energy  tail.  Note 
that  since  the  potential  depression  at  the  center  of  the  drift 
tube  is  approximately  —  500  kV,  the  ion  spectrum  that 
would  be  observed  by  a  grounded  detector  would  be  the 
same  as  the  spectrum  shown  in  Fig.  3  shifted  t  a  the  left  by  0.5 
MeV.  In  the  following,  we  present  a  detailed  picture  showing 
how  the  ions  in  the  high-energy  tail  are  accelerated  to  ener¬ 
gies  exceeding  the  electron  beam  energy. 

In  Fig.  4,  the  potential  6  on  axis  near  the  anode  plane  is 
pii'Uvd  as  a  function  of  distance  z  from  the  anode  plane  at 
several  different  time  steps.  About  40  ps  after  the  start  of  the 
run.  a  virtual  cathode  forms  about  2  mm  from  the  anode 
plane,  and  the  current  measured  downstream  from  the  vir¬ 
tual  cathode  falls  to  a  value  approximately  equal  to  the 
space-charge  limiting  current.  After  formation  of  the  virtual 
cathode,  the  potential  remains  essentially  the  same  for  sever¬ 
al  thousand  time  steps  ( the  position  and  amplitude  of  the 
virtual  cathode  oscillate  in  time  at  the  local  plasma  frequen- 
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cy,  but  the  time-averaged  position  and  amplitude  are  ap¬ 
proximately  the  same).  The  potential  for  time  step  4600. 
which  is  shown  in  Fig.  4,  is  typical  for  this  stage.  During  this 
stage,  the  beam  electrons  ionize  the  neutral  gas,  producing  a 
plasma  in  the  vicinity  of  the  virtual  cathode.  The  intense 
electric  fields  near  the  virtual  cathode  quickly  expel  the  plas¬ 
ma  electrons  from  the  plasma  region.  The  plasma  ions  are 
also  accelerated  by  these  electric  fields  but  leave  the  plasma 
region  at  a  much  slower  race  than  the  plasma  electrons.  They 
also  tend  to  be  accelerated  downstream,  whereas  the  plasma 
electrons  tend  to  be  accelerated  upstream,  where  they  exit 
from  the  drift  tube  through  the  anode  plane.  As  plasma  ions 
drift  through  the  gas.  they  also  contribute  *o  the  ionization. 
As  the  ion  density  in  the  plasma  region  increases,  the  virtual 
cathode  moves  downstream  to  the  end  of  the  plasma  region 
leaving  behind  a  positive  potential  region  ( time  steps  4600- 
5000 ) .  Ions  lying  to  the  left  of  the  peak  of  the  positive  poten¬ 
tial  are  accelerated  upstream  and  exit  the  drift  tube  through 
the  anode  plane.  Ions  lying  to  the  right  of  the  peak  of  the 
positive  potential  are  accelerated  downstream.  As  the  ions 
move  downstream,  the  virtual  cathode  also  moves  down¬ 
stream  with  the  positive  potential  region  widening  and  ex¬ 
tending  downstream  behind  it  ( time  steps  5000-5600 ) .  The 
virtual  cathode  ceases  to  move  when  it  reaches  the  end  of  the 
gas  region  at  z  —  1  cm  ( time  step  5600 ) .  We  will  see  that  ions 
are  accelerated  to  energies  of  nearly  1.5  MeV  dunng  the  time 
that  the  virtual  cathode  and  the  positive  potential  region 
move  downstream. 

The  formation  of  a  positive  potential  region  has  been 
observed  expenmentally  and  in  other  simulations. 14  The 
reason  the  potential  in  the  plasma  region  becomes  positive 
can  be  understood  as  follows.  Because  the  electron  impact 
ionization  cross  section  peaks  at  low  energies,  the  majority  of 
ions  are  formed  near  the  virtual  cathode.  The  buildup  of  ions 
in  the  vicinity  of  the  virtual  cathode  causes  the  beam  elec¬ 
trons  to  spend  more  time  there,  and  as  a  result  the  electron 
density  downstream  from  the  virtual  cathode  decreases. 
Hence,  when  the  ion  density  is  sufficiently  high,  the  position 
of  the  virtual  cathode  can  no  longer  be  maintained  and  the 
virtual  cathode  begins  to  move  downstream.  (In  fact,  our 
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base  of  the  trajectory  indicate  the  location  of  the  peak  electric  field  at  the 
same  time  steps  See  Fig  3  for  system  parameters 


simulations  show  that  immediately  preceding  the  onset  of 
the  acceleration  process,  the  electron  current  measured  im¬ 
mediately  downstream  from  the  plasma  region  falls  abrupt¬ 
ly.  )  As  the  virtual  cathode  moves  through  the  plasma  re¬ 
gion,  the  electrons  behind  the  virtual  cathode  travel  through 
the  plasma  with  essentially  the  same  speed  with  which  they 
are  injected.  The  increase  in  the  speed  of  the  electrons 
through  the  plasma  region  produces  a  sudden  drop  in  the 
electron  density.  Hence,  the  virtual  cathode  leaves  a  net  posi¬ 
tive  space  charge  in  its  wake  and  the  potential  becomes  posi¬ 
tive  in  the  plasma  region. 

In  Fig.  5.  the  phase-space  trajectory  of  a  collectively 
accelerated  proton  is  shown  during  the  collective  accelera¬ 
tion  process.  The  ion  is  created  with  zero  velocity  at  a  dis¬ 
tance  of  about  1.3  mm  from  the  anode  plane  about  2.6  ns 
after  the  start  of  the  run  (time  step  3500).  It  is  rapidly  accel¬ 
erated  by  the  electric  field  near  the  virtual  cathode  to  a  veloc¬ 
ity  of  about  0.008c  ( corresponding  to  an  energy  of  30  keV ) 
and  begins  to  drift  downstream.  When  it  has  drifted  about  2 
mm  downstream,  the  potential  begins  to  go  positive  (at 
around  time  step  4700)  and  the  virtual  cathode  begins  to 
move  downstream.  As  the  virtual  cathode  moves  down¬ 
stream  and  the  positive  region  extends  downstream  behind  it 
( time  steps  4800-5600 ) .  the  proton  is  rapidly  accelerated  to 
a  final  velocity  of  about  0.054c  ( corresponding  to  an  energy 
of  1.4  MeV). 

The  boxes  plotted  along  the  proton  phase-space  trajec¬ 
tory  indicate  points  along  the  trajectory  which  are  separated 
by  equal  time  intervals  of  200  time  steps  or  about  150  ps. 
Each  box  is  labeled  with  the  corresponding  time  step  starting 
with  time  step  4800.  Below  the  phase-space  trajectory,  we 
have  also  plotted  the  location  of  the  peak  electric  field  at  the 
same  time  steps.  The  peak  electric  field  is  always  located  in 
the  region  where  the  potential  drops  sharply  from  the  edge  of 
the  positive  potential  region  to  the  virtual  cathode.  As  the 
virtual  cathode  moves  downstream  and  the  positive  poten¬ 
tial  region  extends  downstream  behind  the  virtual  cathode 
the  peak  electric  field  also  moves  downstream.  As  can  be 
seen  from  the  figure,  the  ton  follows  closely  behind  the  peak 


electric  field  for  time  steps  4800  to  5600.  It  is  this  coherent 
motion  of  the  ion  and  the  peak  electric  field  through  the 
plasma  region  which  enables  the  ion  to  be  accelerated  to 
energies  several  times  the  energy  it  would  gain  if  the  poten¬ 
tial  were  stationary. 

Note  that  the  trajectory  of  the  proton  in  Fig.  5  is  similar 
to  the  trajectory  of  an  ion  being  accelerated  by  a  constant 
electric  field  v. (z)  =  v  -t-  2 (eEJmp)z,  where  v#  is  the 
initial  ion  velocity.  Using  values  obtained  from  the  figure  in 
the  above  formula,  we  find  that  the  electric  field  E.  is  about 
200  MV/m,  which  is  near  the  peak  virtual  cathode  electric 
field  of  250  MV/m. 

The  reason  for  the  coherent  motion  of  the  ion  and  the 
peak  electric  field  during  the  motion  of  the  potential  front 
can  be  understood  by  considering  the  free  expansion  of  a  slab 
of  ions,  which  initially  has  uniform  density.  As  the  slab  ex¬ 
pands,  an  ion  at  the  surface  of  the  slab  always  feels  the  same 
electric  field,  that  is.  the  ion  and  the  electric  field  travel  at  the 
same  velocity.  We  believe  that  the  same  effect  occurs  during 
motion  of  the  potential  front.  In  the  simulations,  the  ions  at 
the  front  of  the  accelerated  ion  beam  are  accelerated  not  only 
by  the  space  charge  produced  by  the  ions  following  behind 
but  also  by  the  electric  field  produced  by  the  virtual  cathode 
which  always  moves  in  front  of  the  ion  beam — the  ions  are 
pushed  ahead  by  the  ions  which  follow  and  simultaneously 
pulled  ahead  by  the  virtual  cathode. 

The  reason  the  potential  from  stops  moving  at  the  end  of 
the  gas  region  can  also  be  understood  by  considering  the  free 
expansion  of  a  slab  of  ions,  as  the  slab  expands,  the  ion 
density  and  the  potential  decrease  uniformly.  In  the  gas  re¬ 
gion.  ionization  of  the  gas,  particularly  by  ions,  can  supply 
ions  at  a  rate  sufficiently  high  to  prevent  total  collapse  of  the 
potential  and  reformation  of  a  virtual  cathode  near  the  an¬ 
ode.  Beyond  the  gas  region,  no  ions  are  produced  and  the  ion 
density  is  too  low  to  permit  the  virtual  cathode  to  move  for¬ 
ward.  (As  Faehl  and  Peter15  observed,  the  maximum  ion 
current  which  can  be  extracted  from  the  plasma  is  space- 
charge  limited  and  hence  ions  can  never  be  extracted  at  a  rate 
sufficient  to  significantly  neutralize  the  electron  beam  space 
charge  downstream. ) 

Although  all  of  the  collectively  accelerated  ions  are  ac¬ 
celerated  by  the  moving  potential  from,  some  ions  continue 
to  accelerate,  although  at  a  much  lower  rate,  after  they  have 
left  the  gas  region.  This  additional  acceleration  is  probably 
produced  by  the  interaction  of  the  collectively  accelerated 
ions  among  themselves — interactions  between  tons  that 
have  been  collectively  accelerated  in  a  bunch  will  cause  ions 
with  velocities  greater  than  the  mean  velocity  to  be  acceler¬ 
ated  and  ions  with  velocities  less  than  the  mean  velocity  to  be 
decelerated. 

Several  bursts  of  accelerated  ions,  which  follow  one  an¬ 
other  at  short  intervals,  are  often  observed  in  simulations. 
The  multiple  bursts  are  produced  when  the  positive  potential 
region  collapses  after  the  tons  have  been  accelerated  out  of 
the  gas  region,  the  virtual  cathode  retreats  to  a  location  in¬ 
side  the  gas  region  and  the  acceleration  process  repeats  itself. 
If  the  pressure  is  increased,  with  all  other  system  parameters 
kept  constant,  the  multiple  bursts  are  suppressed,  probably 
because  the  rate  of  ion  production  in  the  gas  region  is  suffi- 
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cient  to  prevent  collapse  of  the  positive  potential  region. 
Such  multipie  bursts  have  been  observed  in  some  experi¬ 
ments  in  evacuated  systems. 12 

IV.  RESULTS  OF  PARAMETER  STUDIES 
A.  Dependence  on  beam  energy  and  current 

In  Fig.  6,  the  peak  proton  energy  Em,%  measured  at 
2  =  7.5  cm  in  the  first  10  ns  of  the  run  is  plotted  as  a  function 
of  the  injected  beam  current  70  for  beam  energies  Va  of  0.5, 
1.0,  1.5,  and  3.0  MeV.  In  each  case,  the  gas  pressure  is  600 
mTorr  and  the  gas  cloud  width  is  1  cm.  Generally,  £m4J 
increases  with  beam  current  and  beam  energy.  In  Fig.  7,  the 
same  data  are  plotted  with  £mis  normalized  to  Vn  and  70 
normalized  to  the  space-charge  limiting  current  IL ,  given  by 
Eq.  (1).  Notice  that  all  the  plotted  points  lie  on  approxi¬ 
mately  the  same  curve.  For  all  values  of  V0  which  were  used, 
no  collectively  accelerated  protons  are  observed  for 
70 /IL  <  1.  As  Iq/Il  increases  above  unity,  £m„  increases 
sharply  and  then  continues  to  increase  at  at  a  much  slower 
rate.  (Note  that  for  fixed  /o.  the  ratio  ■  JIL  decreases  with 
increasing  V0.  Hence  for  the  same  current  70,  the  ratio 
£m«»/*o  wth  V0  =  0.5  MeV  is  much  larger  than  the  ratio 
for  F0  =  3  MeV.) 

The  observed  dependence  of  the  peak  proton  energy  on 
beam  energy  and  current  can  be  understood  quite  simply  by 
recalling  that  the  strength  of  the  peak  electric  field  increases 
with  F0and7(/7t .  In  Fig.  8(a),  the  phase-space  trajectories 
of  collectively  accelerated  protons  are  shown  for  three  differ¬ 
ent  values  of  the  injected  current — 10,  20,  and  30  kA — and 
the  same  beam  energy  of  0.5  MeV.  In  each  case,  the  gas 
pressure  is  600  mTorr  and  the  gas  cloud  width  is  1  cm.  Re¬ 
calling  that  the  accelerating  electric  field  is  approximately 
proportional  to  the  product  of  the  velocity  „  cr.d  the  dope 
du/dz  of  the  trajectory,  it  can  be  seen  that  as  the  current 
increases  so  does  the  accelerating  electric  fielF  In  Fig.  8(b), 
the  phase-space  trajectories  of  collectively  derated  pro¬ 
tons  are  shown  for  three  different  values  of  the  beam  ener¬ 
gy — 0.5.  1.0,  and  1.5  MeV — and  the  same  beam  current  of 
30  kA.  It  can  be  seen  that  as  the  beam  energy  increases  so 
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current  l0  normalized  to  space-charge  limiting  current  /,  .  See  Figs.  3  and  6 
for  system  parameters. 


does  the  accelerating  electric  field.  (This  dependence  is  not 
so  clear  when  the  trajectories  for  0.5  and  1.0  MeV  are  com¬ 
pared,  which  illustrates  the  complex  nature  of  the  accelera¬ 
tion  process,  that  is,  the  final  ion  energy  depends  not  only  on 
the  strength  of  the  electric  field  but  also  on  the  location  and 
velocity  of  the  ion  when  the  virtual  cathode  begins  to  move 
downstream  and  on  the  velocity  with  which  the  virtual  cath¬ 
ode  moves. ) 
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FIG  8  Phase-space  trajectories  of  collectively  accelerated  protons  for  l  a  l 
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B.  Dependence  on  gas  cloud  pressure 

In  Fig  9.  the  peak  proton  energy  £ma,  measured  at 
c  =  ~.5  cm  in  the  first  20  ns  of  the  beam  pulse  is  plotted  as  a 
function  of  gas  pressure  p„  for  injected  beam  currents  of 
20.  25.  and  30  kA  and  for  beam  energies  I’,  of  (a)  0.5.  (b) 

1  0.  and  ( c  )  15  MeV.  In  each  case,  the  gas  cloud  width  is  1 .0 
cm.  At  50  mTorr.  no  collectively  accelerated  protons  are 
observed  in  the  first  20  ns  of  the  beam  pulse.  Above  50 
mTorr.  £mj<  tends  to  increase  with  increasing  pressure  and 
level  off  at  energies  several  times  the  electron  beam  energy 
Generally,  the  higher  the  beam  energy,  the  higher  the  energy 
and  pressure  at  which  the  energy  levels  off 


The  observed  dependence  of  the  peak  proton  energy  on 
the  pressure  can  be  explained  qualitatively  as  follows.  At  low 
pressures,  the  amount  of  ionization  at  the  end  of  20  ns  is  too 
small  to  neutralize  the  electron  beam  space  charge  and  the 
virtual  cathode  remains  at  its  initial  position;  no  collective 
acceleration  occurs  and  the  peak  proton  energies  which  are 
measured  are  less  than  the  electron  beam  energy.  As  the 
pressure  is  increased,  the  amount  of  ionization  which  is  pro¬ 
duced  during  the  beam  pulse  becomes  large  enough  to  neu¬ 
tralize  the  electron  beam  space  charge  and  permit  the  virtual 
cathode  to  move  downstream.  As  it  moves  downstream,  the 
virtual  cathode  electnc  fields  accelerate  ions  to  energies  sev¬ 
eral  times  the  electron  beam  energy.  Once  the  acceleration 
process  has  begun,  the  final  energy  of  a  collectively  acceler¬ 
ated  ion  depends  only  on  the  magnitude  of  the  accelerating 
electric  field  and  the  distance  over  which  the  ion  travels  co¬ 
herently  with  the  electric  field.  Now,  the  magnitude  of  the 
electric  field  will  be  determined  mainly  by  the  electron  beam 
energy  and  current  and  the  density  of  ions  in  the  plasma 
region  at  the  start  of  the  acceleration  process  (because  the 
electric  field  at  the  surface  of  a  slab  of  ions  is  proportional  to 
the  ion  density  in  the  slab ) .  But  the  ion  density  at  the  start  of 
the  acceleration  process  is  independent  of  the  pressure,  be¬ 
cause  it  is  proportional  to  the  electron  beam  density,  which  is 
determined  by  the  electron  beam  energy  and  current  only. 
The  only  effect  of  increasing  the  gas  pressure  will  be  to  de¬ 
crease  the  time  required  for  the  ion  density  to  build  up  before 
the  acceleration  process  begins  (which  is  in  fact  observed  in 
simulations ) .  Hence,  the  peak  ion  energy  will  depend  only 
weakly  on  the  neutral  gas  pressure,  once  the  gas  pressure  is 
sufficiently  high  to  cause  collective  acceleration. 

The  above  picture  is  somewhat  too  simplistic  because 
the  gas  pressure  also  affects  the  motion  of  the  virtual  cathode 
electric  fields  through  the  gas  cloud.  Generally,  as  the  pres¬ 
sure  increases,  the  virtual  cathode  moves  through  the  gas 
cloud  more  quickly  and  higher  peak  proton  energies  are  ob¬ 
served  because  the  virtual  cathode  electnc  fields  are  able  to 
accelerate  the  protons  over  a  longer  distance. 

Although  the  dependence  of  collective  acceleration  on 
pressure  as  observed  in  our  simulations  is  consistent  with  the 
model  for  collective  acceleration  which  we  have  presented 
above,  it  disagrees  with  the  expenmental  results  obtained  at 
the  University  of  Maryland.  In  particular.  Floyd  et  alh  re¬ 
ported  that  when  a  1.5-MeV.  35-kA.  30-ns  beam  pulse  was 
injected  into  a  25-cm-diam  dnft  tube  with  a  well-localized 
gas  cloud  located  near  the  anode,  collective  acceleration  was 
observed  only  in  a  narrow  pressure  window.  In  our  simula¬ 
tions.  we  observe  no  such  pressure  window.  In  particular,  we 
observe  collective  acceleration  at  all  pressures  up  to  1  Torr 

The  disagreement  may  have  the  following  causes.  First, 
in  our  simulations  we  have  neglected  charge  exchange  be¬ 
tween  protons  and  H:  molecules.  Since  the  mean  free  path 
for  charge  exchange  for  protons  in  H  -  is  approximately  1  cm 
at  100  mTorr  at  approximately  100  keV.1'  at  pressures  of  a 
few  1  (X)  m  1  orr  charge  exchange  probably  prevents  accelera¬ 
tion  of  protons  to  energies  much  greater  than  a  few  100  keV. 
Second,  in  our  simulations  we  have  neglected  sources  of  ioni¬ 
zation  such  as  photoionization  and  prepulse  ionization 
which  could  produce  sufficient  ionization  in  the  gas  region 
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to  prevent  formation  of  a  vtrtuai  cathode  near  the  injection 
plane.  In  simulations  of  injection  of  an  IREB  into  vacuum 
with  a  localized  plasma  located  next  to  the  injection 
plane, 15  '6-21  it  was  found  that  for  plasma  densities  compara¬ 
ble  to  or  exceeding  the  beam  density,  the  virtual  cathode 
formed  near  the  downstream  edge  of  the  plasma.  Finally,  the 
disagreement  with  the  reported  experiments  may  be  the  re¬ 
sult  of  changes  in  the  electron  beam  diode  characteristics  as 
the  pressure  is  increased.  In  the  experiments,  the  peak  gas 
cloud  pressure  at  the  time  of  beam  injection  was  varied  by 
simply  varying  the  time  delay  between  firing  of  the  gas  puff 
valve  and  firing  of  the  diode.  It  was  also  found  that  when  the 
gas  puff  was  fired  too  early,  the  diode  shorted  out.4 


tOO  mTorr 
1 50  mTorr 
200  mTorr 


°  100  mTorr 

5  A  150  mTorr 
^  200  mTorr 


C.  Dependence  on  gas  cloud  width 

In  Fig.  10.  the  peak  proton  energy  £ma,  measured  at 
z  =  7.5  cm  in  the  first  20  ns  of  the  beam  pulse  is  plotted  as  a 
function  of  gas  cloud  width  z0  for  electron  beam  energies  V„ 
of  (a)  0.5.  (b)  1.0,  and  (c)  1.5  MeV  and  several  neutral  gas 
pressures />„.  In  each  case,  the  injected  electron  beam  current 
is  30  kA. 

Figure  10  shows  that  £ma,  tends  to  increase  with  gas 
cloud  width  and  then  level  off.  and  that  this  tendency  be¬ 
comes  more  pronounced  as  the  pressure  increases.  In  Fig 
10(c),  for  example,  which  shows  data  obtained  for  Vn=  1.5 
MeV,  £m>,  increases  approximately  linearly  with  increasing 
gas  cloud  width  before  leveling  off  at  about  6.0  MeV  at 
z„  =  5.0  cm  for p„  =  200  mTorr.  (Whenp,,  =  100mTorr.no 
collectively  accelerated  protons  are  detected  in  the  first  20  ns 
of  the  beam  pulse. ) 

The  tendency  of  £ma,  to  increase  with  increasing  gas 
cloud  width  is  consistent  with  the  simple  picture  of  collective 
acceleration  presented  above.  Protons  are  accelerated  to 
higher  energies  as  the  gas  cloud  width  is  increased  because 
they  are  accelerated  over  longer  distances.  This  trend  is 
more  pronounced  at  higher  pressures  because,  as  noted  in 
the  last  section,  higher  pressures  enable  the  virtual  cathode 
to  move  through  the  gas  cloud  more  quickly.  Closer  exami¬ 
nation  of  the  simulation  results  shows  that  £ma>  tends  to 
increase  as  the  gas  cloud  width  increases  because  the  protons 
are  indeed  accelerated  over  longer  distances.  Figure  1 1 
shows  the  phase-space  trajectories  of  two  co'lectively  accel¬ 
erated  ions  for  the  same  beam  energy  and  current  ( 30  kA ) 
and  the  same  gas  cloud  pressure  (200  mTorr)  but  different 
gas  cloud  widths.  It  is  clear  from  the  figure  that  one  proton  is 
accelerated  to  a  higher  energy  because  it  is  accelerated  over  a 
longer  distarrc.  Th;  figure  also  illustrates  the  reason  for  the 
tendency  of  the  peak  proton  energy  to  level  off  as  the  gas 
cloud  width  is  increased — the  ions  tend  to  outrun  the  mov¬ 
ing  potential  front  as  it  moves  through  the  gas  cloud  (seen  in 
the  3-cm  case). 

Results  of  simulations  performed  with  different  values 
of  the  injected  beam  current  /„  indicate  that  the  value  at 
which  £max  levels  off  also  tends  to  increase  with  beam  cur¬ 
rent,  although  the  correlation  is  relatively  weak.  In  Fig.  12. 
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FIG  10  Peak  proton  energy  measure  at  :  =  "  5  cm  after  20  ns  vs  gas 
cloud  width  r„  for /.  —  30kA.p,,  =  100.  150.  and  200  mTorr.  and  f,  =  lai 
05.  (hi  I0.and(c>  I  5  MeV  See  Fig  .1  for  other  system  parameters 


FIG  1 1  Phase-space  iraiectory  of  collectively  accelerated  proton  for 
I  ,  -  1  5  MeV. /„  =  30kA.p„  =  200  mTorr.  andr,,  =  2  and  3  cm.  See  Fig  3 
for  other  system  parameters 
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FIG  12  Peak  proton  energy  £'rn_tx  measured  ai  z  —  ~.5  cm  alter  20  ns  vs  gas 
cloud  width  r,  for  1  =15  MeV.  /?,.  =  2(X)  mTorr.  and  /.,  =  20.  25.  and  30 
kA  See  Fig  5  for  other  system  parameters. 


for  example,  the  peak  proton  energy  measured  at  r  =  7.5  cm 
in  the  first  20  ns  of  the  beam  pulse  is  plotted  as  a  function  of 
gas  cloud  width  fora  beam  energy  of  1.5  MeV.  a  gas  pressure 
of  200  mTorr.  and  injected  beam  currents  of  20.  25  and  30 
kA 

V.  SUMMARY 

In  order  to  investigate  collective  ion  acceleration  ob¬ 
served  in  experiments  in  which  an  intense  relativisMc  elec¬ 
tron-beam  is  injected  into  an  evacuated  drift  tube  with  a 
localized  gas  cloud  located  near  the  anode,  w  e  have  written  a 
relativistic,  electrostatic  PIC  code  that  realistically  models 
ionization  of  the  gas 

Using  the  code,  we  have  developed  a  simple  model  for 
the  collective  acceleration  mechanism  in  these  experiments. 
In  this  model,  ionization  of  the  gas  causes  the  virtual  cathode 
to  move  downstream  through  the  gas  cloud  As  the  virtual 
cathode  moves  through  the  cas  cloud,  ions  are  accelerated  to 
energies  several  times  the  electron-beam  energy  by  coherent 
motion  of  the  ions  and  the  intense  virtual  cathode  electric 
fields 

Using  the  code,  we  have  also  determined  the  depen¬ 
dence  of  the  peak  ion  energy  on  the  system  parameters,  in¬ 
cluding  the  beam  energy  1,,.  beam  current  /  .  the  gas  cloud 
pressure p,,.  and  the  gas  cloud  width  c,,.  We  find  that  the  peak 
energy  £mj,  tends  to  increase  with  increasing  values  of  beam 
energy  U„  and  beam  current  /,,.  In  particular  we  find  that 
£.„,,  /  U,  tends  to  increase  as  l,,/ 1,  .  where  l ,  is  the  space- 
charge  limiting  current  given  by  Eq.  (  1  ).  and  that  the  func¬ 
tional  dependence  is  nearly  the  same  for  a  wide  range  of 
parameters.  We  aiso  find  that  as  the  pressure  is  increased. 
£,.,,  jumps  from  a  value  approximately  equal  to  the  beam 
energy  to  a  value  several  times  the  beam  energy.  and  that 


generally,  the  higher  the  beam  energy,  the  higher  the  pres¬ 
sure  at  which  £max  reaches  the  final  energy.  We  also  find 
that,  generally,  £max  tends  to  increase  with  gas  cloud  width 
z„  and  then  level  off  when  zn  is  a  few  cm  and  that  this  tenden¬ 
cy  is  more  pronounced  for  higher  pressures.  Our  results  indi¬ 
cate  that  the  mechanism  for  collective  acceleration  which  we 
have  observed  in  our  simulations  can  produce  protons  with 
energies  of  at  least  5-6  times  the  electron  beam  energy  when 
the  beam  and  the  gas  cloud  properties  are  optimized. 
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MEASUREMENTS  AND  SIMULATION  OF  CONTROLLED  BEAMFRONT  MOTION 
IN  THE  LASER  CONTROLLED  COLLECTIVE  ACCELERATOR* 


R.L.  Yao.  WAV.  Destlcr.  C  D.  Striffler,  J.  Rodgers,  Z.  Scgalov 
University  of  Maryland,  College  Park,  N1D  20742 


Abstract 


II.  Experiments 


In  the  Laser  Controlled  Collective  Accelerator,  an  intense 
electron  beam  is  injected  at  a  current  above  the  vacuum  space 
charge  limit  into  an  initially  evacuated  drift  tube.  A  plasma 
channel,  produced  by  time-sequenced,  multiple  laser  beam  ion¬ 
ization  of  a  solid  target  on  the  drift  tube  wall,  provides  the  nec¬ 
essary  neutralization  to  allow  for  effective  beam  propagation. 
By  controlling  the  rate  of  production  of  the  plasma  channel 
as  a  function  of  time  down  the  drift  tube,  control  of  the  elec¬ 
tron  beamfront  can  be  achieved.  Recent  experimental  measure¬ 
ments  of  controlled  beamfront  motion  in  this  configuration  are 
presented,  along  with  results  of  ion  acceleration  experiments 
conducted  using  two  different  accelerating  gradients.  These 
results  are  compared  with  numerical  simulations  of  the  system 
in  which  both  controlled  beamfront  motion  and  ion  accelera¬ 
tion  is  observed  consistent  with  both  design  expectations  and 
experimental  results. 

I.  Introduction 

The  Laser  Controlled  Collective  Accelerator  concept1-3  rep¬ 
resents  an  attempt  to  extend  the  promising  results  from  “natu¬ 
rally  occurring”  collective  ion  acceleration  experiments  to  prac¬ 
tical  accelerators  in  which  the  accelerating  gradient  and  dis¬ 
tance  can  be  systematically  controlled.  The  concept  is  sim¬ 
ilar  to  that  employed  in  the  IFA-1  and  IFA-2  experiments 
of  Olson4'5,  although  the  actual  experimental  configuration  is 
quite  different.  The  basic  concept  behind  the  experiment  is 
shown  in  Fig.  1.  An  intense  relativistic  electron  beam  is  in¬ 
jected  through  a  localized  gas  cloud  into  an  evacuated  drift 
tube  at  a  current  well  above  the  vacuum  space  charge  limit.  A 
virtual  cathode  then  forms  immediately  downstream  of  the  in¬ 
jection  point  and  ions  produced  within  the  localized  gas  cloud 
are  accelerated  to  modest  energies  in  a  manner  similar  to  more 
conventional  collective  accelerators.  At  this  point,  a  channel 
of  plasma  is  produced  in  a  time  sequenced  manner  down  the 
drift  tube  by  laser  ionization  of  a  CH2  target  stnp  located  on 
the  drift  tube  wall.  The  time  sequencing  of  the  plasma  channel 
is  achieved  by  dividing  a  Q-switched  ruby  laser  pulse  into  ten 
approximately  equal  energy  beams  and  using  optical  delays  to 
ionize  sequentially  ten  target  spots  equally  spaced  down  the 
drift  tube.  In  this  manner,  the  virtual  cathode  at  the  beam- 
front  can  be  carefully  accelerated  down  the  drift  tube  and  ions 
trapped  by  the  strong  electnc  fields  at  the  virtual  cathode  can 
be  accelerated  to  high  energies  in  a  controlled  manner. 

In  tins  paper  we  present  in  section  II  results  of  experiments 
in  which  controlled  beamfront  motion  has  been  confirmed  for 
two  different  accelerating  gradients.  Results  of  ion  acceleration 
experiments  are  also  presented.  Numerical  simulations  of  the 
experiments  presented  in  section  III  confirm  both  controlled 
beamfront  motion  and  the  controlled  acceleration  of  ions  by  the 
moving  virtual  cathode  over  significant  distances.  Conclusions 
are  drawn  in  section  IV. 
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As  shown  in  Fig.  1,  an  intense  relativistic  electron  u 
(900  keV.  20  kA,  30  ns)  is  emitted  from  a  4  mm  diamett ;  t,.: 
sten  cathode  located  1  cm  upstream  of  a  stainless  steci  ;u„, 

A  14  mm  diameter  hole  in  the  anode  plate  on  axis  uiiou- 
most  all  of  the  beam  current  to  pass  into  a  downstream 
tube  10  cm  in  diameter.  Seed  protons  for  acceleration  ;u<  . 
vided  by  beam  ionization  of  a  localized  gas  cloud  prounr. ,, 
a  fast  gas  puff  valve.  The  Q-switched  ruby  laser  pulse  ,  (,  j 
ns)  is  divided  into  ten  approximately  equal  energy  beam-  .. 
optically  delayed  to  provide  a  time-sequenced  source  of  ,. 
down  the  50  cm  length  of  the  drift  tube.  Design  consi(ict.,i. 
for  the  experiment  and  results  of  tests  of  the  optical  s\->. 
have  been  reported  previously. 2,3 

Five  current  collecting  wall  probes  were  installed  to  n..; 
sure  beam  current  deflected  to  the  drift  tube  wall  at  the  in- 
front  as  a  function  of  time.  These  probes  were  located  at  ai. 
positions  10,  20,  30,  40,  and  50  cm  downstream  of  the  im,... 
plane.  Total  current  reaching  the  downstream  end  of  im  i;: 
tube  was  measured  using  a  Faraday  cup,  and  accelcrmr,.  . 
energies  were  measured  using  stacked  foil  activation  tcciim.. 
Titanium  (Ti47(p,n)V47)  and  Copper  (Cu^fpnjZn''  i  n  ...  i„ 
were  used  having  threshold  energies  of  3.7  and  4.2  McY. . 
tively.  A  silver  activation  neutron  counter  recorded  nm::. 
produced  by  these  reactions  and  by  accelerated  prom,,-  , 
ing  the  drift  tube  walls. 

Experimental  data  has  been  obtained  for  two  uiifrii:.:  .-. 
celerating  gradients;  one  at  40  MeV/m  over  a  50  cm  ..... 
atmg  distance  and  one  at  90  MeV/m  over  the  same  di-m:  • 
Data  from  the  five  wall  current  probes  for  the  smi.lc  :  r: 
ent  are  shown  in  Fig.  2  for  a)  the  case  when  the  la-n  i  :. 
200  ns  in  advance  of  the  beam,  b)  the  case  where  inc  • 
not  fired  at  all,  and  c)  the  case  where  the  laser  tinum  - 
that  the  plasma  is  produced  by  laser-target  interaction-  '  ■ 
same  time  as  the  beam  is  being  injected  (optimal  l',:i:::.  *■ 

can  be  readily  seen  from  these  results,  good  control  «.v 
beamfront  motion  has  been  achieved  when  the  laser  o.  : 
ing  is  optimal.  Measurements  of  the  accelerated  ion  < ■■.  :  r 
propagated  current,  shown  in  Fig.  3.  clearly  show-  th.v  ' 
can  be  accelerated  up  to  the  designed  output  enerry  v  :.. 
laser  -beam  timing  is  such  as  to  allow  control  of  the  u<  n  •:  ■ 

Experimental  data  obtained  at  the  higher  accci'  -  .m.  r  . 
dient  of  90  MV/m  also  indicate  good  control  over  1 1-<- 
front  motion,  but  peak  ion  energies  observed  are  n>  in.  •  ■'» 
than  those  observed  for  the  lower  gradient  experiment  -  1 : 
undoubtedly  due  to  a  reduction  in  the  electric  field  sti  *■ 
the  virtual  cathode  below  90  MV/m  at  some  point  u;  > '  ■ 

eration  process.  The  following  results  of  numerical  mi . ' 

of  the  experiments  shed  additional  light  on  this  re- 

III.  Numerical  Simulations 

A  particle-in-cell  code  was  used  to  simulate  the  in-  :  ■ 
led  acceleration  experiments.  In  the  simulation-,  n.'  1 
beam  is  assumed  to  be  focussed  by  an  infinitely  st  i .  '  r 

magnetic  field  so  that  the  particles  mo”e  only  alone  tt"1  " 
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^  drift  tube.  The  radius  of  the  beam  is  also  assumed  to  be 
gptch  loss  than  the  radius  of  the  drift  tube  so  that  the  charge 
(jeOSIty.  current  density,  and  axial  electrical  field  are  approx- 
jgUttely  constant  across  the  beam  cross-section.  Ionization  of 
glie  neutral  gas  is  modelled  by  keeping  track  of  the  amount  of 
Kjaization  produced  by  collisions  with  electrons  and  ions  and 
introducing  electrons  and  ions  appropriately. 

In  the  simulations,  the  laser-produced  plasma  is  assumed 
to  completely  neutralize  the  space-charge  on  the  axis  of  the 
jnfi  tube  once  it  reaches  the  electron  beam  from  the  wall, 
^be  lime  required  for  the  laser- produced  plasma  to  reach  the 
fleet ron  beam  from  the  wall  is  assumed  to  be  given  by  the  time 
Inquired  for  a  proton  to  fall  through  a  logarithmic  potential 
deop  r„  from  r  =  /?„.  to  r  =  /?6.  where  l*0  is  the  electron  beam 
toltage,  R.j,  is  the  wall  radius  and  /?(,  is  the  beam  radius. 

Results  are  shown  for  a  900-kY.  20-kA,  1-cm-radius  elec¬ 
tion  beam  which  is  injected  into  a  5-cm-radms,  50-cm-long 
drift  tube  with  a  2-cin-widc,  100-inTorr  hydrogen  gas  cloud  lo¬ 
cated  next  to  the  anode  plane.  The  front  of  the  laser  beam  is 
atfumed  to  travel  down  the  drift  tube  at  a  velocity  which  in¬ 
creases  linearly  from  J,  =  0.04  to  Jj  =  0.2  over  a  distance  of  45 
cm,  corresponding  to  an  accelerating  gradient  of  40  McY/m. 

Figure  4  shows  the  peak  proton  energy  measured  at  45  cm 
versus  the  time  delay  between  the  start  of  the  laser  pulse  and 
the  start  of  the  beam  pulse  rt  -  tb.  In  plotting  this  data,  we 
iMtuned  that  the  laser  requires  10  ns  to  produce  plasma  after 
Striking  the  target  on  the  wall.  For  a  wide  range  of  tl  —  tb. 
the  peak  proton  energy  which  was  measured  actually  exceeds 
the  design  value  of  18.76  McY. 

Figure  5  shows  the  velocity  versus  position  for  an  accel¬ 
erated  proton  for  tl  -  tb  =  -8  ns.  Also  shown  in  the  figure 
it  :Le  vciocuy  versus  position  for  the  front  of  the  laser  beam. 
The  proton  is  accelerated  relatively  smoothly  from  an  initial 
velocity  of  0.04c  to  a  final  velocity  of  0.2c. 

In  all  runs  the  peak  electric  field  £,  fell  by  an  order  of  mag¬ 
nitude.  e  g-  from  3x  108  V/m  to  5x  107  V/m,  as  the  be  0.111  fiulJt 
moved  downstream.  Figure  6  shows  the  magnitude  versus  lo- 
ration  of  the  peak  E,  at  a  number  of  time  steps  approximatelv 
2  ns  apart  for  one  run.  The  decrease  in  the  strength  of  the  peak 
elcctnc  field  observed  in  the  simulations  is  consistent  with  the 
experimental  observation  that  the  peak  ion  energy  actually  fell 
when  a  higher  accelerating  gradient  vas  used. 

IV.  Conclusions 

In  conclusion,  both  experiments  and  numerical  simulations 
bow  indicate  that  the  laser -controlled  collective  accelerator  con¬ 
cept  is  a  promising  one.  Effective  control  over  the  propagation 
cf  a  virtual  cathode  at  the  front  of  an  mtense  relativistic  elec¬ 
tion  beam  has  been  achieved,  and  pre  as  have  been  acceler¬ 
ated  at  a  rate  of  40  MV/m  over  a  distance  of  about  50  cm. 
Furthermore,  numerical  simulations  indicate  that  significantly 
kigher  ion  energies  can  be  achieved  by  either  using  longer  ac¬ 
celerating  distances  land  consequently,  longer  injected  electron 
beam  pulse  durations)  and/or  by  injecting  higher  energy  elec¬ 
tion  beams  to  maintain  higher  electric  field  strengths  at  the 
urtual  cathode. 
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Fig.  1.  Experiment?’  Configuration. 


Fig.  2.  Data  from  5  axially  spaced  current 
collectors  located  on  drift  tube  wall 
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Abstract 


In  this  paper  we  consider  the  propagation  of  an  intense 
electron-ion  beam  across  an  applied  magnetic  field.  In  the  ab¬ 
sence  of  the  applied  field,  the  beam  system  is  in  a  Bennett 
equilibrium  state  that  involves  electrons  with  both  large  axial 
and  thermal  velocities  and  a  cold  stationary  space-charge  neu¬ 
tralizing  ion  species.  Typical  parameters  under  consideration 
are  V„  ~  1  MY.  I  ~  5  kA.  Te  ~  100  keV.  and  beam  radii  ~ 
1  cm.  We  find  that  in  the  intense  beam  regime,  the  propaga¬ 
tion  is  limited  due  to  space-charge  depression  caused  by  the 
deflection  of  the  electron  beam  by  the  transverse  field.  This 
critical  field  is  of  the  order  of  the  peak  self-magnetic  field  of 
the  electron  beam  which  is  substantially  higher  than  the  singie 
particle  cut-off  field. 


Introduction 

The  objective  of  the  present  investigation  is  to  provide 
a  theoretical  understanding  of  the  experiments  performed  at 
the  University  of  Maryland  on  electron-ion  beam  propagation 
across  a  transverse  magnetic  field[l],(2j.  The  configuration  in¬ 
vestigated  is  shown  in  Fig.  1.  Nominal  diode  parameters  are 
1  MV,  20  kA.  and  a  FWHM  of  30  nsec.  The  beam  is  injected 
through  a  hole  into  an  evacuated  drift  tube.  In  order  to  ob¬ 
tain  effective  beam  propagation  to  the  end  of  the  drift  tube, 
a  localized  gas  cloud  near  the  injection  hole  is  employed  and 
collectively  accelerated  ions  from  the  gas  cloud  provide  a  chan¬ 
nel  of  neutralization  that  allows  electrons  injected  late  in  the 
puise  to  propagate.  Under  optimum  conditions  20  kA  can  lie 
propagated  to  the  Faraday  cup  with  no  applied  axial  magnetic 
field. 

In  order  to  investigate  the  effects  of  a  transverse  magnetic 
field  on  this  co-movmg  electron-ion  beam,  a  Helmholtz  coil  is 
placed  about  half-way  down  the  drift  tube.  If  no  gas  is  puffed 
into  the  system,  that  is,  no  ions  are  supplied,  about  1  kA  is 
observed  at  the  Faraday  cup  with  zero  transverse  magnetic 
field.  As  the  transverse  field  is  increased,  the  observed  current 
steadily  decreases  to  zero  at  about  200  gauss.  When  the  gas 
cloud  is  present,  about  5  kA  are  detected  at  the  Faraday  cup 
with  zero  transverse  field.  As  the  transverse  field  is  increased, 
the  entire  5  kA  is  propagated  up  to  about  200  gauss  with  a 
steady  but  slow  decrease  in  collected  current  above  200  gauss. 
The  detected  current  goes  to  zero  around  500  gauss.  The  pur¬ 
pose  of  these  studies  is  to  evaluate  the  properties  of  co-moving 
electron  and  ion  beams  for  use  in  collective  field  accelerators 
and  for  beam  propagation  into  free-spacc. 

Previous  theoretical  studies  have  shown  the  existence  of 
a  self-consistent  downstream  Bennett  equilibrium  for  the  elec¬ 
trons  and  ions  when  no  applied  magnetic  field  is  present(3|.  We 
have  related  these  downstream  properties  to  the  diode  voltage, 
the  transmitted  electron  beam  current,  and  the  ion  properties 


in  the  localized  gas  cloud  region.  For  our  experimental  param¬ 
eters.  V'0  =  1  MV,  /  =  5  kA.  the  downstream  self-pinched  equi¬ 
librium  state  is  composed  of  cold  ions  with  low  axial  speed  ( < 
0.05  c),  and  electrons  with  a  temperature  of  about  60  keV  and 
an  axial  speed  of  about  0.84  c.  The  nearly  charge-neutral  beam 
system  has  been  shown  to  effectively  propagate  up  to  the  diode 
current  of  20  kA  where  the  electron  temperature  is  predicted 
to  be  about  200  keV.  The  5  kA  case  was  chosen  for  examining 
the  effects  of  a  transverse  magnetic  field  on  the  electron-ion 
beam  system,  mainly  because  of  experimental  reproducibility. 

Other  work  in  the  area  of  beam  propagation  across  trans¬ 
verse  magnetic  fields  usually  consider  a  charge  and  current  neu¬ 
tralized  ion  beam  [4]— [8] .  These  systems  have  application  to 
ion  beam  heating  of  dense  plasmas,  where  an  intense  ion  beam 
is  charge  and  current  neutralized  by  electrons  so  that  it  can 
penetrate  the  large  applied  fields  in  a  fusion  plasma.  In  these 
cases  a  two  layer  polarization  model  is  adopted  in  the  rigid 
beam  limit.  The  polarization  electric  field  allows  the  beam  to 
propagate  across  the  magnetic  field  by  means  of  an  £  x  B  drift. 


Model,  Basic  Equations,  and  Results 


Consider  a  rigid  beam  model  of  the  Bennett  equilibrium 
system  discussed  in  the  previous  section.  Electron  and  ion 
beams  of  equal  radii  are  incident  along  the  z-axis  onto  a  trans¬ 
verse  magnetic  field  region.  To  simplify  the  analysis  the  profile 
of  the  Helmholtz  coil  is  replaced  by  a  constant  magnetic  field 
Bt  in  the  y-direction  over  an  axial  distance  Lt  ~  20  cm.  In  the 
parameter  regime  of  interest,  the  electron  gyroradius  is  com¬ 
parable  to  Lj  which  is  much  larger  that  the  effective  Bennett 
beam  radius  a  and  the  ion  beam  is  assumed  to  be  stationary 
and  centered  on  the  z-axis.  The  electron  beam  is  confined  by 
the  self-magnetic  field.  The  rigid  electron  beam  is  deflected  by 
the  transverse  field  Bt  but  feels  a  restoring  force  due  to  the  ion 
beam.  In  this  model,  we  examine  only  the  effects  of  transverse 
forces  on  the  rigid  electron  rod.  Longitudinal  forces  and/or 
effects  that  arise  as  the  beam  front  end  enters  the  transverse 
magnetic  field  region  are  ignored. 

When  there  is  no  transverse  field.  Bt  =  0.  the  electron  and 
ion  components  of  the  beam  system  arc  centered  on  the  z-axis, 
are-charge  neutral,  n,„  =  nto  =  n„,  and  have  identical  Bennett 
profiles  of  the  same  effective  radius,  a,  =  ac  =  a.  That  is.  the 
density  profiles  are 


no(r) = 


(1  +  (r/o)2] 


-,  a  =  t.  e 


(1) 


and  the  electron  beam  self-magnetic  field  due  to  the  mean  axial 
electron  velocity  V„  is 


Os\>lr) 


2  (r/n)2' 


(2) 


The  ions  are  stationary  and  assumed  infinitely  massive.  The 
equations  of  motion  for  the  center  of  the  electron  beam  are 
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dp  r 
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^  =  -eBTVt.  (4) 

at 

where  pz  =  m-v\'r,  p,  =  my\',,  and  7  =  [1  — (Vr/c)2—  (l,/c)2]  1/2 
The  averaged  ion  electric  force  acting  on  the  rigid  electron 
beam  when  the  electron  rod  is  displaced  a  distance  x  is 


F„{x'<  -  [  Er,(r)nt(r)dS' 

-e2n0a  f  x/2  1 

S<0  IT  +  (i/2)2  "  (172) 
sinh-'i  -i/2)  -  sinh-*(3x/2  4-  i3/2) 

x2[l  -r  (x/2)2j3/- 


where  x  =  xja. 

In  general.  Eqns.  (3)  and  (4)  have  to  be  solved  numerically. 
However,  if  the  displacement  distance  x  is  small  compared  to 
the  beam  radius.  i<l.  the  ion  electric  restoring  force  can  be 
approximated  as 


Fe,(x) 


-e‘n<, 

- j 

6c„ 


(6) 


Furthermore,  assuming  the  relativistic  factor  y  remains  con¬ 
stant,  7  =  y0  =  (1  -  ( V„/c)2]"1/2.  we  find  the  solutions  to 
Eqns.  (3)  and  (4)  to  be 


and 


witn 


iff)  =  "~1(\  -  cosud) 

[  .  vJ 

V'*(<)  =  X’,0  |1  -  ^  ~J  (1-cosu :<)  , 

2  _  2  ,  “’p  _eBT  2  _  ***<> 

vJ  =  W_  T"  ~T~  j  —  i  * 

6  my0  eom-r0 


(<) 
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where  are  the  electron  cyclotron  and  plasma  frequencies. 

These  solutions  are  similar  to  the  remits  obtained  by  Peter 
and  Rostoker(8],  who  used  the  two-layer  polarization  model. 
Since  the  oscillatory  behavior  of  the  electron  rod  is  due  to  the 
deflecting  applied  transverse  magnetic  field  and  the  restoring 
ion  beam  electric  field  force,  we  see  that  the  frequency  of  these 
oscillations  is  a  function  of  both  wc  and  up.  Averaging  out  the 
oscillations  in  Eq.  (8).  we  obtain  for  the  time-averaged  mean 
axial  velocity  of  the  electron  beam 

v;  -  r„  (\  - 

where  e  =  1  +  w2/6w2.  This  result  is  identical  to  that  in  refer¬ 
ence  [S]  with  the  modified  definition  of  c. 

Returning  to  Eqns.  (3)  and  (4).  we  have  solved  these  equa¬ 
tions  numerically  for  a  given  beam  current  as  a  function  of  the 
applied  transverse  magnetic  field  amplitude.  The  results  for  a 
5  kA  beam  are  displayed  in  Fig.  2.  From  reference  [3).  the  self- 
consistent  ;  an  electron  axial  velocity  is  VI0  =  .84c  for  a  beam 
radius  of  1.0  cm.  The  peak  self-magnetic  field  for  this  beam 
system  is  about  500  G  [from  Eqn.  (2)).  The  results  displayed 
are  typical  for  beams  in  the  "intense  '  regime.  Specifically,  we 
have  plotted  in  Fig.  2a  the  maximum  beam  displacement  ver¬ 
sus  the  transverse  field  amplitude.  Dj,  shown  as  the  solid  line 
In  Fie  2b.  the  minimum  axial  velocity  of  the  electron  beam, 
winch  occurs  at  maximum  displacement,  is  plotted  versus  Br 
l  solid  line  I.  As  shown,  the  maximum  displacement  continu¬ 
ally  increases  with  Bf  until  a  value  is  reached  such  that  the 
axial  velocity  goes  to  zero.  This  vaiue  of  the  magnetic  field 


-7  =r„(l-l/e). 


(10) 


at  which  — *  0.  is  defined  as  Btc  and  the  displacement 

amplitude  for  this  field  xc.  For  values  of  Bt  above  Btc.  the 
beam  will  blow  up  in  its  own  space-charge.  We  see  that  this 
critical  value  is  very  close  to  the  peak  seif-magnetic  field  of  the 
Bennett  beam. 

The  values  of  Btc  and  x.  can  be  related  from  the  constants 
of  motion.  From  Eqns.  (3)  and  (4).  these  constants  are  the 
total  energy  and  the  axial  canonical  momentum,  which  can  be 
written  as 


and 


(7-1  ) me1  -  eo,  =  i  7„  -  1  )mc~ 


(11) 


P:  +  eBjx  =  p!0  =  my„X',0,  (12) 

where  the  effective  potential  is  given  as  <j>e  =  (1/e)  /*  Fe,(x')dx‘. 
From  Eqn.  (12).  when  pt  =  0,  then  x  — *  ic  and  Bf  — >  Btc. 
that  is 

Btc*c  =  —  =  m~r0e\’.,.  (13) 

c 

In  the  intense  beam  regime,  when  J*  goes  to  zero  at  Btc 
we  also  have  Jr  =  0.  That  is.  at  this  critical  magnetic  field 
~t(imax  =  xc)  =  1,  all  the  electron  beam  energy  goes  into 
potential  energy.  From  Eqn.  (11),  with  7  =  1,  we  obtain 


COc(Xc) 


4  "  f -j  1  to 


=  (*V0  -  1  )mc‘ 


i  14, 


Thus,  fror.,  Eqn.  (14)  we  can  determine  xc  for  a  given  injected 
beam,  and  determine  the  critical  magnetic  field  from  Eqn.  ( 13). 
That  is. 


Btc 


W7,'» 

exc 


my, ,  r. 

e\/3  a 


r 

|  exp 


the  critical  magnetic  field  for  beam  propagation  in  the  intense 
beam  regime. 

If  the  beam  current  is  sufficiently  low.  /  <  4/tA.  we  find 
that  when  J.  —  0.  Jr  is  not  zero.  In  the  single  particle  limit, 
this  is  obvious  since  as  J.  —  0.  Jr  —•  J„,  since  7  =  y0  along  the 
entire  trajectory.  In  this  case  xc  =  p,„/eBTc  =  l's»/—’e  —  B,  = 
Lt.  the  electron  Larmor  radius.  For  beam  currents  between 
the  single  particle  limit  and  the  intense  beam  regime  discussed 
above,  Eqn.  (11)  tells  us  that  part  of  the  initial  beam  kinetic 
energy  goes  into  potential  energy  as  the  beam  is  deflected.  A 
summary  of  these  results  is  shown  in  Fig.  3  where  we  have 
plotted  the  critical  magnetic  field  Btc  versus  beam  current 
I.  This  graph  is  generated  for  the  case  of  a  1  MV  diodc(3j 
and  the  self  consistent  downstream  Bennett  equilibrium.  The 
mean  axial  velocity  X*ta  depends  on  the  transmitted  current 
as  shown  in  refercnce(3|  and  this  is  included  in  the  plot  in 
Fig.  3.  We  see  that  the  single  particle  cutoff  field  is  about 
250£7[l„  =  l.l/V.  Lt  =  IScm  =  R,}.  As  the  current  is  increased 
the  critical  field  remains  at  the  single  particle  limit  until  t lie- 
beam  current  readies  a  value  such  that  its  peak  self-magnetic 
field  is  near  the  single  particle  cut-off  field.  Further  increases 
in  the  beam  current  result  111  entering  the  intense  beam  regime 
where  the  critical  field  is  a  result  of  potential  depression  (7  —  1 

Jr>  « 

t  —  u  /• 

In  order  to  explicitly  compare  our  results  with  those  of  the 
experiment,  we  must  normalize  our  uniform  field  model  with 
a  Helmholtz  coii  profile.  We  have  displayed  in  Fig.  2a  and 
2  b  the  numerical  results  of  solving  Eqns.  (31  and  (4)  for  a 
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5  kA  electron  beam  transversing  a  Helmholtz  coil  configura¬ 
tion.  The  coii  field  amplitude  is  the  peak  on-axis  value.  From 
conservation  of  axial  canonical  angular  momentum  and  energy 
conservation,  the  relationship  between  the  critical  uniform  held 
Bt~  and  that  ot  the  Helmholtz  coii  field  profile  is 


BjC 


Discussion 


(1C) 


For  a  given  injected  electron  beam  that  is  charge  neutral¬ 
ized  but  not  current  neutralized  by  a  stationary  ion  beam  and 
has  a  Dennett  equilibrium  profile,  we  have  determined  the  crit¬ 
ical  transverse  magnetic  field  above  which  the  beam  will  not 
propagate.  In  the  intense  beam  regime,  this  limit  is  a  result 
of  space-charge  depression  caused  by  the  deflection  of  the  elec¬ 
tron  beam  from  the  center  of  the  ion  channel.  The  value  of  the 
critical  field  is  approximately  given  by  the  peak  self-magnetic 
field  of  the  electron  beam. 

A  concern  about  the  ngid  beam  model  is  the  individual 
particle  confinement  in  the  Bennett  equilibrium.  The  appli¬ 
cation  of  a  transverse  field  results  in  a  non-symmetnc  force 
about  the  beam  axis.  We  have  examined  the  individual  par¬ 
ticle  trajectories  for  typical  Bennett  particles  and  have  found 
the  non-symmetry  has  substantial  individual  particle  effects  as 
Bt  approaches  Bto\ai.  We  believe  that  this  will  lead  to  frac¬ 
tional  ioss  of  beam  current  as  the  beam  transverses  the  field 
region.  This  fractional  loss  is  observed  in  the  experiment  as 
Btc  is  approacned. 
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Fig.  2.  Typical  results  in  the  Intense  Beam  Regime.  ia) 
Maximum  electron  beam  displacement  and  (In  minimum  beam 
axial  velocity  versus  the  applied  transverse  macnetic  field  am¬ 
plitude.  The  "solid"  line  represents  the  constant  field  model 
and  the  "dashed"  line  a  Helmholtz  coii  configuration  l  peax  on- 
axis  amplitude).  Fixed  beam  parameters  arc;  I  =  6  kA,  Jto  = 
0.S4.  a  -  1  cm.  and  I?,*(a)  —  500  G. 
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Fig.  1.  Experimental  configuration  used  for  beam  propaga¬ 
tion  studies  across  a  magnetic  field. 
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Fig.  3.  Critical  transverse  magnetic  field,  defined  by  d,  = 
0.  versus  electron  beam  current.  The  "dashed"  line  represents 
the  peak  self-magnetic  field  oi  the  electron  beam.  The  beam 
radii  arc  1  cm  and  the  diode  voltage  is  1MV.  See  reference  (3|. 
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ABSTRACT 

We  report  the  first  measurement  of  the  rms  emittance  of  a  pseudospark-produced 
electron  beam.  A  six  gap  pseudospark  chamber  filled  with  argon  gas  was  operated 
at  ~25  kV.  producing  a  ~10  Hz  repetitive  pulse  train  of  electron  beams.  Typically, 
a  beam  of  average  energy  20  keY.  peak  current  ~50  A  and  pulse  duration  ~10 
ns  FWHM  was  extracted  from  the  chamber.  The  rms  emittance  was  evaluated  by 
using  the  slit-hole  method.  A  typical  value  of  rms  emittance  was  found  to  be  e  ~  55 
mm-tnrad.  yielding  a  normalized  emittance  of  en  ~  15  mm-mrad.  The  normalized 
brightness  of  the  beam  was  then  estimated  as  Bn  =  //e2  ~  2  x  10n  A/(  nrrad* ). 


I.  INTRODUCTION 

Interest  in  high  quality,  high  current  electron  beams  has  been  stimulated  by 
their  application  with  regard  to  such  field  as  advanced  accelerators  and  free  elec¬ 
tron  lasers.12  Pseudospark  discharge  phenomenon  with  interesting  charged  particle 
emission  characteristics  was  reported  by  Christiansen  and  Schultheiss3  in  1979.  At 
the  University  of  Maryland,  an  experiment  was  performed  to  measure  for  the  first 
time  the  emittance  of  an  electron  beam  produced  by  a  pseudospark  discharge. 

In  this  paper,  we  shall  describe  the  operation  of  a  six  gap  pseudospark  chamber 
which  is  characteristically  similar  to  the  devices  reported  by  other  laboratories."-' 
Repetitively  pulsed  ( ~10  Hz)  electron  beams  were  extracted  through  an  anode  hole 
into  a  drift  region.  The  beam  currents  were  measured  by  a  Rogowski  coil  at  the 
anode  and  a  Faraday  cup  placed  downstream.  A  simple  emittance  meter8  consisting 
of  a  series  of  thin  slits  and  detector  film  was  placed  downstream  of  the  anode. 
Angular  distributions  of  sheet  beamlets  formed  after  passing  through  the  slits  were 
recorded  on  radiachromic  film.1'  The  optical  density  distribution  of  the  film  after 
exposure  to  the  beam  provide  us  with  the  distribution  of  the  transverse  components 
of  electron  velocities  from  which  the  rms  emittance  is  evaluated.  Subsequently,  we 
describe  the  emittance  analysis  and  estimation  of  brightness  of  the  beam.  Finally. 


the  experimental  results  are  discussed. 


II.  EXPERIMENTAL  ARRANGEMENT 

The  experimental  setup  is  shown  in  Fig.  1.  The  discharge  chamber  consists 
of  a  hollow  cathode,  five  modules  of  intermediate  electrodes  and  insulators,  and 
an  anode  with  a  center  hole  for  electron  beam  extraction.  Intermediate  electrodes 
made  of  3.2  mm  thick  brass  have  outer  diameter  of  6.35  cm  and  center  hole  diameter 
of  3.2  mm.  The  3.2  mm  thick  plexigias  insulator  washers  have  outer  and  inner 
diameters  of  7  cm  and  2.54  cm  respectively.  In  addition  to  the  chamber  capacitance 
of  11  pF  between  anode  and  cathode,  a  low  inductive  external  capacitor  of  420 
pF  was  added  for  the  present  work.  A  50  cm  long  drift  chamber  is  attached  to 
the  anode  side  to  accommodate  diagnostics  such  as  the  emittance  meter  or  the 
Faraday  cup.  A  capacitance-manometer  type  vacuum  gauge  was  used  to  measure 
the  gas  fill  pressure  which  was  almost  statically  balanced  by  a  needle  valve  and  a 
throttle  vacuum  valve  as  shown  in  Fig.  1.  Two  Rogowski  coils  were  molded  into 
axisymmetric  grooves  milled  into  both  side  of  the  anode  flange  so  that  the  azimuthal 
component  of  the  magnetic  field,  which  arises  from  the  axisymmetric  current,  is 
predominantly  supported  and  other  components  (noises)  are  suppressed.  The  rise 
time  of  the  Rogowski  coil  system  may  be  approximately  given  by  F/Zo,  where  L  is 
the  inductance  of  the  coil  and  Z0  is  the  characteristic  impedance  of  the  transmission 
line  used,  and  is  found  to  be  less  than  0.5  ns.  A  Faraday  cup  consisting  of  a  3  cm 
diam  graphite  beam  collector  and  a  10  mfi  current  viewing  resistor  placed  on  axis 
of  the  downs; i earn  drift  chamber  was  employed  to  measure  the  l" earn  current  at 
various  axial  positions.  The  chamber  voltage  was  measured  by  a  homemade  high 
impedance  resistive  voltage  probe  of  division  ratio  1:20.000  into  a  50  ft  load.  The 
response  of  the  probe  was  greatly  improved  by  a  proper  compensation  for  the  stray 
capacitance  m  the  Coxbon  resistors  used.  The  resulting  rise  time  is  less  than  0.5 
ns.  and  the  RC  droop  time  constant  due  to  the  blocking  capacitor  is  ~  120  ns. 
For  emittance  measurement,  a  slit-hole  type  emittance  meter  was  employed  and 
was  placed  in  the  drift  chamber  9  cm  downstream  of  the  anode.  The  emittance 
meter  consists  of  a  series  of  parallel  thin  slits  of  200  /jib  width  and  2  mm  spacing 
constructed  from  0.6  mm  thick  stainless  steel  plate;  2  mil  thick  radiachromic  film, 
used  as  a  beam  detector  was  placed  12  mm  downstream  of  the  slit  plane. 


Downstream'  Upstreom 


Fig.  1.  Experimental  setup. 
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The  system  was  operated  typically  at  30  kY.  with  an  argon  gas  pressure  of  30 
mTorr. 

For  the  eiruttance  measurement,  exposures  of  approximately  rive  consecutive 
beam  pulses  to  the  raaiachromic  film  through  the  slits  in  the  emittance  meter 
were  enough  to  produce  an  appropriate  density  profile,  with  the  peak  value  of  the 
optical  density  not  exceeding  0.5.  This  ensure  that  the  measured  optical  density 
distribution  is  linear  to  the  beam  intensity.9  The  density  profile  was  obtained  by 
scanning  the  film  using  an  optical  microdensitometer,  and  is  shown  in  Fig.  2. 


III.  ROOT-MEAX-SQUARE  EMITTANCE  ANALYSIS 

The  root-mean-square  ( rms  I  emittance10  /or  effective  emittance11'.  which  has 
been  widely  used  as  a  measure  of  beam  quality,  was  defined  as 

c  =  4  |<  .r"  ><  x'1  >  —  <  xx'  >‘j  .  ill 

where  x'  is  the  gradient  of  the  particle  trajectory  given  by  r'  =  dx Id:  =  }>r/p:.  and 
the  angular  brackets  denote  average  values  over  the  two-dimensional  trace  space  as 

<  o  >=  J  op(x.  x')dxdx'.  :2) 

wiiere  p  is  th?  projected  density  in  two-dimensional  trace  space,  and  is  assumed  to 
be  normalized,  i.e..  J  pdxdx1  =  1. 

It  is  very  reasonable  to  assume  that  as  in  Ref.  S.  the  beam  produced  in  this 
experiment  is  axisymmetric  and  of  Maxwellian  transverse  velocity  distribution.  This 
allows  us  to  use  the  simple  slit-hole  type  emittance  meter  mentioned  above,  whose 
results  can  be  easily  analyzed.8  We  find  empirical  functions  o(r).  J(r),  and  crt r )  as 
functions  of  radial  position  r  from  the  density  profile  shown  in  Fig.  2.  where  a(  r )  is 
the  mean  diverging  angle.  3(r)  represent  the  peak  values,  and  rr(r)  is  the  rms  width 
of  the  individual  distributions.  Numerical  integrations  were  then  performed  using 
the  empirical  functions  a.  3.  and  a  to  find  p  and  the  moments  <  r*  >.  <  x':  >. 
and  <  xx'  >  i  see  Ref.  5  for  details).  Several  isodensity  contours  of  the  resultant 
p tx.x'i  were  constructed  in  x  -  x'  space  (known  as  an  ennttance  plot)  as  shown 


Fig.  2.  Optical  density  distribution  on  the  radiachromic 
film  after  exposure  to  the  beamlets. 


in  Fig.  3.  The  measured  rms  emittance.  given  by  Eq.  ( 1 )  with  the  obtained  moments, 
was  found  to  be 


f  ~  55mm  —  mrad.  (  3) 

It  is  interesting  to  note  that  the  emittance  plot  in  Fig.  3  shows  a  slightly  S-shaped 
contour;  this  is  indicative  of  the  existence  of  nonlinear  focusing!  defocusing)  in  the 
beam  indicating  the  possible  existence  of  space  charge  field  of  substantial  strength. 
The  normalized  rms  emittance  en  =  3a e.  which  is  invariant11,12  when  there  is  ac¬ 
celeration  in  the  axial  direction,  is  conveniently  used  for  comparing  beam  qualities 
at  different  energies.  With  an  average  beam  energy  of  20  keV  inferred  from  the 
voltage  probe  signal,  the  normalized  emittance  is  estimated  to  be  en  ~  15  mm- 
mrad.  Another  useful  invariant  of  the  beam  associated  with  the  emittance  is  the 
brightness."  We  use  here  a  simpler  mathematical  expression  of  normalized  rms 
brightness  Br,  =  //c2.  so  as  to  compare  with  other  results  by  other  laboratories 
using  different  definitions.13  With  a  measured  current  of  ~50  A.  we  obtained  the 
brightness 


=  //e2  ~  2  x  10"  A/m2rad2.  (4) 

which  is  an  order  of  magnitude  higher  than  that  of  other  high  brightness  sources.13 
It  should  be  noted  here  that  the  results  in  this  work  are  all  time  integrated  over  5 
shots. 


IV.  CONCLUSIONS 

We  have  measured,  for  the  first  time,  the  emittance  of  a  pseudospark-produced 
electron  beam.  A  six  gap  pseudospark  chamber  with  argon  gas  was  operated  at  ~ 
25  kV  and  produced  a  ~  10  Hz  repetitive  pulse  train  of  electron  beams.  Typically, 
a  beam  of  average  energy  20  keV.  peak  current  ~50  A  and  pulse  duration  ~10 
ns  FWHM  was  extracted  from  the  chamber.  The  rms  emittance  was  evaluated  by 
using  the  slit-hole  method.  The  typical  value  of  rms  emittance  was  found  to  be  e  ~ 
55  mm-mrad.  yielding  a  normalized  emittance  of  en  ~  15  mm-mrad.  These  values 


x(mm) 


Fig.  3.  The  emittance  plot  in  x-x'  trace  space. 
The  isodensitv  contours  correspond  to  0.1.  20. 
40.  60,  and  80%  of  the  peak  value. 


are  time  integrated  over  5  shots:  thus  they  represent  upper  limits.  The  normalized 
rms  brightness  of  the  beam  was  then  estimated  as  Bn  =  I  j  ei  ~  2  x  1 0 1 1 A  / 1  nr  rad J ). 
which  is  an  order  of  magnitude  higher  than  existing  high  brightness  electron  beam 
sources. 
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A  comprehensive  analysis  of  pulse  formation  in  inductive-energy  pulsed  power  circuits  is 
presented.  The  output  waveforms,  with  relevant  circuit  parameters,  are  derived  for  four  basic 
systems,  including  two  new  circuit  configurations  useful  for  practical  applications. 


I.  INTRODUCTION 

Inductive-energy  pulsed  power  systems  (PPS’s)  have 
recently  attracted  a  great  deal  of  attention,  mainly  because 
of  their  compactness  in  size  as  compared  to  the  more  often 
used  capacitive  PPS’s.  Due  to  relatively  slow  progress  in 
the  development  of  fast  opening  switches,  most  of  the 
inductive-energy  PPS’s  now  in  use  have  been  limited  to 
slow  pulse  systems'  in  which  a  simple  lumped  inductor  is 
employed  a»  the  energy-storage  element.  It  has  been 
shown2  that,  in  lieu  of  the  lumped  inductor,  a  transmission 
line  can  be  used  to  store  the  inductive  energy  in  order  to 
produce  a  square  output  pulse.  Such  square-pulse  genera¬ 
tion  has  b..en  experimentally  demonstrated  by  using  a 
light-activa  ed,  fast  semiconductor  opening  switch.3  As  the 
developme  t  of  opening  switches  progresses,  one  may  ex¬ 
pect  to  see-  -as  in  the  case  of  capacitive  PPS’s — a  demand 
for  a  vari  ty  of  circuit  configurations  in  order  to  meet  the 
requiremen  j  of  the  many  and  diverse  applications  that  will 
naturally  arise. 

It  is  realized  that  there  exists  a  “dual”  relationship4 
between  tne  capacitive-  and  the  inductive-energy  PPS’s. 
One  can  construct  an  inductive-energy  PPS  circuit  from  a 
capacitive  mergy  PPS  circuit  by  replacing  capacitance 
with  inductance,  inductance  with  capacitance,  resistance 
with  cond  dance,  closing  switch  with  opening  switch,  and 
series  corn  action  with  parallel  connection.  The  behavior  of 
voltage  in  the  dual  circuit  is  the  same  as  that  of  current  in 
the  origini ;  circuit,  and  vice  versa. 

In  this  paper  we  present  a  comprehensive  analysis  of 
pulse  formation  in  four  basic  inductive-energy  PPS  cir¬ 
cuits,  including  two  new  circuits.  Most  .ther  possible  sys¬ 
tems  may  be  considered  as  a  variation  and/or  combination 
of  the  above  systems.  The  circuits  described  in  the  follow¬ 
ing  sections  are  ( 1 )  the  lumped  inductor  system,  (2)  the 
current  charged  transmission  line,  ( 3 )  the  dual  of  the  LC 
generator,3  and  (4)  the  dual  of  the  Blumlein  line.5  Also 
discussed  are  the  features  and  results  of  each  system  in 
conjunction  with  possible  applications. 

II.  LUMPED  INDUCTOR  SYSTEM 

This  may  be  the  simplest  form  of  inductive-energy 
PPS’s.  It  consists  of  a  current  charged  inductor  and  an 
opening  switch  with  a  capacitor  connected  in  parallel,  as 
shown  in  Fig.  1(a).  The  stored  energy  in  the  inductor, 
\LIq  is  released  to  the  load  resistor  by  the  opening  switch. 
This  circuit  is  the  well  known  parallel  RLC  circuit,  for 
which  the  governing  equation  is  given  by 


dlv  1  dv  v 
dr+RCdt~LC 


(1) 


The  switch  is  opened  at  f  =  0,  interrupting  the  initial 
charging  current  J0.  Three  different  types  of  solutions  for 
the  ouptut  waveform  exist,  depending  on  the  circuit  pa¬ 
rameters:  (a)  fora>w0  (overdamped), 

v(t)  =  -  2  °  2  |-/;e~atsinh (or  —  <a5)l/:r,  (2) 

c  t  a  —  ct)0 ) 

where  a  =  1/(2 CR)  and  o>0  =  (LC)  "  l/2;  (b)  for  a  =  co0 

(critical  damping), 

v{t)  =  (I0t/C)e-a‘,  (3) 

(c)  for  a<co0  (underdamped), 

v(t)  =  rl,S  ~LzT i /z  sin^o  ~  a‘)ult.  (4) 

c  (a)0  —  a  ) 

These  solutions,  with  typical  parameters,  are  plotted  in 
Fig.  1(b). 

Due  to  the  presence  of  the  parallel  connected  capaci¬ 
tance,  the  ouptut  waveforms  have  a  finite  rise  time,  which, 
in  turn,  results  in  a  reduction  of  peak  amplitude.  The  rise 
time  Tp  (at  which  the  amplitude  peaks)  is  found  by  solving 
dv/dt(,Tp)  sO,  yielding  the  corresponding  peak  amplitude 
Vp  —  v(Tp).  The  expressions  for  Tp  and  Vp  for  the  above 
three  cases  are  listed  in  Table  I. 

It  should  be  mentioned  that  this  simple  system  is  the 
most  popular  electrical  PPS  found  in  practical  applica¬ 
tions.  The  ignition-coil  system  used  in  most  automobile 
engines  is  a  slight  variation  of  this  type.  The  role  of  the 
parallel  capacitance  in  this  case  is  to  optimize  the  rise  time 
in  such  a  way  that  the  finite  rate  of  output  voltage  nse 
matches  the  voltage-holding  characteristics  of  the  mechan¬ 
ical  opening  switch  (known  as  “breaker  points")  as  it 
opens.' 


III.  CURRENT  CHARGED  TRANSMISSION  LINE 

As  shown  in  Ref.  2,  a  transmission  line  can  be  used  in 
lieu  of  the  lumped  inductor  and  capacitor  as  a  pulse- 
forming  network.  The  transmission  line  of  characteristic 
impedance  Z0  is  initially  current  charged  to  /0  as  shown  in 
Fig.  2(a).  In  this  case  the  stored  magnetic  field  energy. 
ISliH2  dv,  is  again  found  to  be  IL/q  (just  as  for  the  lumped 
inductor  system),  where  L  is  the  total  inductance  of  the 
transmission  line. 
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FIG.  1.  Lamped  inductor  system  (a)  schematic  diagram,  and  (b)  typical 
output-vohage  waveforms  showing  overdamped,  critically  damped,  and 
underdamped  cases. 


The  current  /0  T«y  be  considered  the  superposition  of 
a  positively  traveling  wave  (PTW)  V* 

=  ‘jZoIo)  and  a  negatively  traveling  wave  (NTW) 
(/"  —'ilo,  V~  =  —  jZo/q),  both  of  which  are  constantly 
reflected  back  and  forth  by  the  shorted  end  (T  =  —  1 )  and 
the  closed  switch  (T  =  —  1)  until  the  switch  opens  [see 
Figs.  2(b)  and  2(c)].  F  is  the  reflection  coefficient0  (the 
ratio  of  reflected  voltage  wave  to  incident  voltage  wave), 
and  is  given  by  T  =  (R  —  Zq )/(R  +  Zq).  After  the  switch 
opens  (T  =  0)  at  t  =  0,  the  charging  current  is  then  inter¬ 
rupted,  allowing  the  PTW  (J+,K+  )  to  travel  towards  the 
resistive  load,  while  the  shorted  end  converts  the  rest  of  the 
NTW  (I-.V-)  to  a  PTW  [see  Fig.  2(d)].  As  the  PTW 
arrives  at  the  matched  load  (T  =  0),  the  wave  is  com¬ 
pletely  dissipated,  resulting  in  the  generation  of  a  square 
output  pulse,  this  resultant  output  pulse  at  the  matched 
load  has  parameters  /<*,,  =  j /<>  =  j/qZq,  and  r  =  2f/v, 


FIG.  2.  Current  charged  transmission  line:  (a)  schematic  diagram:  spa¬ 
tial  distribution  of  (b)  current  and  (c)  voltage,  shown  as  the  superposi¬ 
tion  of  a  positively  '.'...rL.:  j  wave  and  a  negatively  traveling  wave  before 
the  switch  is  opened  ( r  <  0 ) .  and  ( d )  spatial  distribution  of  voltage  dunng 
the  period  0  <  t  <  f/v. 


where  Sis  the  length  of  the  transmission  line  and  t;  is  the 
wave  velocity  in  the  transmission  line  [see  Fig.  3(a)]. 

If  the  load  is  mismatched  (i.e..  R=£Z0)  the  PTW  ar¬ 
riving  at  the  load  experiences  a  partial  reflection  due  to  the 
mismatch  \  T  =^0);  this  reflected  wave  is  then  again  re¬ 
flected  at  the  shorted  end  and  sent  back  to  the  load.  This 
cycle  repeats  with  gradually  reduced  amplitude,  resulting 
in  a  senes  of  postpulses  after  the  main  pulse  as  shown  in 
Figs.  3(b)  and  3(c).  Note  that  the  two  distinctive  types  of 
resulting  waveforms,  corresponding  to  R  >  Z0  and  R  <Zo, 
respectively,  have  a  resemblance  to  those  of  the  over- 
damped  (a>u>c)  and  underdamped  (a<u)0)  cases  in  the 


TABLE  I.  Riae  time  T,  and  peak  amplitude  Vp  of  the  output  pulse  in  the  lumped  inductor  system  for  three  different  damping  conditions. 
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FIG.  4.  Dual  of  the  LC  generator:  (a)  schematic  diagram,  and  (b) 
typical  output-voltage  waveforms  normalized  to  R/0  for  which  ( LC)  *  1  : 
=  I06,  and  L/R  =  10 ' 4,  10  ’  \  and  10  -  6 


FIG  3  Typical  output-volume  waveforms  of  the  current  charged  trans¬ 
mission  line  with  (a)  matched  load  ( R  =  Zo),  (b)  R  >  Zo  (R  =  2Zo), 
and  (c)  R  <  Zo  (R  =  Zo/2). 

lumped  inductor  system.  The  amplitude  of  the  nth  post¬ 
pulse  for  both  cases  is  given  by 

^u.^Zo/oHZo  — /n7(Z0 +  /?)"* ’].  (5) 

where  n  =  0  corresponds  to  the  main  pulse. 

This  current  charged  transmission  line  system  may  be 
considered  the  dual  of  the  simple  voltage-charged  trans¬ 
mission  line  system;  both  systems  produce  a  square  output 
pulse.  In  the  case  of  lumped  energy-storage  systems,  the 
peak  amplitude  of  the  output  is  restricted  by  the  finite  rise 
time,  which  is  inevitable  due  to  stray  capacitance,  espe¬ 
cially  for  very  short  pulse  generation.  For  the  case  of  the 
current-charged  transmission  line  system,  a  very  high- 
power  pulse  may  be  easily  generated  by  simply  shortening 
the  pulse  duration  for  a  given  stored  energy,  just  as  with  its 
counterpart  of  the  voltage-charged  transmission  line 

IV.  DUAL  OF  THE  LC  GENERATOR 

We  introduce  here  a  new  circuit  as  shown  in  Fig.  4(a). 
Two  inductors,  which  are  current  charged  in  senes,  and  an 
opening  switch  are  utilized  From  the  circuit  point  of  view, 
this  is  a  dual  of  the  LC  generator,  which  is  a  well-known 
capacitive  energy  PPS.  The  circuit  may  be  described  by  the 
following  third-order  differential  equation: 

d3v  2 R  d2v  1  dv  Rv 

d?^Td?"IcJt^Uc=0  (6) 


With  appropriate  initial  conditions  u(0)  =  0,  dv/ 
dt(0)  =  0,  and  d2v/dt2( 0)  =  IqR/LC,  the  solution  (see 
Appendix)  is  found  to  be 

VU)  LC[(a-p)2  +  co2} 

X(e~al  —  e~e'{cos  cot  —  l  (a  —  P)/co] sin  cot}), 

(7) 

where  a  =  2f?(3L)_l  —  (r  +  d)Wi  —  (r  —  d)wi, 
0  =  2  -'[(/■  -t.  d)wl  +  (r  -  d),/3]  +  2 R(3L)-\ 
co  =  31/22  -  '[(r  4-  d)wy  -  (r  -  d)v\d  =  (qy 
i-  r)W2,q  =  (3  LC)-'  -  (2R)2(3L)  - 2,  and  r 
=  —  R(6L2C)  ~  1  -  ( 2R )  3(  31)  “3.  The  behavior  of  this 
solution  depends  largely  on  the  circuit  parameters  R,  L, 
and  C;  three  typical  cases  are  plotted  in  Fig.  4(b). 

The  main  advantage  of  this  system  over  the  simple 
lumped  inductor  system  is  that,  as  with  its  dual  counter¬ 
part,  the  LC  generator,  the  output  current  amplitude  can 
be  as  high  as  twice  the  initial  charging  current. 

V.  DUAL  OF  THE  BLUMLEIN  LINE 

Here  we  introduce  another  new  circuit  that  is  the  dual 
of  the  Blumlein  line,  a  popular  capacitive  energy  PPS.  This 
circuit  may  also  be  considered  as  a  distributed  circuit  ver¬ 
sion  of  the  above  mentioned  dual  LC  generator.  As  shown 
in  Fig.  5(a),  two  identical  transmission  lines  of  character¬ 
istic  impedance  Zq  are  initially  current  charged  in  senes  to 
1q-  This  I0  may  be  considered  a  superposition  of  a  PTW 
(I*  o.  T*=iZo/0)  and  a  NTW  (/-«$/„, 
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V~  =  —  5Z0/0).  which  are  continuously  reflected  by  both 
ends  (T  —  —  1),  maintaining  a  steady  state  as  shown  in 
Figs.  5(b)  and  5(c).  Note  that  both  the  PTW  and  NTW 
experience  a  discontinuity  at  the  midpoint  where  the  load 
resistance  R  is  connected  in  parallel;  it  can  be  shown,  how¬ 
ever,  that  the  voltage  (current)  across  the  load  resistance 
produced  by  both  the  PTW  and  NTW  are  of  equal  ampli¬ 
tude  but  of  opposite  sign.  Thus  no  energy  is  dissipated  in 
the  load,  just  as  if  there  were  no  discontinuity. 

At  t  =  0,  the  current  /0  is  disrupted  by  an  opening 
switch  (T  =  1 );  the  PTW  (  V+  )  arriving  at  the  switch  end 
is  reflected  back  with  equal  amplitude  and  same  sign 
(  V~  =  jZq/o)  as  shown  in  Fig.  5(d).  This  newly  generated 
NTW  ( V~  )  arrives  at  the  midpoint  and  “sees”  the  dis¬ 
continuity;  for  a  matched  load  ( R  =  \ZQ)  the  effective  re¬ 
sistance  and  the  reflection  coefficient  at  this  point  are  Zq/3 
and  T  =  —  j,  respectively.  The  resulting  reflected  and 
transmitted  waves  are  shown  in  Fig.  5(e).  At  this  point, 
the  voltage  across  the  load  resistance  due  to  this  new  NTW 
is  of  the  same  sign  as  that  of  the  PTW  from  the  shorted 
end,  thus  generating  the  output  voltage  and  current 
Kout  =  /oZq  and  =  I0  until  time  r  =  2f/v  when  both  the 
reflected  wave  (towards  the  switch)  and  the  transmitted 
wave  (towards  the  shorted  end)  reflect  back  again  to  the 
midpoint;  after  this  time,  both  voltage  and  current  are 
completely  canceled  along  the  entire  length  of  both  trans¬ 
mission  lines  {see  Fig.  5(f)].  It  is  straightforward  to  show 
that  for  the  cases  with  mismatched  loads  ( R^={Z0)  a  senes 
of  postpulses  of  gradually  decreasing  amplitude  with  either 
the  same  sign  or  alternating  sign  (depending  on  the  con¬ 
dition  R  >  {Zq  or  R  <  jZq)  are  generated,  as  shown  in  Fig. 
6.  The  amplitude  of  the  nth  postpulse  (again  n  =  0  corre¬ 
sponds  to  the  main  pulse )  is  given  by 


V(xi 


r=.os  r.i 


FIG.  S.  Dual  of  the  Btumiem  line:  (a)  schematic  diagram;  spatial  distri¬ 
bution  of  (b)  current  and  (c)  voltage,  shown  as  the  superposition  of  a 
positively  traveling  wave  and  a  negatively  traveling  wave  before  the 
switch 's  opened  (t< 0);  and  spatial  distributions  of  voltage  during  the 
periods,  (d)  0<  t  <  f/v.  (e)  t/v  <  t  <  If/ u,  and  (f)  2 f/v  <  l  <  i//v.  respec¬ 
tively. 

roat  =  2RZ0I0nZ0-2R)y(Z0+2R)n^].  (8) 

Note  that,  in  this  case,  each  pulse  is  separated  by  time 


V(U 

4  n=0 


0  0  5  t  1,5: 


V(ti 

*  n=0 


0  5:  1.5:  2.5:  3  5  :  4  5:  5.5: 


FIG.  6.  Typical  output  waveforms  of  the  dual  of  the  Blumiein  line  with 
(a)  matched  load  ( R  =  5Z0).  (b)  R  >jZo  (R  =  4,).  and  ( c )  R  <:Z„ 
1/1  =  4/  4) 
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2f/v.  This  is  in  contrast  to  the  current-charged 
transmission-line  case  in  which  no  such  interpulse  separa¬ 
tion  exists,  as  described  previously. 

As  with  its  dual  counterpart,  the  Blumlein  line,  this 
system  is  capable  of  producing  a  square  pulse  of  equal 
amplitude  to  the  charging  current;  this  is  twice  that  of  the 
simple  current  charged  transmission-line  case.  Another  no¬ 
table  feature  of  this  system  is  that  the  opening  switch  is 
positioned  at  the  end  of  the  transmission  line;  this  allows 
direct  connection  of  the  load  to  the  output  of  the  system, 
thus  removing  possible  circuit  interference  by  the  switch. 

VI.  CONCLUSIONS 

We  have  systematically  analyzed  pulse  formation  in 
basic  inductive-energy-storage  PPS  circuits.  Exact  expres¬ 
sions  for  output  waveforms  with  relevant  circuit  parame¬ 
ters  are  derived  for  four  basic  circuits:  the  lumped  inductor 
system,  the  current  charged  transmission  line,  the  dual  of 
the  LC  generator,  and  the  dual  of  the  Blumlein  line.  The 
last  two  systems  are  newly  proposed  and  may  be  useful  for 
practical  applications. 
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APPENDIX:  DERIVATION  OF  EQ.  (7) 

For  the  circuit  of  Fig.  4,  using  KirchhofFs  current  law 
at  the  node  where  v  is  measured,  we  obtain 


This  is  a  standard  cubic  equation  having  three  roots — z,t 
z2,  and  z3.  The  calculation  of  these  roots  may  be  found  in 
any  mathematical  handbook.  In  this  case,  one  real  and  two 
complex  conjugate  roots  are  found  as 

z{=—a,  z2=  —  (5  +  ico,  z}=  —  0  —  ico,  (A5) 

where  a  =  2R(3L)~l  —  (r  +  d)'n  —  (r  —  d)wi , 
0  =  2 -'[(/•  +  d)Wi  +  (r  -  d)l/}]  +  2R(3L)-\ 
co  =  31/22  ~'[(r  +  d)wi  -  (r  -  d),/3],  d  =  (q2 
+  r2)172.  q  =  (3 LC)~]  -  (2R)2(3L)-\  and  r 
=  —  /?(6Z.2C)  -  1  -  (2/?)3(3Z.)  ~  3. 

The  final  form  of  the  solution  to  Eq.  (A3)  may  now  be 
wntten  as 

v(t)  =Ae~at  -f  e  ~8'{B  cos  cot  4-  C  sin  cot).  (A6) 

To  evaluate  the  unknown  constants  A.  B,  and  C,  we  need 
three  initial  conditions.  These  initial  conditions  are  found 
as  follows 

(i)  initial  voltage  u(0)  =0; 

(ii)  taking  the  derivative  of  Eq.  (Al),  eliminating  the 
term  di/dt  in  Eq.  (A2),  and  setting  u(0)  =0,  we  obtain 
dv/dt(0)  =  0;  and 

(iii)  taking  the  second  derivative  of  Eq.  (Al )  and  the 
derivative  of  Eq.  (A2),  eliminating  the  term  d2ic/dt1,  and 
setting  u(0)  =  0,  dv/dt{0)  =  0,  and  ic(0)  =  /q,  we  obtain 
d2v/dt2(0)  =  RIq/LC. 

We  now  obtain  the  final  solution  using  the  above  initial 
conditions  as 


V  1  n 


v  dt  4-  it  —  0, 


where  ic  is  the  charging  current  for  the  capacitor.  The 
voltage  v  may  be  expressed  in  terms  of  ic  as 


dic  1  C 

v=Ld,^c)«‘<d‘ 


Elimination  of  Ic  in  Eqs.  ( Al )  and  ( A2)  yields  the  follow¬ 
ing  third-order  differential  equation: 

d-v  2  Rdh  1  dv  Rv 

d?  +  Td?  +  Lc7t  +  Uc=0-  (A3) 

Assuming  the  exponential  solution  v  =  7.A„  exp (z„f), 
the  characteristic  equation  for  Eq.  (A3)  is  wntten  as 


X(e  —  e  e'{coswf  —  [(a  —  0)/a>]sin <wf}). 

(A7) 
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EXPERIMENTAL  INVESTIGATION  OF 
PSEUDOSPARK-PRODUCED  ELECTRON  BEAMS 


K.  K.  Jain  and  M.  J.  Rhee 
Laboratory  for  Plasma  Research 
University  of  Maryland 
College  Park,  MD  20742  U.  S.  A. 

Abstract 

The  pseudospark  discharge  which  produces  pulsed  electron  beams  of  inter¬ 
esting  characteristics  has  been  studied.  The  multigap  pseudospark  chamber 
filled  with  argon  gas  is  operated  at  voltages  up  to  50  kV.  The  electron  beams 
of  peak  current  ~150  A  and  pulse  duration  ~20  ns  for  Cext  =  980  pF  at  repet¬ 
itive  frequency  up  to  1  kHz  are  extracted.  A  study  of  the  empirical  scaling 
law  for  the  electron  beam  with  various  experimental  parameters  is  carried  out. 

The  normalized  rms  emittance  of  the  electron  beam,  measured  by  a  slit-hole 
technique,  is  ~15  mm-mrad,  yielding  a  normalized  brightness  of  Bn  ~  2x  1010 
A/m2rad2. 

1  INTRODUCTION 

In  recent  years,  pseudospark  discharge  has  received  a  significant  attention  because  of 
its  ability  to  sustain  anomalous  current  density  in  low  pressure  region  of  the  Paschen 
curve1*,  gas  breakdown  (conduction)  in  nanosecond  time  scale2)  and  production  of 
high  brightness  electron  beams.3,4)  As  a  result,  the  pseudospark  opens  a  wide  field  of 
applications  such  as  in  accelerator,  free  electron  laser  and  high  voltage,  high  current 
switch  technology. 

The  pseudospark  discharge  may  be  characterized  as  a  gas  discharge  between  a 
planar  anode  and  a  hollow  cathode  in  the  lower  pressure  side  of  the  characteristic 
breakdown  curve,  which  is  very  similar  to  the  Paschen  curve  for  parallel  electrodes. 
The  pseudospark  chamber  consists  of  a  hollow  cathode  separated  from  a  planar 
anode  by  an  insulator  in  the  outer  region.  The  key  features  of  the  pseudospark  are 
that  although  it  operates  in  low  pressure  region  (unlike  spark),  gas  breakdown  takes 
place  in  fast  time  scale  and  sustains  large  current  density,  like  in  a  spark.  During 
the  breakdown  phase,  generation  of  a  short  duration,  pinched  electron  beam  has 
been  observed. 


PLEXIGLAS  INSULATORS 


Figure  1:  Schematic  of  the  experimental  setup 

We  present  here  experimental  results  of  an  investigation  of  electron  beams  pro¬ 
cured  by  a  pseudospark  with  reference  to  its  ability  to  produce  high  current,  high 
quality  beams.  Discharge  and  electron  beam  parameters  and  their  scaling  behavior 
are  determined. 

2  EXPERIMENTAL  SYSTEM 

The  experimental  setup3,4*  is  schematically  shown  in  Fig.  1.  The  pseudospark  sys¬ 
tem  used  in  this  study  is  a  modular,  multigap  discharge  chamber,  consists  of  a  hollow 
cathode,  intermediate  brass  electrodes  and  insulating  plexiglas  disks,  and  an  anode. 
The  intermediate  electrodes  have  a  center  hole  of  3.2  mm  diam  which  provides  a 
channel  for  the  discharge  and  also  passage  for  the  electron  beam.  The  electron  beam 
is  extracted  through  anode  hole  of  3.2  mm.  Experiments  have  been  carried  out  with 
1,  3,  6,  and  10  gap  systems.  Negative  high  voltages  up  to  50  kV  are  applied  to  the 
cathode  via  a  200  Mfi  current  limiting  resistor.  In  addition  to  the  self-capacitance  of 
the  chamber,  the  total  discharge  capacitance  was  increased  by  adding  various  values 
of  external  capacitance  to  investigate  the  discharge  and  beam  characteristics  with 
different  stored  energies.  A  6.3  cm  diam,  50  cm  long  drift  chamber  is  attached  to 


the  anode  side  of  the  discharge  chamber  to  accomodate  diagnostics  such  as  a  Fara- 


Pressure  (  mTorr  > 


Figure  2:  Measured  characteristic  breakdown  curve  for  four  different  gap  systems. 

day  cup  and  an  emittance  meter.  The  pseudospark  chamber  and  the  drift  tube  are 
evacuated  to  less  than  10“5  Torr  pressure  by  a  diffusion  pump.  The  argon  working 
gas  is  fed  through  a  needle  valve  and  nearly  static  gas  pressure  is  maintained  in  the 
experimental  system  by  adjusting  the  needle  valve  and  the  throttle  vacuum  valve. 
The  operating  gas  pressure  range  is  between  20  and  250  mTorr. 

Two  calibrated  miniature  Rogowski  coils  around  the  center  hole  are  built  into 
the  anode  flange.  The  coil  facing  the  cathode  side  (upstream)  measures  the  total 
discharge  current  while  the  other  coil  (downstream)  measures  electron  beam  current 
ejected  into  the  drift  chamber.  The  cathode  voltage  is  measured  by  a  capacitively 
coupled  resistive  voltage  divider.  A  movable  100  mft  Faraday  cup  of  5  cm  diam  is 
inserted  on-axis  to  determine  the  electron  beam  current  at  various  axial  positions. 
The  rms  emittance  of  the  electron  beam  is  obtained  using  a  slit-hole  type  emittance 
meter.  It  consists  of  an  array  of  seven  parallel  slits  of  300  /im  width  and  a  50  n m 
thick  radiachromic  film  as  a  beam  detector,  placed  47  mm  downstream  of  the  slit 
plane.  The  slits  sample  the  beam  at  the  given  slit  locations  producing  flat  beamlets. 
The  beamlets  are  then  allowed  to  disperse  with  their  transverse  velocity  after  passing 
through  the  slits.  The  spatial  current  density  profiles  of  dispersed  beamlets  at  the 
detector  plane  are  recorded  on  the  radiachromic  film,  from  which  the  rms  emittance 
of  the  electron  beam  is  analyzed. 
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Figure  3:  Temporal  evolution  of  the  six-gap  pseudospark  discharge:  (a)  discharge 
current  (upstream  coil)  356  A/div,  (b)  Faraday  cup  output  (1  cm  away  from  the 
anode)  50  A/div,  (c)  injected  electron  beam  current  (dowmstream  coil)  32  A/div, 
(d)  cathode  voltage  signal  10  kV/div. 

3  EXPERIMENTAL  RESULTS  AND  DISCUSSIONS 

The  breakdown  voltage  characteristics  of  the  pseudospark  chamber  is  determined 
for  various  number  of  gaps.  In  Fig.  2  the  breakdown  voltage  Vb  is  plotted  as  a 
furcticn  of  pressure  for  four  differed,  number  of  gaps.  It  is  evident  from  this  result 
that  the  breakdown  voltage  js  an  approxiamte  function  of  pd,  V*  «  }{pd),  where  d. 
which  is  proportional  to  the  number  of  gaps,  is  the  distance  between  the  cathode 
and  anode.  Repetitive  pulse  discharges  are  easily  attained  with  a  time  period  r  of 
which  is  approximately  determined  by  r  =  RC£n{\ o/(Fo  ~  14)).  By  choosing  an 
RC  value  properly,  periodic  pulse  discharges  with  a  frequency  of  1  Hz  to  1  kHz  are 
obtained. 

The  temporal  behavior  of  the  single  pulse  discharge  of  the  six-gap  pseudospark 
chamber  with  Ctxl  =  380  pF  and  14  =  —  24  kV  is  shown  in  Fig.  3.  The  discharge 
process  occurs  in  two  stages.  Up  to  ~15  ns,  the  cathode  voltage  drops  slowly,  and 
the  currents  measured  by  Rogowski  coils  and  Faraday  cup  are  small.  At  ~15  ns,  the 
main  breakdown  occurs.  The  cathode  voltage  drops  to  almost  zero  and  the  currents 
rise  very  fast  to  large  values  within  a  fast  time  scale.  The  discharge  current  shows 
a  damped  oscillating  waveform  with  a  peak  value  of  570  A.  The  measured  electron 
beam  current  is  ~  100  A  in  Fig.  3.  The  beam  current  density,  estimated  assuming 
beam  diam  equals  to  the  anode  hole  opening,  is  1.4  kA/cm2.  Typically,  10  to  20% 
of  the  peak  discharge  current  appears  as  the  electron  beam  current.  The  beam- 
induced  plasma  current  appears  to  be  contributing  to  the  dowmstream  Rogowski 


Breakdown  Voltage  (kV) 


Figure  4:  Dependence  of  electron  beam  current  measured  by  the  downstream  Ro- 
gowski  coil  and  Faraday  cup  with  breakdown  voltage  for  six-gap  pseudospark  ( Cezl 
=  380  pF). 

coil  and  Faraday  cup  signals  later  in  time,  thus  leading  to  their  slow  decay. 

The  dependence  of  the  electron  beam  current  measured  by  the  downstream  Ro- 
gowski  coil  as  well  as  Faraday  cup  (located  at  1  cm  away  from  the  anode  )  as  a 
function  of  breakdown  voltage  is  illustrated  in  Fig.  4  for  Cext  =  380  pF  and  for 
six-gap  chamber.  The  electron  beam  current  increases  linearly  with  the  breakdown 
voltage  up  to  25  kV  and  has  a  slope  of  ~  5  A/kV.  For  a  fixed  breakdown  voltage, 
the  electron  beam  current  increases  with  Cex<  and  the  beam  currents  up  to  a  few 
kA  are  extracted  out  of  the  anode  hole.  Higher  values  of  Cext  also  lead  to  increase 
of  duration  of  the  electron  beam  current  pulse. 

The  rms  emittance5)  (or  effective  emittance),  which  is  a  measure  of  beam  quality, 
is  determined  using  a  slit-hole  type  emittance  meter  with  the  analysis  technique 
described  in  Refs.  4  and  6.  The  rms  emittance  of  the  electron  beam  is  found  to  be 
t  ~  65  mm-mrad.  The  calculated  normalized  emittance  tn  =  /3-ye  at  beam  energy  of 
20  keV  is  ~  18  mm-mrad.  Thus,  the  normalized  beam  brightness  can  be  obtained 
from  the  expression  Bn  =  2//7T2e2  which  is  ~  2  x  lO10  A/m2rad2  and  compares 
favorably  with  that  of  other  high  brightness  electron  beam  sources. 

4  CONCLUSIONS 

The  behavior  of  a  multigap  pseudospark  system  and  electron  beam  produced  are 
studied.  The  pseudospark  discharge  occurs  reproducibly  at  the  pressure  and  voitage 


defined  by  the  lower  pressure  side  of  the  characteristic  breakdown  curve.  The  break¬ 
down  voltage  is  a  function  of  pd.  where  d  is  the  distance  between  the  cathode  and 
anode.  During  the  breakdown  phase,  the  pseudospark  ejects  high  current  density 
electron  beams.  An  electron  beam  of  average  energy  of  20  keV.  pulse  duration  ~10 
ns  and  of  peak  current  ~  1 50  A  is  obtained  for  Cext=  980  pF.  The  beam  current 
increases  with  the  breakdown  voltage  up  to  25  kV  and  with  the  external  capaci¬ 
tance.  Approximately  20%  of  the  stored  electrical  energy  appears  as  the  total  beam 
energy.  The  time-integrated  rms  omittance  of  the  electron  beam,  measured  using 
a  slit-hole  method,  is  ~  65  mm-mrad.  The  normalized  brightness  of  the  beam  is 
estimated  to  be  Bn  ~  2  x  1010  A/m2rad2.  which  compares  favorably  with  that  of 
existing  high-brightness  sources.  Thus,  the  pseudospark  appears  to  be  a  simple  , 
compact,  high-brightness  electron  beam  source. 
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ABSTRACT 

A  high-brightness  electron  beam  produced  by  a  six-gap  pseudospark 
chamber  operated  in  5  -  25  kV  voltage  and  30  -  70  raTorr  pressure 
ranges  is  experimentally  investigated.  The  electron  beam  of 
current  150  A  at  an  average  energy  of  20  kev  and  pulse  duration 
of  10  ns  is  extracted  with  a  repetitive  frequency  of  -10  Hz.  The 
electron  beam  current  scales  linearly  with  the  breakdown  voltage 
and  about  20  %  of  the  stored  energy  is  converted  into  the  total 
beam  energy.  The  time-integrated  rms  emittance  of  the  electron 
beam  is  measured,  and  a  typical  value  is  found  to  be  65  nun  mrad, 
yielding  a  normalized  brightness  of  the  beam  Bn  ~  2  x  10iO 
A/(m‘rad2 ) . 

PACS  numbers:  07.77.+p,  41.80.Dd,  52. 80. Mg 

a)  Present  address:  Max-Planck  Institut  fur  Quantenoptik , 

D-8046  Garching;  GSI-Darmstadt ,  D-6100,  Darmstadt,  F.R.G. 
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INTRODUCTION 


I . 


High-current  and  low  transverse-energy-spread  electron-beam 
sources  are  required  in  advanced  accelerators  such  as  e~e-  linear 
colliders1-  and  f ree-electron  lasers.2'3  Thermionic  cathodes 
produce  electron  beams  of  relatively  low  transverse  velocity- 
spread;  however,  the  current  density  is  limited*  to  less  than  100 
A/cm2  and  nonuniformity  of  the  charge  distribution  may  cause 
significant  eraittance  growth.5  On  the  other  hand,  field  emission 
diodes  are  capable  of  producing  high  current-density  electron 
beams  but  usually  have  high  energy  spread.  Strong  magnetic 
fields  are  often  used  to  guide  such  a  beam*  to  reduce  the  beam 
spreading.  Use  of  strong  magnetic  fields,  however,  has 
disadvantages,7  such  as  rippling  of  the  beam  surface,  undesirable 
beam  instabilities,  etc. 

The  limitations  of  presently  available  sources  have 
motivated  a  search  for  new  types  of  electron  beam  sources  such  as 
photocathodes9  that  can  produce  a  beam  of  high  brightness.  In 
recent  years , 1 1  the  pseudospark  device  has  been  found  to  be 
capable  of  generating  electron  beams  of  high  current  density 
(  >  1  kA/cm‘)  and  high-brightness.  The  breakdown  voltage  between 
two  parallel  electrodes  as  a  function  of  gas  pressure  is 
described  by  the  well-known  Paschen  curve.  However,  if  the 
planar  cathode  is  replaced  by  a  hollow  cathode,  the  discharge 
characteristics  are  changed;  the  breakdown  voltage  curve  is 
slightly  modified  from  the  Paschen  curve  in  the  region  between 
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the  Paschen  minimum  and  vacuum  breakdown.11-13  The  discharge 
mechanism,,  in  this  region,  is  different  from  that  of  a  high 
pressure  spark,  but  the  breakdown  time  is  comparable;  this  is  why 
it  is  called  a  "pseudospark."  A  single-gap  pseudospark  chamber 
consists  of  a  hollow  cathode  and  an  anode.  Normally,  for 
breakdown  in  the  lower  pressure  side  of  the  Paschen  minimum,  the 
longer  discharge  path  is  preferred,  which  in  the  case  of 
pseudospark  geometry  is  along  the  axis.  The  hold-off  voltage  of 
a  single-gap  is  limited  to  -50  kV  due  to  onset  of  surface 
f lashover ; 14  however  this  can  be  increased  by  additional  stacks 
of  intermediate  electrodes  and  insulator  disks  between  the  anode 
and  the  cathode.  Thus,  a  multi-gap  pseudospark  chamber  can  be 
operated  at  higher  voltages. 

During  the  breakdown  phase,  the  generation  of  a  short 
duration,  pinched  electron  beam  has  been  observed.  Until  now, 
the  understanding  of  the  mechanism  of  electron  beam  generation  has 
been  merely  gualitative.  During  the  pre-breakdown  phase,  electrons 
and  ions  are  created  by  ionization.  The  electrons  are 
immediately  collected  by  the  anode  while  the  positive  charges 
move  toward  the  cathode  and  mainly  into  the  hollow  region  of  the 
cathode,  forming  a  virtual  anode.  In  the  breakdown  phase,  a  fast 
carrier  multiplication  process3  inside  the  virtual  anode  leads  to 
a  high  plasma  density,  which  is  followed  by  a  very  fast 
breakdown,  discharging  the  system.  The  field  distribution 
between  virtual  anode,  cathode,  and  anode  favors  the  formation  of 
a  pinched  electron  beam13  directed  toward  the  anode. 
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In  this  paper,  we  report  results  of  an  investigation  of  an 


electron  beam  produced  by  a  pseudospark  with  reference  to  its 
ability  to  produce  high-current,  high-quality  beams.  Detailed 
experimental  studies  are  performed  to  determine  the  electrical 
parameters  of  the  beam  and  its  scaling  behavior.  The  rms 
emittance  of  the  electron  beam  is  measured  and  then  the 
brightness  is  analyzed.  In  Sec.  II,  the  experimental  device  and 
associated  diagnostics  used  are  described.  Results  and 
discussions  are  presented  in  Sec.  III.  The  main  results  of  the 
experimental  study  are  concluded  in  the  last  section. 


II.  EXPERIMENTAL  ARRANGEMENT 

Figure  1  shows  a  schematic  of  our  experimental  setup.10  The 
pseudospark  chamber  consists  of  an  easy-to-change ,  O-ring  seal, 
modular,  six-gap  electrodes,  and  insulator  disks  stacked  between  a 
hollow  cathode  and  an  anode.  The  anode  has  a  3.2  mm-diam,  on- 
axis  hole  for  the  beam  extraction.  Intermediate  electrodes  made 
of  3.2  mm-thick  brass  have  outer  diara  of  6.35  cm  and  a  center 
hole  of  diam  3.2  mm.  The  3.2  mm  thick  plexiglas  insulator  disks 
have  outer  and  inner  diam's  of  7.0  cm  and  2.54  cm,  respectively. 

In  addition  to  the  self  capacitance  of  11  pF  of  the  chamber 
between  cathode  and  anode,  a  low  inductance,  external  capacitance 
of  380  pF,  760  pF,  or  980  pF  is  added  to  investigate  the  beam 
behavior  with  different  stored  energies.  Negative  high  voltage 
up  to  30  kV  is  applied  to  the  hollow  cathode  via  a  200  Mfl 
charging  resistor.  A  50  cm-long  drift  tube  is  attached  to  the 
anode  side  for  beam  propagation  and  to  accomodate  diagnostics, 
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such  as  a  Faraday  cup  and  an  emittance  meter. 

A  mechanical  pump  evacuates  the  experimental  system  down  to 
about  1  mTorr.  The  pressure  is  measured  by  a  capacitance 
manometer  type  vacuum  gauge.  The  argon  working  gas  is  fed 
through  a  needle  valve,  from  the  cathode  side,  into  the  system. 
The  operating  pressure  range  is  between  30  to  70  mTorr.  Two 
calibrated  Rogowski  coils  are  molded  into  axisymmetric  grooves 
milled  into  both  sides  of  the  anode  flange.  The  coil  facing  the 
cathode  side  (upstream  coil)  measures  the  total  discharge 
current,  while  the  other  coil  facing  the  drift  chamber 
(downstream  coil)  measures  the  electron  beam  current  (or  net 
current)  passing  through  the  anode  hole  into  the  drift  chamber. 
The  cathode  voltage  is  measured  by  a  high- impedance  resistive 
voltage  probe.  A  movable  100  mfl  Faraday  cup  of  5  cm-diam 
graphite  is  inserted  on-axis  through  the  end  flange  of  the  drift 
tube  to  measure  the  electron  beam  current  at  various  axial 
positions  as  it  propagates  in  the  drift  tube. 

The  rms  emittance  of  the  pseudospark-produced  electron  beam 
is  measured  using  a  slit-hole  type  emittance  meter.15  It 
consists  of  a  series  of  parallel  thin  slits  of  300  urn  width  and  2 
mm  spacing,  constructed  from  a  0.6  mm  thick  stainless  steel 
plate.  These  slits  produce  sheet  bearalets  which  are  detected  by 
a  radiachromic  film16  placed  47  ram  farther  downstream.  The  film 
material  gives  a  linear  response  up  to  an  absorbed  dose  of  10a 
rad.  By  use  of  an  applied  magnetic  field  it  is  found  that  the 
response  pattern  on  the  film  is  completely  swept  away,  assuring 
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that  UV  or  other  light,  for  which  the  film  is  also  sensitive, 
produced  in  the  experiment  does  not  contribute  to  the  film 
response.  The  resulting  spatial  density  distribution  on  the  film 
is  scanned  by  an  optical  microdensitometer  and  the  emittance  is 
analyzed . 


III.  MEASUREMENTS,  RESULTS,  AND  DISCUSSION 

A.  Breakdown  and  beam  characteristics 

The  electrical  breakdown  characteristic  of  the  pseudospark 
chamber  is  measured  by  fixing  the  charging  voltage  to  a  certain 
value,  and  slowly  increasing  the  gas  pressure,  uriL.il  breakdown 
occurs.  Nearly  static  gas  pressure  is  maintained  in  the 
experimental  system  by  adjusting  the  needle  valve  and  throttle 
vacuum  valve,  as  shown  in  Fig.  1.  The  measured  breakdown-voltage 
Vto  is  plotted  as  a  function  of  pressure  in  Fig.  2,  which  shows 
that  the  breakdown  voltage  is  a  sensitive  function  of  pressure  in 
the  30  -  45  mTorr  range.  Having  an  RC  time  constant  of  a 
fraction  of  a  second,  a  periodic  pulse  discharge  at  1  -  15  Hz  can 
be  easily  attained.  The  frequency  of  discharge  can  be  adjusted 
by  either  changing  the  charging  voltage  V0,  the  gas  pressure,  or 
both . 


The  temporal  evolution  of  the  discharge  is  shown  in  Fig.  3 
for  a  breakdown  voltage  of  Vto  =  -24  kV  and  C.xt  =  380  pF.  All  the 
signals  show  that  the  discharge  process  occurs  in  two  stages. 

For  the  first  15  ns,  the  discharge  voltage  drops  slowly  from  -24 
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kV  to  -22  kV,  and  the  discharge  and  electron  beam  currents  are 
small  (less  than  10  %  of  the  peak  value).  At  around  15  ns,  the 
main  breakdown  occurs.  The  discharge  voltage  drops  from  -22  kV 
to  almost  zero  in  the  next  20  ns.  During  this  period,  the 
discharge  and  the  electron  beam  current  start  rising  to  a  large 
value.  The  discharge  current  shows  a  damped  oscillating 
waveform.  The  downstream  Rogowski  coil  and  the  Faraday  cup 
(located  at  1-cm  away  from  the  anode)  register  only  one  negative 
peak  at  the  time  of  discharge  current  maximum.  The  peak 
discharge  current  in  Fig.  3(a)  is  570  A.  The  measured  electron 
beam  current  is  -100  A  at  Vto  =  24  kV  and  Cmlcx,  =  380  pF.  The  beam 
current  density  at  the  source  exit  is  estimated:  assuming  a  beam 
diara  of  3.2  mm  (  anode  hole  opening),  the  observed  beam  current 
of  100  A  yields  the  beam  current  density  of  1.4  kA/cm2 . 

Typically,  10  to  20  %  of  the  total  discharge  current  appears  as 
electron  beam  current.  Perhaps  beam-induced  plasma  current  is 
contributing  to  the  downstream  Rogowski  coil  and  Faraday  cup 
signals  later  in  time,  thus  leading  to  their  slow  decay. 
Measurements  of  the  beam  current  as  it  propagates  in  the  drift 
tube  are  done  with  a  movable  Faraday  cup.  In  Fig.  4(a),  the  peak 
value  of  the  beam  current  is  plotted  as  a  function  of  axial 
distance  from  the  anode,  showing  exponential  decay. 

The  dependence  of  electron  beam  current  as  a  function  of 
breakdown  voltage  is  determined  and  shown  in  Fig.  5  for  C.xt  - 
380  pF.  The  electron  beam  current  increases  linearly  with  the 
voltage  and  has  a  slope  of  -5  A/kV.  It  is  observed  that  the 
higher  values  of  lead  to  increase  in  both  the  duration  and 
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amplitude  of  the  beam  current.  The  beam  current  roughly  scales 


as  the  square  root  of  The  beam  current  scales  linearly 

with  the  breakdown  voltage  for  all  values  of  C.^'s  used.  The 
efficiency  of  the  electron  beam  generation  by  the  pseudospark  may 
be  defined  here  as 

-  =  f  vtoi  dt/hcmxxy»,  (1) 

and  is  found  to  be  nearly  20  %  in  the  case  of  =  380  pF . 

B.  Emittance  measurement 

The  rms  emittance17'1®  or  effective  emittance,5''  defined 
as 


e  =  4 ( <x2  ><x ' 2 >  -  <xx'>2)1/2,  (2) 

is  a  figure  of  merit  for  the  beam  quality.  In  Eq.  (2),  x'  is  the 
gradient  of  the  particle  trajectory  given  by  x'  =  dx/dz  =  P^/P* 
and  the  angular  brackets  denote  average  values  over  the  two- 
dimensional  phase-space  as 


<  <p  > 


o  (  x  ,  x ' )  dxdx ' 


(3) 


where  . (x,x')  is  the  projected  density  in  the  two-dimensional 
phase-space.  The  second  moments  of  this  distribution,  <x2>. 
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<x,=>,  and  <xx'>,  are  related  to  the  beam  width,  the  velocity 
spread,  and  the  beam  divergence,  respectively. 

Various  experimental  methods18'20  have  been  used  to  measure 
the  beam  emittance,  such  as  two-slit,  hole-slit,  and  hole-hole. 
Under  the  assumption  of  a  Maxwellian  transverse  velocity 
distribution  and  an  axisymmetric ,  nonrotational  beam,  a  simple 
slit-hole  type  emittance  meter  has  recently  been  suggested. 15 

In  our  experiment,  we  have  used  a  slit-hole  type  emittance 
meter  with  the  analysis  technique  described  in  Ref.  15  to 
determine  the  pseudospark-produced  electron  beam  emittance.  The 
electron  beam  is  allowed  to  pass  through  an  array  of  seven  thin 
slits.  The  slits  are  long  (in  the  y  direction)  as  compared  to 
the  beam  dimension  and  thin  (in  the  x  direction).  Thus,  the 
slits  integrate  over  all  velocity  distributions  in  the  y 
direction  while  the  velocity  distribution  in  the  x  direction  is 
determined.  The  transmitted  sheet  beamlets  formed  by  the  slits 
are  then  allowed  to  traverse  a  field-free  region  of  4.7  cm  in 
length  before  they  strike  a  thin  radiachromic  film.  \  shutter 
(which  is  not  shown  in  Fig.  l)  blocks  the  beam  path  after  a  few 
discharges  when  the  darkening  profile  on  the  film  due  to  the  beam 
is  just  visible:  the  peak  value  of  optical  density  of  the 
profile  is  kept  under  0.5,  which  insures  that  the  film  response 
is  linear. 

To  obtain  the  rms  emittance  of  the  beam,  we  scan  the  optical 
density  profile  of  sheet  beamlets  as  recorded  on  the  film  along 
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the  line  of  y  =  0  by  a  microdensitometer.  A  typical  measured 
density  profile  as  a  function  of  radial  position  r  (or  x)  is 
shown  in  Fig.  6.  In  order  to  simplify  the  emittance  analysis  as 
described  in  Ref.  15,  the  measured  density  profile  of  each 
beamlet  is  approximated  as  a  Maxwellian.  Ideally,  infinitely 
thin  slits  should  be  used  in  the  emittance  meter  to  directly 
obtain  the  mean  diverging  angle  a1(  the  peak  density  &±,  and  the 
rms  width  a ±  of  the  individual  profile,  where  i  denotes  the  ith 
slit.  In  the  experiment,  however,  a  finite  slit  width  of  0.3  mm 
is  used  to  obtain  the  beamlet 's  profile,  which  is  subject  to 
deviation  from  the  true  profile  of  the  beam  distribution.  We 
have  derived  a  method  of  correcting  the  rms  width  and  tho  teak 
height  of  the  density  profile  produced  by  the  finite-width  slit. 
In  this  method,  the  density  profile  on  the  detector  plane  is 
considered  as  a  superposition  of  profiles  of  beamlets  which  are 
produced  by  infinitely  thin  slits  distributed  continuously  across 
the  finite  width  of  the  slit.  If  one  assumes  that  the  beam  has  a 
Maxwellian  velocity  distribution  and  its  intensity  remains 
constant  (or  a  very  weak  function  of  x)  within  the  slit  width  2w, 
the  resulting  profile  for  the  ith  sheet  beamlet  is  found  as 

'  C  x  w  >  /  l 

Di(x)  =  Hi  expl-^x'/aj2]  dx' 

J  <  X-W ^  /L 

X+W  x*W 

=  (IT/2)  1/2Bia![erf  (  — )  -  erf( - )],  (4) 

3 1  Si 

where  x  is  the  distance  measured  from  the  ith  slit  position, 
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erf(x)  =  2/fn  j^expC -u2  )  du ,  and  st  =  /2atL.  It  is  noted 

that  for  an  ideal  case  as  w  ->  0,  Eq.  (4)  reduces  to 
2w6iL  Lexp[ ( x '/ffi ) 2 ] /  which  is  2w/L  times  the  original 
distribution.  Given  the  peak  value  of  measured  profile  Dt(0), 
one  can  find  the  as 

Bi  =  D  (  0  )  /  [  (  2n )  ^  y  2  a  erf(— — - )].  (5) 

J2a  tL 

The  difference  of  the  two  error  functions  inside  the  bracket  [  j 
in  Eq.  (4)  determines  the  shape  of  the  resulting  profile  as  a 
function  of  x  for  two  parameters  w  and  st  =  fla y L.  The  full- 
width  at  half -maxi mums  (FWHM's)  of  Eq.  (4)  are  numerically 
calculated  for  various  values  of  w/L  and  a.  The  a  as  a  function 
of  measured  FWHM  (normalized  to  L)  with  a  parameter  w  (normalized 
to  L)  is  plotted  as  shown  in  Fig.  7.  Thus,  one  can  find  aL  in 
terms  of  measured  FWHM  of  the  ith  beamlet,  knowing  the  slit  width 
w  used.  The  height  of  the  beam  distribution  6t  is  then  found 
from  Eq.  (5)  using  the  obcained  value  of  at.  A  general  method 
which  can  be  applicable  to  other  beam  distributions  is  derived 
and  will  be  treated  elsewhere.  Analytical  expressions  of  a(r), 
3(r),  and  c(r)  are  found  by  curve  fitting  the  discrete  values  of 
at,  fit,  and  a L ,  respectively. 

Numerical  integrations  are  then  performed  using  the 
expressions  of  a(r),  6(r),  and  o(r)  to  find  and  the  moments  • 
<x2>,  <x'2>,  and  <xx'2>  (see  Ref.  15  for  details).  The  rms 
emittance  given  by  Eq.  (3),  with  the  obtained  moments  at  an  axial 
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distance  of  9  cm  from  the  anode,  is  found  to  be 


e  =  65  mm  mrad  .  (6) 

The  normalized  rms  emittance  =  B  v  e  (where  6  =  v/c,  ■  = 
(l-fl2)_1/2,  v  is  the  beam  velocity,  c  is  the  velocity  of  light) 
is  used  for  comparing  the  beam  quality  at  different  energies. 
With  an  average  beam  energy  of  20  keV,  the  calculated  normalized 
emittance  is  -18  mm  mrad.  Once  the  emittance  and  the  beam 
current  are  determined,  the  normalized  beam  brightness  can  be 
calculated  from  the  expression2'3'21 


For  measured  beam  current  of  32  A  at  9  cm  axial  distance 
from  the  anode  for  this  set  of  experiments, 

Br,  =  2  x  1010  A/(m2  rad2),  (8) 

which  compares  favorably  with  that  of  other  high-brightness 
sources.  Further,  assuming  a  Maxwellian  transverse  velocity 
distribution  with  an  rms  velocity  vrm.  =  (kT/m0)1/2  in  the 
transverse  direction,  the  normalized  emittance  at  the  source  can 
be  written  as  en  =  2r*  ( kT/nioC2 ) 1  /2 ,  where  r„  is  the  electron  beam 
source  radius  and  T  is  the  temperature.  Assuming  the  same  value 
of  €„  =  18  mm  mrad  at  the  source  and  with  beam  radius  of  1.6  ram, 
an  approximate  electron  beam  temperature  of  15  eV  is  obtained. 

The  electron  beam  profile  and  its  emittance  are  determined 
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at  different  axial  distances  from  the  anode  with  the  help  of  a 
movable  eraittance  meter.  By  fitting  the  optical  density  peaks  of 
the  beamlets  recorded  on  the  film  to  a  Gaussian  profile,  the  rms 
radius  of  the  beam  is  obtained  as  shown  in  Fig.  4(b).  The  beam 
radius  increases  from  2.2  mm  to  5.0  mm  during  its  propagation 
over  an  axial  distance  of  15  cm.  The  measured  dependence  of 
emittance  on  axial  distance  from  the  anode  is  plotted  in  Fig. 
4(c).  This  observed  increase  of  the  beam  emittance  as  it 
propagates  may  be  due  to  the  nonlinear  space-charge  effect. 

The  calcu) i . ed  value  of  -  from  the  measurements  is  used  to 
obtain  isodensity  contours  of  c(x,x7)  in  x-x7  space  (known  as  an 
emi trance  plot)  and  is  shown  in  Fig.  8  for  the  breakdown  voltage 
of  20  keV  and  at  an  axial  distance  of  15  cm  from  the  anode.  It 
is  interesting  to  note  from  Fig.  8  that  the  emittance  plot  shows 
slightly  S-shaped  contours  (non-elliptical)  which  is  indicative 
of  the  existance  of  nonlinear  forces  on  the  beam,  again  perhaps 
due  to  the  strong  nonlinear  space-charge  field5  of  the  beam. 


IV.  CONCLUSIONS 

We  have  studied  the  characteristics  of  a  pseudospark-produced 
electron  beam.  The  pseudospark  discharge  takes  place  at  the 
lower  pressure  side  of  the  Paschen  minimum  in  the  characteristic 
breakdown  curve.  The  main  pseudospark  discharge  occurs  after 
some  time-delay.  During  the  breakdown  phase,  the  pseudospark 
discharge  ejects  an  electron  beam  of  average  energy  20  keV,  peak 
current  150  A,  and  pulse  duration  -10  ns.  The  observed  current 
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corresponds  to  beam  current  density  greater  than  10 1  A/cm“  (at 
the  anode  exit) ,  which  is  an  extremely  high  value  at  such  a 
modest  beam  energy  of  20  keV.  Approximately  20  %  of  the  stored 
electrical  energy  appears  as  the  total  beam  energy.  The  time- 
integrated  rms  eraittance  of  an  electron  beam  is  measured  using  a 
slit-hole  method.  An  rms  emittance  of  value  -  65  mm  mrad  is 
found  at  an  axial  distance  of  9  cm  from  the  anode.  The 
normalized  brightness  of  the  beam  is  then  estimated  as 
B„  »  2  x  1010  A/(m2  rad2),  which  compares  favorably  with  that  of 
existing  high-brightness  sources.  The  electron  beam  current 
scales  linearly  with  the  breakdown  voltage  as  5  A/kV.  This  is  a 
very  encouraging  result,  suggesting  that  the  electron  beam 
brightness  may  further  increase  at  higher  voltage  operation  of 
the  pseudospark.  Thus,  the  pseudospark  appears  to  be  a  simple, 
compact,  high-brightness,  electron-beam  source. 
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FIGURE  CAPTIONS 


Fig.  1 
Fig.  2. 
Fig.  3. 


Fig.  4. 

Fig.  5 

Fig.  6 

Fig.  7 


Schematic  of  the  experimental  setup. 

Measured  characteristic  breakdown  curve. 

Temporal  evolution  of  the  pseudospark  discharge: 

(a)  discharge  current  (upstream  coil)  356  A/div, 

(b)  Faraday  cup  output  ( 1-cra  away  from  the  anode)  50 
A/div,  (c)  injected  electron  beam  current  (downstream 
Rogowski  coil)  32  A/div,  (d)  cathode  voltage  signal 
10  kV/div. 

Dependence  of  beam  parameters  on  the  axial  distance 
from  the  anode:  (a)  the  electron  beam  current,  (b)  rms 
beam  radius,  (c)  rms  emittance. 

Variation  of  electron  beam  currents  measured  by 
the  downstream  Rogowski  coil  and  Faraday  cup  with 
breakdown  voltage  (C.xt  =  380  pF) . 

Optical  density  distribution  on  the  radiachromic  film 
after  exposure  to  the  sheet  bearalets  at  a  9  cm  distance 
from  the  anode. 

Plot  of  measured  FWHM  of  the  beamlet's  optical  density 
distribution  against  the  rms  width  of  the  original 
distribution  for  various  values  of  slit-width  w.  Here 
FWHM  and  w  are  normalized  to  L,  the  distance  between  the 
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film  and  the  slit  plane. 


Fig.  8 


Two-dimensional  isodensity-contour  plot  of  phase-space 
density  - ( x , x ' )  . 
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ABSTRACT 

A  brief  review  on  the  status  of  pseudospark  research  is  given.  A  modular  pseu¬ 
dospark  chamber  was  built  to  investigate  this  new  discharge  type  as  a  source  of  intense 
electron  beams.  Typical  experimental  parameters  are  a  breakdown  voltage  of  25  kV 
and  a  discharge  current  of  about  500  A,  with  an  electron  beam  of  ~100  A  and  ~10 
ns  FWHM  emitted  through  the  anode  hole  into  a  drift  chamber.  Electrical  param¬ 
eters  of  the  circuit,  including  the  plasma  channel,  were  evaluated  by  monitoring  the 
discharge  current  waveform.  First  results  of  emittance  measurements  of  the  electron 
beam  are  presented.  Based  on  the  measured  rms  emittance  of  ~  55  mm-mrad.  a  nor¬ 
malized  brightness  of  about  2  x  10u  A/m2-rad2  is  estimated.  This  characterizes  the 
pseudospark  as  a  high  brightness  electron  beam  source. 

REVIEW  OF  PSEUDOSPARK  MECHANISM  AND  APPLICATIONS 

The  “pseudospark”  phenomenon  was  first  reported  by  Christiansen  and  Schultheiss 
in  1978*  as  a  fast,  gas  discharge  which  occurred  at  low  pressure  in  a  special  device, 
called  a  pseudospark  chamber.  The  breakdown  time  q!  this  discharge  is  found  to  be 
similar  to  that  of  a  high  pressure  spark,  although  their  mechanisms  appear  to  be  totally 
different.2-3  A  single-gap  pseudospark  chamber  consists  of  an  anode  with  a  center  hole 
and  a  hollow-cathode;  the  electrodes  are  separated  by  a  few  mm  thick  insulating  washer. 
The  pseudospark  chamber  may  be  described  as  a  combination  of  a  restricted  linear 
discharge  geometry  and  a  hollow  cathode.  The  breakdown  characteristics  depend  on 
the  operating  point  on  the  breakdown  curve:  fast  breakdown  occurs  at  higher  voltage 
(lower  pressure).  A  characteristic  breakdown  curve  similar  to  the  Paschen  curve  lor 
parallel  electrodes  is  obtained.  The  pseudospark  is  initiated  on  the  left  side  of  this 
breakdown  curve.  In  this  regime,  the  gas  discharge  occurs  along  the  longest  possible 
path  in  the  gap  which,  in  the  case  of  a  pseudospark  geometry,  is  found  on  the  axis. 
The  holdolf  capability  of  one  gap  is  limited  to  40  to  50  kV  due  to  the  onset  of  surface 
flashover  and  field  emission;  however,  this  holdofT  voltage  can  be  enhanced  bv  adding 
additional  stacks  of  intermediate  electrodes  and  insulator  discs  forming  a  muitigap 
rhamber.  Thus  the  field  is  graded  by  the  capacilively  divided  voltage,  resulting  in  a 
reduced  voltage  across  each  gap  by  the  increased  number  of  stacks. 

Although  the  pseudospark  pitenomcnon  was  discovered  more  than  a  decade  ago. 
ns  mechanisms  are  still  not  fuliv  understood.  Some  of  the  interesting  information  was 
obtained  through  spectroscopic  observations.2,3 

The  discharge  starts  with  a  high  impedance  Townsend  predischarge,  which  lasts 
for  a  few  microseconds  before  the  main  breakdown  occurs.  During  this  period  a  light 
emitting  zone  has  been  observed,  starting  from  the  anode  and  moving  towards  the 
cathode  with  a  constant  velocity  which  is  measured  to  be  1  cm/ps  in  150  mTorr  of 
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argon  eas.  The  current  rises  from  10  jrA  to  a  few  hundred  jiA  until  a  stationary 
high  voltage  glow  discharge  with  constant  current  is  achieved.  At  this  point  the  light 
emitting  zone  becomes  stationary  at  a  distance  of  a  few  mm  in  front  of  the  cathode. 
During  the  predischarge,  a  high  reduced  electrical  field  strength  E/n  of  over  10~14 
V  cm2  is  produced,  resulting  in  generation  of  high-energy  runaway  electrons.  This 
explains  the  electron  emission,  through  the  anode  hole  during  the  predischarge.  At  the 
same  time  positive  ions  drift  to  the  cathode  accumulating  a  positive  space  charge. 

When  the  positive  space  charge  on  axis  reaches  a  critical  value  at  the  cathode  hole, 
the  hollow  cathode  discharge  is  ignited  and  the  main  discharge  takes  place.  Ionization 
waves  moving  into  the  hollow  cathode  at  velocities  of  approximately  100  cm/ps  are 
observed.  The  light  intensity  in  the  gap  between  the  anode  and  cathode  rises  and  an 
intense  electron  beam  is  ejected  at  the  anode.  When  the  voltage  drops,  within  ~10  ns. 
current  rise  rates  of  up  to  8  x  1011  A/s  and  current  densities  exceeding  100  kA/cm;  are 
observed.4  These  extremely  high  current  densities  in  a  glow  discharge  are  inconsistent 
with  a  cold  cathode  emission  process,  and  are  attributed  to  field  enhanced  emission 
from  the  melted  electrode  surface.5 

It  is  possible  to  trigger  the  main  discharge  in  the  hollow  cathode  region  with  sub¬ 
nanosecond  jitter,  without  the  long  predischarge,  when  a  sulficient  amount  of  charge 
carriers  is  provided  in  the  hollow  cathode  region.  This  can  be  achieved  either  by  ini¬ 
tiating  a  surface  flash  over6,  by  providing  a  pulsed  glow  discharge',  or  by  illumination 
with  U V-light.8  First  numerical  simulation  results  of  the  pre-  and  early  main  discharge 
phases  were  reported  recently.9 

As  was  introduced  initially,  the  pseudospark  discharge  produces  well  pinched  elec¬ 
tron  beams  with  current  densities  up  to  10s  A/cm2  and  power  densities  up  to  109 
W/cin;.  Such  high-power  electron  beams  can  be  used  for  material  processing,  for  ex¬ 
ample  drilling  holes  in  metal  targets  or  evaporation  of  semiconductors  and  isolators 
to  produce  layers  of  high  temperature  superconducting  materials.10  With  the  emis¬ 
sion  of  the  mtense  electron  beam,  pulsed  microwave11  and  X-rav  generation  are  also 
observed.212  possibly  allowing  further  applications  in  these  fields  in  the  future. 

The  pseudospark  is  easily  triggered  with  high  precision  in  the  hollow  cathode  region. 
Initially,  this  behavior  favored  the  development  of  high  power  pseudospark  switches, 
first  began  at  CERN.6  Since  then,  pseudospark  switches  tor  various  applications  have 
been  developed.  Current  pulses  of  up  to  '200  kA  with  a  titter  of  —  1 00  ns  were  switched 
at  0.3  1 1  z  for  over  500.000  shots  without  significant  electrode  deterioration.13  With  a 
pseudospark  switch  designed  for  a  laser  system,  a  jitter  of  4  ns  at  25  kA  peak  current 
and  100  Hz  repetition  rate  were  measured.4  Multichannel  pseudospark  switches11  have 
achieved  current  rise  rates  of  up  to  2.4  x  1012  A/s.  with  a  peak  current  of  15  kA  and 
a  jitter  of  1  ns. 

EXPERIMENTAL  STUDIES  AT  THE  UNIVERSITY  OF  MARYLAND 

An  easv  to  change,  o-ring  sealed,  modular,  multigap  discharge  chamber  system  with 
brass  eiectrodes  and  plexiglas  discs  was  constructed.  For  the  experiments  reported  here 
a  ti  gap  pseudospark  chamber  with  a  3.2  mm  diam  center  hole  on  axis  was  used.  The 
experimental  setup  is  shown  in  Fig.  1. 
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Fig.  1.  Experimental  set  up. 

Negative  high-voltage  up  to  35  kV  is  applied  to  the  hollow-cathode  via  a  200  MQ 
current-limiting  resistor.  In  addition  to  the  self-capacitance  of  1 1  pF  of  the  chamber, 
the  total  discharge  capacitance  was  increased  by  adding  external  capacitors  of  420  pF 
and  760  pF  to  investigate  the  behavior  with  different  stored  energies.  The  chamber 
was  mounted  on  a  plexiglas  drift  tube  of  50  cm  length  and  6.3  cm  diameter.  A  brass 
flange  which  terminated  the  tube  was  connected  via  6  brass  bars  equally  spaced  on  a 
radius  of  4.4  cm  to  the  grounded  anode,  thus  providing  a  current  return  path. 

A  two-stage  mechanical  pump  evacuated  the  assembly  to  about  1  mTorr.  The 
working  gas  was  usually  argon,  which  was  injected  through  a  needle  valve  from  the 
cathode  side  into  the  system.  The  operating  pressure  range  varied  from  20  to  60 
mTorr. 

Two  calibrated  miniature  Rogowski  coils  around  the  center  hole  were  built  into  the 
anode  flange  as  can  be  seen  in  Fig.  1.  The  coil  facing  the  cathode  side  (upstream  coil) 
measured  the  total  current  in  the  discharge  channel.  The  other  coil  facing  the  drift 
chamber  (downstream  coil)  measured  the  current  which  passed  througn  the  anode  hole 
into  the  drift  chamber.  The  chamber  voltage  was  measured  by  a  catiacitivelv  connled 
resistive  voltage  probe  of  division  ratio  1:20.000  into  a  50  U  load.  The  nsciime  of  this 
probe  after  proper  compensation  was  less  than  ~400  ps  and  its  RC  droop  time  constant 
is  ~  120  ns. 

A  movable  1  mil  graphite  Faraday  cup  of  3  cm  diam  was  inserted  on  axis  through  the 
end  flange  of  the  drift  tube  to  measure  the  electron  beam  current  at  certain  distances 
downstream  of  the  anode. 
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A  dual  beam  oscilloscope  (Tektronix  7844)  allowed  for  time  correlated  monitoring 
of  two  signals  simultaneously. 

The  emittance  of  the  electron  beam  was  analyzed  using  a  slit- hole  type  emittance 
meter.15  It  consisted  of  an  array  of  seven  200  /im  wide  and  2mm  long  slits  which  were 
cut  into  a  0.6  mm  thick  stainless  steel  plate.  These  slits  produce  sheet  beamlets  which 
were  detected  by  a  radiachromic  film16  placed  12  mm  further  downstream.  The  film 
material  gives  linear  coloration  response  up  to  an  absorbed  dose  of  10®  rad.  By  use 
of  an  applied  magnetic  field,  it  was  assured  that  UV  or  other  light  produced  in  the 
experiment  did  not  contribute  to  the  film  response.  The  resulting  spatial  distribution 
pattern  on  the  film  was  scanned  by  a  microdensitometer,  and  the  emittance  and  beam 
profile  were  analyzed.15 

The  breakdown  characteristic  of  the  discharge  chamber  was  determined  by  fixing  the 
voltage  to  a  certain  value  and  slowly  increasing  the  pressure  until  breakdown  occurred. 
In  Fig.  2  the  breakdown  voltage  Uq  is  plotted  as  a  function  of  pressure. 


Fig.  2.  Measured  characteristic  breakdown  curve. 

When  the  power  supply  with  voltage  Uo  is  connected  via  a  resistor  R.  the  chamber 
capacitance  C  is  charged  up.  At  a  rerrvn  pressure  and  voltage,  i'fj,  defined  by  the 
characteristic  breakdown  curve  (Fig.  2),  the  voltage  breaks  down.  The  charging  resis- 
*or  R  virtually  decouples  the  discharge  circuit  from  the  power  supply.  Therefore  the 
low  im;  -'•dance  discharge  cannot  be  sustained,  and  is  extinguished.  Consequently  the 
chamber  capacitance  C  is  changed  again  untii  the  breakdown  voltage  t’s  is  reached 
and  the  cycle  repeats.  Thus  a  periodic,  pulsed  discharge  is  obtained,  the  time  period 
r  of  which  is  approximately  determined  by 

f/0 

r=  RCln— — .  (1) 

to  -  Ob 
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By  choosing  the*  RC  value  adequately,  periods  of  typically  0.1  to  1  second  can  be 
obtained. 


RESULTS 

Figure  3  shows  a  typical  set  of  time  correlated  electrical  signals  obtained  by  Ro- 
gowski  coils,  voltage  probe  and  Faraday  cup.  The  discharge  current  is  measured  by  the 
upstream  coil  and  shows  a  damped  oscillating  waveform.  The  voltage  drops  within  20 
ns  from  90%  to  10%  of  its  maximum.  The  maximum  current  corresponds  to  almost 
zero  voltage  indicating  an  inductance  dominated  discharge.  The  downstream  coil,  as 
well  as  the  Faraday  cup,  register  only  one  negative  peak  at  the  time  of  the  discharge 
current  maximum. 


Fig.  3.  Typical  electrical  signals  of  discharge:  time  scale  10  ns/div.  a)  discharge 
current  (upstream  coil)  300  A/div;  b)  voltage  signal  10  kV/div: 
c)  injected  electron  beam  current  (down-stream  coil)  84  A/div; 
d)  Faraday  cup  signal  (15  cm  behind  anode)  20  A/div. 

All  signals  show  that  the  discharge  process  occurs  in  two  stages.  Up  to  time  tv  the 
voltage  drops  slowly,  and  currents  measured  by  upstream  and  downstream  coils  start 
smoothly.  At  time  tv  the  main  breakdown  occurs.  This  two-step  breakdown  is  also 
observed  in  the  waveform  of  the  downstream  Faraday  cup  signal.  Of  100  A  electron 
current  leaving  the  anode,  as  registered  by  the  downstream  Rogowski  coil.  10  A  are 
measured  by  the  Faraday  cup.  13  cm  downstream  of  the  anode. 

The  measured  discharge  currents  with  two  different  external  capacitors  were  fitted 
ro  an  expression 


no  =  /(0)exp  (  — 


R(t  -  In) 
2  L 


sinivO  t  -  I,-, 


i  2) 


where  ^  = 


ib-Ui)" 


1/2 


Thus,  the  least-squares  fitted  time  averaged  values  for  the  discharge  inductance  L 
and  the  resistance  II  were  determined,  as  summarized  in  fable  1. 
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Table  I.  The  least-squares  fitted  values  of 
inductance  L  and  resistance  R  for  different 
discharge  conditions. 


G(  pF) 

20  kV 

25  kV 

.10  kV 

Rif!)  .  LfnH'l 

R(fi) 

LfnH)  | 

R(fl) 

Linill 

420 

2.2  56.2 

2.1 

57.5  | 

1.7 

55.3 

760 

1.7  40.9 

1.9 

40.2  i 

;  i.s 

44.4 

The  dependence  of  discharge  current  and  beam  current  measured  by  the  upstream 
and  downstream  coil  as  a  function  of  charging  voltage  for  three  capacitance  values  is 
illustrated  in  Fig.  4.  The  current  signals  show  an  increase  in  amplitude  with  increasing 
voltage  L’b  and  capacitance.  The  beam  profile  and  emittance  were  measured  with  an 
external  capacitance  of  390  pF  operated  at  a  breakdown  voltage  of  24  kV.  These  results 
were  obtained  bv  averaging  more  than  20  shots. 


(IV) 


Fig.  4.  Upstream  and  downstream  coil  current  as  a  function 
of  charging  voltage  for  discharge  capacitance  values. 

To  find  the  beam  prolile  at  different  axial  distances,  the  experimental  data  was 
fitted  to  a  Gaussian  distribution.  The  values  for  the  root-mean-square  i  r ms  i  radius  6 
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are  shown  in  Fig.  5.  The  beam  expands  from  an  rms  radius  of  1.5  mm  to  2.5  mm 
measured  at  3  cm  and  18  cm,  respectively,  as  it  propagates  through  the  low  pressure 
gas. 


Fig.  5.  RMS  radius  6  as  function  of  distance  from  anode. 


Fig.  ti.  Density  distribution  (optical  density  O.D.)  found 
behind  a  7  slit  system  at  a  9  cm  distance  from  anode. 

Figure  b  shows  the  profile  of  the  scanned  beamlet  distribution  measured  at  9  cm 
downstream  of  the  anode.  The  rms  emittance,  frTnj,  is  defined  as  '' 


<rrr,J  =  4  (<  I2  ><  x'  2  >  ~  <  lx'  >;)  , 
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where  x'  or  y'  denote  the  gradients  of  the  particle  trajectories,  given  bv  x'  —  dx Id:  and 
/  =  dy/dz.  The  brackets  (<>)  indicate  averaged  values  using  the  four-dimensional 
density  distribution  p^(x.x'.  y,  y')  in  the  four-dimensional  “trace  space".  As  the  quan¬ 
tities  of  interest  in  brackets  (<>!  are  second  moments  of  x  -  x'  coordinates  only, 
the  more  convenient  projected  density  pi  on  two-dimensional  trace  space  can  be  used, 
where  p?(x,x')  =  p^(x,x\  y.y'Ydydy1  .lSl  Assuming  a  Maxwellian  transverse  velocity  dis¬ 
tribution.  a  Gaussian  density  profile,  and  an  axisymmetric.  nonrotational  beam,  the 
rms  emittance  for  the  electron  beam  is  calculated  by  numerical  integration15  and  found 
to  be 


(  =  55  mm  —  mrad 


-8  -4  0  4  8 

x(mm) 

Fig.  7.  Two-dimensional  emittance  contour  plot  of  trace  space 
density  pj(x.x').  The  lines  represent  projected,  constant  emittance 
values,  from  outside  to  inside:  99%.  S0%,  (>0%,  40%.  20%. 

For  an  estimated  mean  beam  energy  of  20  keV,  this  leads  to  a  normalized  rms  emittance 
of  about  15  mtn-mrad.  The  two-dimensional  emittance  plot  is  obtained  as  shown 
in  Fig.  7.  Assuming  a  Maxwellian  veiocitv  distribution  with  an  rms  velocity  e,-,,  = 
(kT / m0)l/2  in  the  transverse  direction,  the  normalized  emittance  at  the  source  can  be 
written  as 

-  2r,\kT/m0c:)i,:  i  .) 

where  r,  is  the  source  radius  and  T  is  the  temperature.  Assuming  the  same  value  of  »n 
=  15  m  in- in  rad  at  the  source,  and  the  center  hole  radius  l  (i  mm  as  the  source  radius, 
an  approximate  temperature  of  11  eV  is  obtained. 
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CONCLUSIONS 

The  electrical  discharge  characteristics  and  the  emittance  of  the  ejected  electron 
beam  of  a  pseudospark  discharge  have  been  investigated. 

The  pseudospark  breakdown  occurs  rather  reproduciblv  at  the  pressure  and  voltage 
defined  by  the  characteristic  breakdown  curve.  The  discharge  takes  place  in  two  stages, 
both  of  which  are  accompanied  by  electron  beam  emission.  The  “lectron  energy  dis¬ 
tribution  has  still  to  be  determined.  After  the  electron  emission  occurs,  the  discharge 
drops  into  a  low  impedance  mode  whose  waveform  was  modeled,  assuming  a  simple 
LRC  circuit  analysis.  The  electron  beam  is  self-focused  as  it  propagates  through  low 
pressure  gas.  The  time  integrated  emittance  'vas  measured  at  an  axial  distance  of  9 
cm.  An  rms  value  of  about  bo  mm-mrad  was  found.  Together  with  the  measured  beam 
current  /  of  bO  A  at  that  distance  and  the  normalized  rms  emittance  tn  —  lb  mm-mrad. 
the  normalized  brightness  B ,  is  obtained  as 

Bn  =  — s=  2  x  10“  ;l-r 
(<„)*  m-rau- 

Tiius  the  pseudospark  disciiarge  appears  to  represent  a  higii  brightness  electron  beam 
source. 
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High-brightness  pseuctosparit-produced  electron  beam 
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The  first  time-integrated  root  mean  square  (nns)  emittance  measurement  of  a  pseudospark  - 
produced  electron  beam  is  presented.  From  a  six-gap  pseudospark  chamber  with  argon 
working  gas,  ~  10  Hz  repetitive  pulsed  electron  beams  of  average  energy  —20  keV,  peak 
current  —  SO  A,  and  pulse  duration  — 10  ns  are  extracted  into  a  drift  tube.  A  typical  value  of 
measured  rms  emittance  is  found  to  be  £— 55  mm  mrad,  yielding  a  normalized  rms  emittance 
of  e„  =:  1 5  mm  mrad.  The  normalized  brightness  of  the  beam  is  then  estimated  to  be 
B„  z;4x  10'°  A/(m2  rad2). 


Interest  in  high-quality,  high-current  electron  beams 
has  been  stimulated  by  stringent  source  requirements1  for 
advanced  accelerators  such  as  e+e~  linear  colliders2  and  for 
high-power  free -electron  lasers. 3  The  pseudospark  discharge 
phenomenon  with  an  interesting  charged  particle  emission 
characteristic  was  reported  by  Christiansen  and  Schulth- 
eiss.4  The  pseudospark  discharge  may  be  characterized  as  a 
gas  discharge  between  a  planar  anode  and  a  hollow  cathode 
in  the  lower  pressure  side  of  the  characteristic  breakdown 
curve,  which  is  very  similar  to  the  Paschen  curve  for  parallel 
electrodes.  During  the  breakdown  phase,  emission  of  highly 
pinched  electron  beams  with  current  densities  of  up  to 
10*  A/cm2  were  observed.4  These  observations  suggested  to 
us  to  investigate  the  quality  and  brightness  of  the  beam  pro¬ 
duced  in  the  pseudospark.  In  this  letter,  for  the  first  time,  the 
measurement  of  time-integrated,  root  mean  square  (rms) 
emittance  of  the  electron  beam  produced  by  a  pseudospark 
discharge  is  described.  The  normalized  beam  brightness  is 
evaluated  from  the  measured  rms  emittance  and  the  average 
beam  energy  which  is  inferred  from  the  voltage  waveform. 

The  experimental  setup  is  schematically  shown  in  Fig. 
1 .  A  six-gap  pseudospark  chamber,  which  is  characteristi¬ 
cally  similar  to  the  devices  reported  by  other  laboratories,4-* 
is  used  in  this  work.  An  easy  to  change,  o-ring  sealed,  modu¬ 
lar  device  consists  of  a  hollow  cathode,  five  modules  of  inter¬ 
mediate  electrodes,  and  Plexiglas  insulating  disks,  and  an 
anode.  All  the  electrodes  are  made  of  brass  and  have  an  outer 
diameter  of  6.35  cm  and  a  center  hole  of  3.2  mm  diameter 
which  provide  a  channel  for  the  discharge  and  also  for  the 
electron  beam  passage.  The  3.2-mm-thick  Plexiglas  insula¬ 
tor  disks  have  outer  and  inner  diameters  of  7  and  2.54  cm, 
respectively.  An  external  capacitor  of  390  pF  is  added  to  the 
system  which  provides  an  adequate  amount  of  stored  energy 
for  the  electron  beam  production.  A  6.3-cm-diam,  50-cm- 
long  drift  chamber  is  attached  to  the  anode  side  of  the  dis¬ 
charge  chamber  to  accommodate  diagnostics  such  as  the 
Faraday  cup  and  the  emittance  meter.  A  capacitance-ma¬ 
nometer-type  vacuum  gauge  is  used  to  measure  the  gas  fill 
pressure  which  is  almost  statically  balanced  by  a  needle 
valve  and  a  throttle  vacuum  valve.  Two  miniature  Rogowski 
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coils  of  rise  time  —0.5  ns  are  molded  into  axisymmetric 
grooves  milled  into  both  sides  of  the  anode  flange  so  that  the 
azimuthal  component  of  the  magnetic  field,  which  arises 
from  the  axisymmetric  current,  is  predominantly  supported 
and  other  components  (noises)  are  suppressed.  A  movable 
Faraday  cup  consisting  of  a  100  mfl  current  viewing  resis¬ 
tor,  and  a  3-cm-diam  graphite  charge  collector  (which  is  not 
shown  in  Fig.  1 )  connected  to  the  end  of  a  rigid  coaxial  cable 
is  placed  on  axis  of  the  downstream  drift  chamber  to  mea¬ 
sure  the  electron  beam  current  at  various  axial  positions. 
This  system  is  also  conveniently  utilized  to  support  a  simple 
emittance  meter  which  is  attached  to  the  end  of  the  cable 
(see  Fig.  1 )  in  lieu  of  the  graphite  charge  collector;  this  al¬ 
lows  not  only  emittance  measurements  at  various  axial  posi¬ 
tions,  but  also  simultaneous  measurement  of  tne  beam  cur¬ 
rent.  The  discharge  chamber  voltage  is  measured  by  a 
homemade  high-impedance  resistive  voltage  probe  of  divi¬ 
sion  ratio  1:20  000.  The  rise  time  is  found  to  be  <0.5  ns,  and 
the  RC  droop  time  constant  due  to  the  blocking  capacitor  is 
— 120  ns.  The  rms  emittance  of  the  electron  beam  is  mea¬ 
sured  by  using  a  simple  slit-hole  emittance-meter  system.9  It 
consists  of  an  array  of  seven  parallel  slits  of  200  /im  width 
and  2  mm  spacing,  constructed  from  a  0.6-mm-thick  stain¬ 
less-steel  plate,  and  a  50-ftm  thick-  (2  mil)  radiachromic 
film, 10  as  a  beam  detector,  placed  1 2  mm  downstream  of  the 
slit  plane.  The  slits  sample  the  beam  at  given  slit  locations 
producing  flat  beamlets.  The  beamlets  are  then  allowed  to 
disperse  with  their  transverse  velocity  after  passing  through 
the  slits.  The  spatial  current  density  profiles  of  dispersed 
beamlets  at  the  detector  plane  are  recorded  on  the  radiachro¬ 
mic  film,  from  which  the  transverse  velocity  distribution  of 
the  beam  can  be  analyzed. 

The  system  is  operated  typically  at  30  kV,  with  an  argon 
gas  fill  pressure  of  40  mTorr.  With  a  200  Mfl  charging  resis¬ 
tor,  slow  repetition  frequencies  ( 1-10  Hz)  of  the  electron 
beam  pulses  are  easily  attained;  such  frequencies  allowed  us 
easier  control  of  the  number  of  beam  pulse  exposures  to  the 
emittance  meter.  The  radiachromic  film  of  the  emittance 
meter,  placed  9  cm  downstream  of  the  anode  plane,  is  ex¬ 
posed  to  approximately  five  consecutive  beam  pulses  to  pro¬ 
duce  an  adequate  density  profile,  with  the  peak  value  of  the 
optica]  density  well  below  0.5.  This  ensures  that  the  mea¬ 
sured  optica]  density  distribution  is  linear  to  the  beam  inten¬ 
sity.10  The  beam  current  measured  at  this  location  is  50  A. 
The  beam  current  density  profiles  of  all  beamlets  recorded 
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FIG.  I.  Schematic  representation  of  the  experimental  setup. 
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FIG.  3.  Eminence  plot  in  x-x'  space.  The  uodensity  contours  correspond  to 
0. 1.  20, 40, 60,  and  80%  of  the  peak  value. 


on  a  radiachromic  film  are  scanned  with  an  optical  micro- 
densitometer;  a  typical  scan  is  shown  in  Fig.  2. 

Therms  eminence 11  or  effective  emittance1'12  which  has 
been  widely  used  as  a  measure  of  beam  quality  is  defined  as 

f  =  4(<*2><x'2>  -  <xx'>2)1/2,  (1) 

where  x'  is  the  gradient  of  the  particle  trajectory  given  by 
x'  =  dx/dz  =  px/px,  and  the  angular  brackets  denote  aver¬ 
age  values  over  the  two-dimensional  trace  space  as 
(q> )  =  /  ~ 1  / p6(xpc')dx  dx',  where p  is  the  projected  cur¬ 
rent  density  in  two-dimensional  trace  space,  and  the  total 
current  is  given  by  I  —  $  p  dx  dx". 

It  is  very  reasonable  to  assume  that,  as  in  Ref.  9,  the 
beam  produced  in  this  experiment  is  axisymmetric  and  of 
Maxwellian  transverse  velocity  distribution.  Such  an  as¬ 
sumption  allows  us  to  use  the  simple  slit-hole-type  emittance 
meter  mentioned  above,  whose  results  can  be  easily  ana¬ 
lyzed.9  We  find  empirical  functions  a(r),0(r),  and  a(r)  as 
functions  of  radial  position  r  from  the  density  profile  shown 
in  Fig.  2,  where  a(r)  is  the  mean  diverging  angle,  0(r)  rep¬ 
resents  the  peak  values,  and  a(r)  is  the  rms  width  of  the 
individual  distributions.  Numerical  integrations  are  then 
performed  using  the  empirical  functions  a,  0,  and  a  to  find 
p(xpc')  and  the  moments  (x2),  (x'2),  and  (xx')  (see  Ref.  9 
for  details).  Several  iaodensity  contours  of  the  resultant 
p(xpc')  are  constructed  in  x-x'  space  (known  as  an  emit¬ 
tance  plot )  as  shown  in  Fig.  3.  The  measured  rms  emittance, 
given  by  Eq.  ( 1 )  with  the  obtained  moments,  is  found  to  be 
e=:55  mm  mrad.  The  error  analysis  on  the  emittance  is 
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FIG.  2.  Typical  optical  density  profile  of  radiachromic  film  produced  by 
electron  beamier*  is  the  emittance  meter. 

1747  Appl.  Ptiys.  LatL,  Voi.  56,  No.  16. 30  April  1900 


somewhat  difficult;  nevertheless,  a  rough  estimate  is  made 
and  found  to  be  ±  15%.  It  is  interesting  to  note  that  the 
emittance  plot  in  Fig.  3  shows  a  slightly  inverse  S-shaped 
contour,  this  is  indicative  of  the  existence  of  nonlinear  focus¬ 
ing  (defocusing)  in  the  beam  possibly  due  to  the  space- 
charge  field  of  substantial  strength.  The  normalized  rms 
emittance  defined  as  e„  =  0ye,  where  0  =  v/c  and 
Y  =  (1  —  01)  "  1/2 ,  is  invariant1 1-19  in  the  presence  of  axial 
acceleration,  and  is  conveniently  used  for  comparing  beam 
qualities  at  different  energies.  With  an  average  beam  energy 
of  20  keV  inferred  from  the  voltage  probe  signal,  the  normal¬ 
ized  emittance  is  estimated  to  be  f,  ~  15  mm  mrad.  Another 
••■sefu'.  invariant  of  the  beam  associated  with  the  emittance  is 
the  brightness,13  for  which  unfortunately,  no  universally 
agreed  exact  definition  exists.  We  use  here  a  typical  expres¬ 
sion'*  of  normalized  rms  brightness  Bm  =  27/(7^ ).  With  a  I 
measured  current  of  50  A,  we  have  obtained  the  brightness 
B„  ~4xl010  A/(m2rad2),  which  is  substantially  higher  I 
than  that  of  other  high  brightness  sources.319  It  should  be  j 
noted  here  that  the  electron  beam  pulses  are  very  reproduc¬ 
ible,  and  the  results  in  this  work  are  all  time  integrated  over  a 
few  shots  of  10  ns  pulses. 

In  conclusion,  we  have  measured,  for  the  first  time,  the 
emittance  of  a  pseudospark-produced  electron  beam.  A  six- 
gap  pseudospark  chamber  with  argon  gas  is  operated  at  30 
kV  and  produced  a  repetitive  pulse  train  of  electron  beams. 
Typically,  a  beam  of  average  energy  20  keV,  peak  current 
~  50  A,  and  pulse  duration  — 10  ns  full  width  at  half  maxi¬ 
mum  is  extracted  from  the  chamber.  The  rms  emittance  is 
evaluated  by  using  the  slit-hole  method.  The  typical  value  of 
the  rms  emittance  is  found  to  be  f~  55  mm  mrad,  yielding  a 
normalized  emittance  of  e„  ~  1 5  mm  mrad.  These  values  are 
time  integrated  over  a  few  shots  of  entire  pulses  and  repre¬ 
sent  the  upper  limit.  The  normalized  rms  brightness  of  the 
beam  is  then  estimated  as  B„  =  11  )  ~4x  1010  A J 

( m2  rad2 ),  which  is  substantially  higher  than  that  of  existing 
high-brightness  electron  beam  sources. 

We  are  grateful  to  M.  Reiser  for  valuable  discussions, 
and  K.  K.  Jain  for  critical  comments  in  preparing  the  final 
manuscript.  The  technical  assistance  of  D.  Cohen  and  C. 
Daniels  is  also  acknowledged.  This  work  is  supported  by  the 
Air  Force  Office  of  Scientific  Research  and  the  U.  S.  Depart¬ 
ment  of  Energy. 
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An  anode  foil  changer  for  pulsed,  high-current,  relativistic  electron  beam  diode  systems  is 
described.  This  device  allows  multiple  shots  to  be  taken  with  a  foil-anode  diode  system  while 
maintaining  a  continuous  vacuum.  Shot-to-shot  reproducibility  is  greatly  improved,  while  at  the 
^>ame  time  a  significant  time  savings  is  realized.  Also,  various  diagnostics  and  options  are  added  to 
the  basic  design  which  further  increase  its  utility  and  flexibility. 


Diode  systems  used  in  pulsed  intense  relativistic  electron 
beam  experiments  have  many  uses  in  current  research.'  and 
their  design  may  be  considered  to  be  centered  on  two  basic 
experimental  configurations2:  those  using  foil  anodes.3  and 
foilless  diodes.3-5  These  systems  have  various  design  needs 
which  must  be  addressed.  Foil  anodes  in  present  use  make 
for  a  small  anode-cathode  gap,  providing  the  low  impedance 
necessary  for  the  production  of  high-current  beams.  When 
the  beam  is  being  injected  into  a  gas-filled  chamber,  some 
type  of  barrier  must  be  used  to  separate  the  vacuum  diode 
area  from  the  gas-filled  area.  Diagnostics  are  often  included 
to  measure  various  beam  and  diode  characteristics,  such  as 
diode  current  or  injected  beam  current.  All  such  factors 
must  be  considered,  while  at  the  same  time  designing  the 
general  experimental  setup  to  be  easy  and  efficient  to  use. 

Unfortunately,  there  are  a  great  many  problems  that 
generally  go  along  with  such  foil  diode  systems.  Diode  sys¬ 
tems  using  a  foil  anode  require  that  the  diode  be  taken 
apart — and  the  vacuum  broken — after  each  shot  to  replace 
the  anode  foil.  Not  only  is  this  a  time-consuming  process,  but 
shot-to-shot  reproducibility  is  reduced  as  the  experimental 
setup  is  tampered  with  after  each  firing;  vacuum  and  mainte¬ 
nance  problems  may  also  be  aggravated  in  this  process. 

A  device  has  been  built  which  satisfies  all  these  require¬ 
ments.  w  hile  at  the  same  time  eliminating  many  of  the  prob¬ 
lems.  This  device — which  we  call  a  foil  changer — greatly 
enhances  experimental  reproducibility,  efficiency,  and  flexi¬ 
bility.  Figure  1  shows  the  experimental  diode  setup  with  the 
foil  changer  in  place.  This  foil  changer  allows  multiple  shots 
to  be  taken  with  a  foil-anode  diode  system  while  maintaining 
a  continuous  vacuum  and  without  disturbing  experimental 
conditions:  with  our  present  version,  over  20  shots  can  be 
taken  before  foil  replacement  is  necessary,  which  is  usually 
enough  to  complete  one  experimental  senes.  This  gives  the 
obvious  benefit  of  speeding  up  research  and  data  aquisition. 
Just  as  important,  the  integnty  of  the  expenmentai  setup  is 


assured  for  the  entire  senes  of  shots,  and  reproducibility 
from  shot  to  shot  is  improved. 

Figure  2  shows  a  diagram  of  the  foil  changer  device. 
This  particular  version  was  designed  to  be  compatible  with  a 
1 5-cm  vacuum  system  and  drift  chamber.  The  outer  diame¬ 
ter  is  16.8  cm.  Vacuum  flanges  compatible  with  Dependex 
vacuum  equipment  are  built  in  on  both  sides  of  the  foil 
changer.  A  4.4-cm  diam  hole  is  bored  through  the  center  of 
the  device  over  which  the  anode  foil  will  be  placed;  the  beam 
will  be  injected  through  the  foil  and  center  hole  into  the  drift 
region  beyond.  Two  7-cmX2.2-cmX  1.75-cm  windows  are 
milled  into  the  body  of  the  foil  changer  on  either  side  of  the 
beam  hole.  In  these  cavities  a  strip  of  anode  foil  material  is 
rolled  up  and  placed  inside;  the  stnp  of  foil  is  wrapped 
around  rods  which  extend  to  knobs  outside  the  device  via 
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Fig.  1.  Schematic  of  diode  system  with  foil  changer  in  place. 
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Fig.  2.  Diagram  of  foil  changer. 


vacuum  feedthroughs.  This  foil  strip  can  then  be  advanced 
across  the  anode  plane  from  one  window  to  the  other  after 
each  shot,  maintaining  vacuum  conditions  at  all  times.  Fig¬ 
ure  3  shows  a  photograph  of  the  foil  changer,  along  with  a 
sample  stnp  of  anode  foil  showing  the  senes  of  holes  created 
by  the  beam  dunng  an  expenmental  senes. 

The  anode  foil  itself  is  used  as  a  barrier  between  the 
diode  and  drift  regions.  Such  a  barrier  is  needed  in  experi¬ 
ments  in  which  the  drift  region  is  filled  with  a  gas,  while  the 
diode  area  remains  evacuated,  for  example  beam  injection 
into  a  neutral  gas.  To  achieve  this  barrier,  an  O  ring  is  placed 
in  a  groove  underneath  the  foil  strip  and  around  the  center 
hole.  A  cover  plate  with  a  central  window  of  the  same  diame¬ 
ter  as  that  of  the  center  hole  is  then  placed  over  the  foil  on  the 
diode  side.  This  cover  plate  may  be  opened  and  closed  by 
means  of  a  cam.  which  extends  to  a  lever  outside  the  body  of 
the  foil  changer  via  a  vacuum  feedthrough.  When  the  cover 
plate  is  "closed.”  the  foil  is  firmly  pressed  against  the  O  ring, 


Fig  .'  Photograph  of  foil  changer  and  sample  strip  of  anode  foil  after  an 
experimental  senes  o!  beam  shots  i  the  cover  plate  is  not  shown  t 


creating  a  hermetic  barrier  between  the  diode  and  drift  re¬ 
gions.  When  "opened"  after  tiring  a  shot,  the  pressure 
against  the  foil  is  released,  and  the  foil  stnp  can  be  advanced 
for  the  next  shot.  This  can  all  be  done  while  maintaining  a 
continuous  vacuum.  It  should  be  mentioned  that  this  cover 
plate  serves  another  useful  purpose;  it  protects  the  body  of 
the  foil  changer  from  beam  damage.  While  the  cover  plate  is 
a  separate  unit  and  is  easily  replaced,  the  foil  changer  itself  is 
not. 

The  basic  concept  of  this  foil  changer  has  been  given 
added  utility  and  flexibility  by  the  addition  of  vanous  design 
options.  For  diagnostics,  two  Rogowski  coils  have  been  built 
in  to  the  device.  The  coil  on  the  upstream  (diode)  side  mea¬ 
sures  the  diode  current,  while  the  coil  on  the  downstream 
(drift  region)  side  measures  net  injected  beam  current. 
These  coils  were  designed  using  a  toroidal  core  of  plexiglas 
wrapped  with  19  turns  of  3.2-mm  wide  copper  tape;  the  core 
diameter  is  1 3.6  cm  and  the  cross-sectional  area  is  2  mm  X  2 
mm.  These  coils  were  permanently  set  into  a  groove  cut  just 
inside  the  vacuum  flanges  by  means  of  transparent  epoxy. 
The  leads  of  the  Rogowski  coils  were  coupled  to  semirigid 
coaxial  cables  outside  the  device  through  small  holes  bored 
into  the  rim  of  the  foil  changer.  The  diode  and  injected  cur¬ 
rent  signals  can  then  be  monitored  by  oscilloscopes.  These 
Rogowski  coils  are  of  low  inductance,  with  a  rise  time  of 
L  /R  <  1  ns. 

An  added  feature  built  into  our  foil  changer  is  the  addi- 
riori  of  a  »)"T-valve  feedthrough  entrance  leading  to  just 
downstiearn  of  the  anode  foil  plane.  This  allows  the  foil 
changer  and  its  associated  diagnostics  to  be  used  for  foilless 
diode  configurations  also,  using  a  puff  valve  in  place  of  a  foil 
anode.  When  needed,  the  puff  valve  is  simply  fit  into  the 
vacuum  feedthrough  entrance  and  the  foil  stnp  is  removed; 
when  not  in  use.  the  entrance  is  simply  blocked  off. 

Taken  together,  the  capabilities  of  this  device  make  for 
much  more  reliable  and  efficient  expenmental  work.  Be¬ 
cause  an  entire  expenmental  senes  can  be  performed  in 
v  acuum,  without  taking  the  diode  apan  after  each  shot,  re¬ 
producibility  is  greatly  enhanced.  The  monolithic  and  com¬ 
pact  size  make  for  easy  expenmental  assembly,  while  conse¬ 
quently  minimizing  mechanical  problems.  The  many 
built-in  features  make  it  simple  to  change  expenmental  con¬ 
ditions.  i.e.,  foil  anode,  putf-gas  system,  etc.  Finally,  because 
it  is  no  longer  necessary  to  take  apan  the  diode  after  each 
shot  to  replace  the  foil,  a  significant  amount  of  time  is  saved; 
instead  of  taking  perhaps  one  series  of  ten  shots  in  a  day,  one 
can  take  several  series  of  twenty  shots  or  more. 

This  work  was  supported  by  the  Air  Force  Office  of 
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It  is  shown  that  for  a  beam  of  arbitrary  distribution  function,  the  root-mean-square  emittance 
is  invariant  under  any  linear  transformation  if  the  determinant  of  transfer  matrix  is  unity.  The 
in  -nance  condition  is  generalized  by  defining  a  normalized  root-mean-square  emittance  to 
include  the  case  with  a  nonunity  determinant  which  is  associated  with  an  axial  acceleration. 


Since  the  root-mean-square  ( rms)  emittance  was  intro¬ 
duced  by  Lapostolle'  and  Sacherer,2  it  has  been  widely  used 
as  a  figure  of  merit  for  beam  quality.  Lawson  et  al.1  have 
established  theoretical  grounds  for  use  of  the  rms  emittance 
in  which  the  entropy  of  the  beam  is  closely  connected  with  its 
rms  emittance.  The  nonlinear  effect  in  the  system  increases 
the  entropy  and  thus  the  rms  emittance.  In  addition,  a  num¬ 
ber  of  authors  have  shown4-6  that  the  rms  emittance  is  an 
invariant  of  the  motion  if  the  focusing  effect  is  linear.  The 
need,  for  reference  purposes,  for  a  general  description  of  in¬ 
variance  properties  became  apparent. 

In  this  Brief  Communication,  we  present  a  general  treat¬ 
ment  of  invariance  properties  of  the  rms  emittance  in  any 
linear  system.  In  the  linear  system,  the  evolution  of  the  dis¬ 
tribution  function  is  governed  by  a  linear  transformation  of 
the  individual  particle  point.  We  show  that  only  if  the  deter¬ 
minant  of  the  transfer  matrix  of  the  linear  transformation  is 
unity,  the  rms  emittance  is  invariant  regardless  of  the  form  of 
the  distribution  function.  This  invariance  condition  is  ex¬ 
tended  by  defining  a  normalized  rms  emittance  to  include 
the  nonunity  determinant  case,  which  is  associated  with  an 
axial  acceleration. 

The  x-plane  rms  emittance  is  defined1 2  as 

ex  =  k((x2)(x'2)  -  Ixx')2)1'2 ,  (1) 

where  x  is  the  gradient  of  particle  trajectory  given  by 
x  =  dx/dz  —  px/p,.  and  the  angular  brackets  ( (  ) )  denote 
an  average  value  over  the  two-dimensional  trace  space6  as 
{g)  =  igp(xjc‘  )dx  dx  .  The  validity  of  the  emittance  calcu¬ 
lation  with  p(xjt)  in  two-dimensional  trace  space  is  not 
simple  in  general  and  is  discussed  in  Refs.  4  and  5.  For  the 
purpose  of  the  present  work,  we  assume  a  simple  beam  with  a 
single  degree  of  freedom,  such  as  a  strip  beam  or  equivalent, 
in  which  the  x  and  y  motions  are  decoupled  so  that  the  beam 
can  be  represented  by  density  distribution  in  two-dimension¬ 
al  x-px  phase  space.  Furthermore,  we  assume  that  the  instan¬ 
taneous  axial  velocity  at  a  given  point  z  is  constant  over  the 
beam  cross  section.  The  p(xjc)  may  be  considered  as  the 
instantaneous  density  distribution  function  of  particle  tracks 
left  on  the  x-x  trace  plane,  at  a  given  value  of  z,  by  a  short 
slice  of  beam  as  it  passes  by.  The  scale  factor  k,  commonly 
taken  as  4,  will  be  omitted  throughout.  It  should  be  pointed 
out  that  the  rms  emittance  is  not  a  mathematical  rms  value  of 
emittance  as  opposed  to  what  the  term  as  well  as  the  symbols 
commonly  used  seem  to  imply. 

We  consider  a  linear  transformation7  in  two-dimension¬ 
al  trace  space. 


x 

x 


a 

c 


(2) 


by  which  values  of  x  and  at  point  z  are  related  to  those  at 
z  | .  The  elements  of  the  transfer  matrices  a,b,c ,  and  d  may  be 
functions  of  z  but  independent  of  x,  and  x'.  This 

linear  transformation  has  properties  that  straight  lines 
transform  into  straight  lines  and  ellipses  transform  into  el¬ 
lipses.  Particle  tracks  within  a  square-shaped  differential 
area  dxt  dx\  at  z,  transform  into  a  parallelogram  dx  dx’  at  z 
as  depicted  in  Fig.  1 .  Since  the  particles  are  conserved  in  this 
linear  transformation,  the  same  number  of  tracks  must  be 
found  within  dxt  dx\  and  dx  dx' .  This  is  expressed  by  the 
relation 


p,(x,jt\  )dxy  dx\  -  p(xX)dx  dx' 

=  p(axt  -f  bx\  ,cx ,  +  dx\  )\J  | <£x,  dx\  , 

(3) 


where  J  =  dUjc'l/dCx,,*; )  =  ad  -  be.  a  Jacobian,  which 
in  this  case  is  equal  to  the  determinant  of  transfer  matrix  |  A  \ . 

Now  we  consider  a  transformation  of  the  rms  emittance 
value  from  z ,  to  z.  Each  term  in  Eq.  ( 1 )  is  calculated  in  terms 
of  variables  x,  and  x\  with  the  aid  of  an  integral  theorem 

f  F(xry)dx  dv  =  I  F  [xt  u.v).y(u.v)  j \J  \du  dv  , 


FIG.  1  The  linear  transformation  is  schematically  illustrated.  The  transfer 
mains  elements  used  in  this  example  are  a  =1.6  =  0.5.  c  =0.5,  and  d  =  2 
and  the  corresponding  determinant  and  Jacobian  are  \A  I  =  J  =  1.75. 
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as 

<x:)  =  J  x:p(xjc')dx  dx' 

=  J  (ax,  +  bx\  ); 

xp(ajc,  +  bx\ ,cx,  -t-  dx\  )\J  \dx,  dx\  .  (4) 

From  the  particle  conservation,  Eq.  (3 ),p  in  Eq.  (4)  is  re¬ 
placed  by  p, (*,.*; ),  yielding 

(x:)  =  J  (ax,  -rbx‘,)2p,(xlX,)dx,dx't 

=  a:<xj )  •+•  lab  <x,x; )  -t-  b  2(x'2)  .  (5) 

Similarly, 

<x'2)  =c:(x])  •+■  led  <X|jc; )  -+-  d  '(x'1)  .  (6) 

and 

(xx  )  =  ac(x2,  >  -f-  (ad  +  bc)(x,x\ )  -+-  bd  (x':)  .  (7) 

Combining  Eqs.  ( 5)— (7),  we  find  the  following  simple 
expression  for  the  transformation  of  rms  emittance: 

(x2)  (x‘2)  -  (xx')2 

-  (ad  -  bc)2((x2,)(x',2)  -  U,*;  );)  ,  (8) 

f,  =  A  f. 

Equation  ( 8 )  clearly  shows  that  the  rms  emittance  is  invar¬ 
iant  under  any  linear  transformation  if  and  only  if  the  trans¬ 
fer  matrix  determinant  is  unity.  Typical  examples  of  such  a 
transformation  are  force-free  drift  and  thin  lens  systems.  It  is 
noted  that  this  result,  Eq.  (8),  is  a  consequence  of  the  linear 
transformation  imposed  but  with  an  arbitrary  distribution 
function p(xjc  ) .  It  should  not  be  construed  that  Liouville's 
theorem  is  not  applicable  to  this  case.  The  trace-space 
area  is  (mn0yc)~'  times  that  in  phase  space,  since 
x  =  p%/ (mn0yc).  Therefore,  the  transformed  area  in  x 
-  x  space  is  subject  to  change  if  the  factor  0y  changes  in  the 
transformation. 

It  is  of  interest  that  the  determinant  A  can  be  related  to 
'he  paraxial  ray  equation/1  The  paraxial  equation  may  be 
reduced  into  a  general  form 

x"  -  g(z)x'  -  h(z)x  =  0  .  (9) 

Since  Eq.  (9)  is  linear  and  second  order,  it  therefore  has  a 
general  solution  which  is  any  combination  of  two  indepen¬ 
dent  solutions.  Supposing  that  Eq.  ( 9 )  governs  the  motion  of 
particles  in  a  linear  system  described  by  the  transformation 
,  Eq.  i  2  ;  J.  the  transformation  x  =  ax,  —  bx‘.  may  be  con¬ 
sidered  to  be  a  general  solution  to  Eq.  (9)  where  .t,  and  x\ 
are  initial  conditions  at  z,.  The  terms  a(z)  and  b(z),  being 
independent  solutions,  must  satisfy  Eq.  (9)  separately, 

j  '  -  ga  —  ha  =  0  .  b  "  —  gb  '  -+-  hb  =  0  .  <  10) 

These  simultaneous  equations  are  solved  for  g(z),  yielding 

g(z)  =  -  W/W=  -  A  17  A  ;  .  (It) 

where  W^is  the  Wronskian.  and  the  prune  denotes  the  deriva¬ 
tive  with  respect  to  z.  If  x  and  x'  in  Eq.  ( 2 )  are  continuous, 
then  a'  =  c  and  b'  =  d,  and  thus  the  Wronskian  is  equal  to 


the  determinant  of  the  transfer  matrix:  W  =  A  |.  Hence,  one 
finds  the  determinant  in  terms  ofg(z)  as 


(12) 


In  the  absence  of  a  magnetic  field,  the  term  g(z).  the  coeffi¬ 
cient  of  x  of  paraxial  equations,  is  reduced'  to 

g(z)  =  y'/(P:y)  ,  ( 13) 


which  is  associated  with  an  axial  acceleration.  Equation 
(13)  is  then  substituted  into  Eq.  ( 12).  and  one  finds  that 


A  =(3sy,/(Py)  . 


(14) 


It  is  evident  that  after  substitution  of  Eq.  (14)  into  Eq.  ( 8 ) ,  a 
new  emittance  defined  as 


=0yet  =0,yxet  (15) 

is  in  variant.  This  may  be  referred  to  as  normalized  rms  emit¬ 
tance  resemblant  to  the  normalized  emittance.  It  is  clear 
from  Eq.  ( 14)  that  since  the  Jacobian  is  equal  to  | A  I,  the 
transformed  rms  emittance,  which  has  a  dimension  of  area, 
is  inversely  proportional  to  0y.  This  is  caused  by  the  fact 
that,4  s  as  in  the  case  with  normalized  emit:ance,  if  an  area  in 
x-px  phase  space  is  invariant  by  Liouville’s  theorem  and p,  is 
changed  by  the  presence  of  axial  acceleration  during  the 
transformation,  the  area  in  x-x‘  trace  space  is  reduced  by  the 
factor  (0y)  ~ '.  For  the  same  reason  *s  normalized  emit¬ 
tance,'  use  of  the  normalized  rms  emittar.ee  may  be  conven¬ 
ient  when  considering  beams  in  which  there  is  acceleration, 
and  also  when  comparing  beams  at  different  energies  and 
with  different  particles. 

In  conclusion,  it  is  shown  that  the  rms  emittance  is,  in 
general,  an  invariant  quantity  under  any  linear  transforma¬ 
tion  if  the  determinant  of  transfer  matnx  is  unity.  By  defin¬ 
ing  normalized  rms  emittance.  the  invariance  condition  is 
generalized  to  include  the  case  with  a  nonunity  transfer  ma¬ 
tnx  determinant  when  axial  acceleration  exists. 
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Root-mean-square  emittance  increase  in  a  nonlinear  lens 
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An  exact  derivation  of  the  root-mean-square  emittance  increase  when  a  beam  of  arbitrary 
distribution  passes  through  a  nonlinear  lens  system  is  presented.  It  is  found  that  the  emittance 
increase  depends  on  the  spherical  aberration  coefficient  of  the  lens  and  characteristics  of  the 
beam  distribution. 


The  root-mean-square  (rms)  emittance  growth  asso¬ 
ciated  with  the  nonlinearity  of  a  focusing  lens  system  has 
been  an  important  issue, 1  as  the  spherical  aberration  is  inher¬ 
ently  associated  with  the  focusing  lenses  used  in  accelerating 
and  transport  channels.  It  is  well  known  that  the  nonlinear¬ 
ity  of  an  array  of  lenses  gives  rise  to  the  filamentation  of 
emittance,  resulting  in  an  increase  of  the  rms  emittance 
while  the  phase-space  area  occupied  by  the  beam  remains 
unchanged.  Given  the  nonlinear  nature  of  the  phenomenon, 
the  studies  on  this  subject  have  often  been  undertaken  by 
heuristic  computational  simulations.2  We  would  like  to  pres¬ 
ent  here  a  general  treatment  of  rms  emittance  increase  when 
a  beam  of  arbitrary  distribution  passes  through  an  ideal  thin 
lens  with  spherical  aberration. 

The  rms  emittance  is  defined3  4  as 

(  =  k((x:)(x2)  -  (xx')2)'12.  (1) 

where  k  ( the  scale  factor  commonly  taken  as  4 )  will  be  omit¬ 
ted  throughout,  x  is  the  gradient  of  the  particle  trajectory 
given  by  x  -  dx/dz  =  px/p:.  and  the  angular  bracket  de¬ 
notes  an  average  value  over  the  two-dimensional  trace  space 
as 

(g(xjc')>  =  N  ~  '  J  g(xpc')p(xjt')dx  dx',  (2) 

where p(x,;c')  =  J”  Xpt(x,y,x'  ,y‘)dy  dy  is  the  projection  of 
the  four-dimensional  distribution  pjU.y.x'.v')  onto  two-di¬ 
mensional  trace  space  and  N  =  J'  *  .  p  dx  dx' . 

We  consider  a  beam  of  arbitrary  distribution  p,( x,.x\  ), 
which  is  passed  through  a  one-dimensional  thin  lens  with 
spherical  aberration.  Then,  the  beam  is  transformed  as 

x  =  x„  (3) 

x  —  x't  +  ax,  +  bx\,  (4) 

where  coefficients  a  and  b  may  be  expressed  in  terms  of  the 
focal  length  /  and  the  spherical  aberration  coefficient  C.  of 
the  lens:  a  =  l//and6  =  C./f*.  It  should  be  noted5  here  that 
there  is  a  coupling  between  x  and  y  motion  in  the  axisymmet- 
ric  lens  with  the  spherical  aberration:  thus  the  transforma¬ 


tion  in  such  a  lens  cannot  be  described  by  Eq.  (4). 

As  shown  in  Ref.  6,  if  the  Jacobian  of  the  transformation 
is  unity,  the  averaged  value  given  by  Eq.  (2)  in  the  trans¬ 
formed  space  coordinates  can  be  conveniently  found  in 
terms  of  the  x,-x;  coordinates  and  distnbutionp,(JC,pc; )  as 

(g(xpc'))  =  A"' J  g(x,x')p(xpc')dxdx' 

=  yV-lJ  g(x,X,  +  ax,  +  bx\ ) 

Xp,(x,j t\  )J  dx,  dx\,  (5) 

viheTcJ  =  d(x.x')/d(x,X,  )  =  1  for  the  transformation  giv¬ 
en  by  Eq.  ( 2 ).  Thus  Eq.  ( 5 )  can  be  reduced  to  a  simpler  form 

as 

(g(xpc'))  =  (g(x,,  x\  +  ax,  +  bx},  )>,,  (6) 

where  the  subscript  1  after  the  angular  bracket  [which  de¬ 
notes  integration  over  the  x,-x\  space  with  distnbution 
p,(x,pt\  )  ]  will  be  omitted  throughout.  Thus  the  moments 
necessary  for  the  present  analysis  can  be  readily  found  in  an 
algebraic  form  as 

(xmx  ")  =  (x'"(x\  -f  ax ,  -rbx',)n).  (7) 

Now  we  evaluate  the  moments  in  Eq.  ( 1 )  in  terms  of  x,- 
coordinates  with  the  aid  of  Eq.  (7)  as 

U:>  =  <xj),  (8) 

(x’:)  =  (u;  +  ax,  -r  bx\  ):) 

—  -t ~  a~  (x~, )  -t -  b  ~ {x , )  +  2a {x ,x , ) 

~  2b  <jc,'x;  )  +  lab  <x, ).  (9) 

{xx‘)  =  <X|(x;  4-  ax,  +  bx\ )) 

=  <x,x; )  -t-  a(x:,  )  -+-  b  <x4 ).  (10) 

Substituting  Eqs.  ( 8 )— (  10)  in  Eq.  (1),  we  find  the  trans¬ 
formed  rms  emittance  as 
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=  <xf><*;J>  -  (x,x\ 

+  62(<x2X*?>  -  <jct>2) 

+  2 6(  < X }  > (x]x\ )  -  (x\X\  )(x\)) 

=  cf  +  6  :a,  4-  2b0t.  (11) 

Equation  (11)  shows  that  the  square  of  emittance  is  in¬ 
creased  by  b  2a,  +  2 b0y,  where  a,  and  0,  are  measures  of 
some  characteristic  quantities  of  beam  distribution  before 
the  transformation.  It  is  noted  that  the  emittance  increase  in 
Eq.  ( 1 1 )  is  dependent  on  the  nonlinear  coefficient  b  in  Eq. 
(2),  but  independent  of  the  linear  focusing  coefficient  a  as 
expected,  since  the  rms  emittance  is  invariant6  under  any 
linear  transformation. 

It  is  very  interesting  to  show  how  the  quantities  a  and  0 
transform  as  the  beam  passes  through  the  nonlinear  lens. 
Again  using  Eq.  (3),  the  moments  and  product  moments 
composing  the  quantities  a  and  0  are  readily  found,  yielding 

a=(x-)(xb)  -(x*)2 
=  <x;><*?>  -  <*;>2 

=  a,  (12) 

and 

0  =  <x2){x}x')  -  (xx')<x4) 

=  (xj)<Jt{x;>  -  <x,x;X.r?> 

+  6(<x2)<x6>  -  (xTX) 

—  0X  +  ba,.  (13) 

Note  that  the  quantity  a  remains  unchanged,  but  the  change 
in  0  is  simply  the  nonlinear  coefficient  b  times  a,.  This  sim¬ 
ple  transformation  of  a  and  0  through  the  nonlinear  lens 
may  be  very  useful  in  the  evaluation  of  rms  emittance  growth 
in  successive  transformations  such  as  in  periodic  focusing 
channels. 

As  an  example,  a  beam  of  distribution 

pixjc')  =  exp[  -  x’/(2<52)  —  x'V(2 a2)  |  ( 14) 

is  passed  through  a  thin  nonlinear  lens  that  transforms  the 
beam  distribution  according  to  Eqs.  (3)  and  (4).  It  is  then 
straightforward  to  find  the  rms  emittance  e  and  the  quanti¬ 
ties  a  and  0  before  the  transformation  by  direct  integration 

[Eq.  (2)]: 

e;  —  I'cr, 

a,  =66".  (15) 

0,  =0. 

After  the  transformation  the  same  quantities  are  found  by 
using  Eq.  (6)  as 

e2  =  f)  +  b  la,  +  2b0,  =  52cr  +  66  ’<5H, 
a  =  66*,  (16) 

0  =  bat  =  6  b6*. 

The  same  results  of  this  particular  case  may  also  be  obtained 
by  direct  integration  using  Eq.  (2),  which  is  lengthy  but 
straightforward.  The  emittance  diagrams  in  x-x'  space  of  the 

a  S3 


FIG  !  Emittance  diagrams  in  x-x 
trace  space  for  (a)  a  distribution 
p(x„ jc',  )  =  exp(  -  *;/(26;)  -  x\'/ 
(2<r)|,  where  <5  =  0  01  m  and 
<7  =  0.01  rad  (b)  The  same  distnbu- 
tion  transformed  by  a  =  x „ 
x  =  x\  i-  ax,  bx),  where  a  =  -  . 
rad/m  and  b=  —  1000  rad/m' 


beam  before  and  after  the  traprf~rmation  in  this  example  are 
illustrated  in  Fig.  1 . 

In  conclusion,  we  have  presented  an  exact  derivation  of 
the  rms  emittance  increase  of  an  arbitrary  beam  in  a  one¬ 
dimensional  nonlinear  lens  system.  The  quantities  a  and  0 
associated  with  the  beam  distribution  transform  linearly 
through  the  nonlinear  lens  system.  This  fact  may  be  useful  in 
determining  the  emittance  growth  in  a  transport  or  acceler¬ 
ating  svxtem  in  which  successive  focusing  lenses  are  em¬ 
ployed.  The  results  may  be  useful  for  a  variety  of  one-dimen¬ 
sional  electrostatic  and  magnetic  lens  systems.  It  may  also  be 
possible  to  extend  the  use  of  these  results  in  order  to  estimate 
the  emittance  increase  in  other  types  of  systems  such  as  those 
using  quadrupoie  and  axisymmetric  lenses. 
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